色散方程用matlab编译,急切求助用matlab曲线拟合色散方程 - 程序语言 - 小木虫 - 学术 科研 互动社区...
或者你仿照这个代码改改
function y = fitme_diheral( c,phi )
% A1=c(1);
% k1=c(2);
% x1=c(3);
% kBT=443*8.3145/1000;
aa=c(1);
bb=c(2);
cc=c(3);
dd=c(4);
ee=c(5);
ff=c(6);
psi=(phi-pi);
% y=0.5*(f1*(1+cosx)+f2*(1-cos(2*x))+f3*(1+cos(3*x))+f4*(1-cos(4*x)));
y=aa-bb.*cos(psi)+cc.*(cos(psi)).^2-dd.*(cos(psi)).^3+ee.*(cos(psi)).^4-ff.*(cos(psi)).^5;
end
clear;
clc;
load diherals.txt;
data=diherals;
%%avge the theta<0 and theta>0
a=data(1:180,2);
b=data(361:-1:182,2);
avg=0.5*(a+b);
data(1:180,2)=avg;
data(361:-1:182,2)=avg;
p=[24.5156799887984,-1.90007004205358,-1.66217137679272,3.06265999012061,2.33333792446086,-2.48454162655902;];
% data(:,2)=acf(data(:,2));
data3=data;
data3(:,1)=data3(:,1)*pi/180;
for i=1:50
% [p,res]=lsqcurvefit('dihe',p,data3(:,1),data3(:,2));
[p,res]=lsqcurvefit('fitme_diheral',p,data3(:,1),data3(:,2));
%注意修改函数名dihe
plot(data3(:,1),data3(:,2),'ro');
hold on
plot(data3(:,1),fitme_diheral(p,data3(:,1)));
%注意修改函数名dihe(p,data3(:,1)
end,
色散方程用matlab编译,急切求助用matlab曲线拟合色散方程 - 程序语言 - 小木虫 - 学术 科研 互动社区...相关推荐
- matlab 摄氏度符号怎么打,【求助】弱弱的问下,Matlab如何在公式中输入希腊字母? - 程序语言 - 小木虫 - 学术 科研 互动社区...
很多时候都要在matlab画图的时候添加一些公式符号之类的,有一些特殊的字符并不能直接从键盘上输入,比如希腊字母等等.但是有想用,因为这样使图看起来漂亮而且容易理解. 例如:我想输入摄氏度的符号,怎么 ...
- matlab实现kuan滤波,用Matlab设计一个0.5-20HZ的带通滤波器! - 程序语言 - 小木虫 - 学术 科研 互动社区...
这是我在网上看到的,希望对你有所帮助. ----------------------------------------------------------------------------- fu ...
- matlab二重指针,VC++中函数返回数组指针或者带指针的结构体的编译方式是否可取? - 程序语言 - 小木虫 - 学术 科研 互动社区...
libralibra 构造函数+析构函数应该是比较正规的做法吧, 你的结构体/类被提早释放的问题,是不是那个结构体是局部变量,当超出作用域时对象被销毁,如果存在析构函数,会被默认调用, 如果你返回结构 ...
- matlab风力机叶片仿真教程,我用matlab编写了一个风力机叶片的优化设计程序,有点问题 - 程序语言 - 小木虫 - 学术 科研 互动社区...
主程序 clear; global r; global R; global phi; global lamb; V1=7;P=400;lamb0=6;B=3; D=2.5; R=D/2; N=60*l ...
- matlab中并行用不了,请教matlab在linux下的并行问题 - 程序语言 - 小木虫 - 学术 科研 互动社区...
刚在网上搜了一下,matlabpool open local 8 设置节点数,但是必须在matlab下面才行,比如我打开linux, >>matlab 进入matlab系统,然后matla ...
- matlab雷诺曲线,多重网格法求解雷诺方程的MATLAB编程 - 程序语言 - 小木虫 - 学术 科研 互动社区...
这是MATLAB有限差分法求解雷诺方程,可以参考下.希望对你有用 function [Fx,Fy]=FDM_circular(X,Y,X_dot,Y_dot) Fx=0; Fy=0; epsilon= ...
- c语言调用dll实例 结构体传参,Matlab调用Dll,对于自定义结构体数组如何传参? - 程序语言 - 小木虫 - 学术 科研 互动社区...
%[例子1] % Matlab调用DLL,整形数组传参 v = [1,22,33;44,-56,67]; pv = libpointer('int16Ptr',v); m = get(pv,'Valu ...
- matlab如何取等距点,matlab编程等距取点 - 程序语言 - 小木虫 - 学术 科研 互动社区...
各位大神,我想用matlab编一个程序,目的是对于一组数据等距取点,但不会写,请各位赐教! 首先,我有一个21×4的矩阵A[21,4],希望对每一列数据每间隔3个点取一个数,然后将这些取出来的点保存在 ...
- matlab imfill holes,到底怎么实现,bwfill和imfill - 程序语言 - 小木虫 - 学术 科研 互动社区...
认真研读一下MATLAB的help文档吧,解释最权威: BWFILL Fill background regions in binary image. BWFILL is a grandfathere ...
- matlab简单的程序,一段简单的matlab程序 - 程序语言 - 小木虫 - 学术 科研 互动社区...
原程序是可以运行的,为全面理解程序内容,我将分以下几个部分进行分析: 1."for x=varx"怎么理解? 请参看Matlab关于for函数的帮助文件: Syntax:for i ...
最新文章
- android 7.0 mi5s,官方回应 小米5s升级安卓7.0还要再等等
- python 项目中的requirements.txt文件
- 老股民经验之谈 这些股票买入必死无疑
- 数据库限制内存使用方法
- 特定场景下取代if-else和switch的方式
- 一张纸折多少次可以变成珠穆朗玛峰那么高?
- linux8安装bbr_CentOS 7安装4.9内核开启BBR
- 秦始皇:请叫我数据治理专家
- JDK8 Stream 操作
- 华为年底将推多款平板产品:顶配版将搭载5nm麒麟9000
- 4399小游戏flash插件怎么下载_Flash即将关闭,但这个小游戏平台,或许可以帮你找回4399的回忆...
- python是一种代表简单主义思想的语言_setuptools-36.6.0
- js ajax 传输list,jQuery ajax请求返回list数据动态生成input标签,并把list数据赋值到input标签...
- 公司打卡少几秒分析--学会拒绝,学会选择
- Android反编译工具dex2jar的使用
- 卡西欧函数计算机广告词,卡西欧函数计算器,终身学习如有神助!
- 威金蠕虫(网吧杀手)肆虐互联网 九千用户十余企业遭攻击
- ClickHouse遇见RoaringBitmap
- 【JAVA长虹键法】第 十二式 外观模式(23种设计模式)
- 北京PMP考试考点在哪里?