或者你仿照这个代码改改

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曲线拟合色散方程 - 程序语言 - 小木虫 - 学术 科研 互动社区...相关推荐

  1. matlab 摄氏度符号怎么打,【求助】弱弱的问下,Matlab如何在公式中输入希腊字母? - 程序语言 - 小木虫 - 学术 科研 互动社区...

    很多时候都要在matlab画图的时候添加一些公式符号之类的,有一些特殊的字符并不能直接从键盘上输入,比如希腊字母等等.但是有想用,因为这样使图看起来漂亮而且容易理解. 例如:我想输入摄氏度的符号,怎么 ...

  2. matlab实现kuan滤波,用Matlab设计一个0.5-20HZ的带通滤波器! - 程序语言 - 小木虫 - 学术 科研 互动社区...

    这是我在网上看到的,希望对你有所帮助. ----------------------------------------------------------------------------- fu ...

  3. matlab二重指针,VC++中函数返回数组指针或者带指针的结构体的编译方式是否可取? - 程序语言 - 小木虫 - 学术 科研 互动社区...

    libralibra 构造函数+析构函数应该是比较正规的做法吧, 你的结构体/类被提早释放的问题,是不是那个结构体是局部变量,当超出作用域时对象被销毁,如果存在析构函数,会被默认调用, 如果你返回结构 ...

  4. 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 ...

  5. matlab中并行用不了,请教matlab在linux下的并行问题 - 程序语言 - 小木虫 - 学术 科研 互动社区...

    刚在网上搜了一下,matlabpool open local 8 设置节点数,但是必须在matlab下面才行,比如我打开linux, >>matlab 进入matlab系统,然后matla ...

  6. matlab雷诺曲线,多重网格法求解雷诺方程的MATLAB编程 - 程序语言 - 小木虫 - 学术 科研 互动社区...

    这是MATLAB有限差分法求解雷诺方程,可以参考下.希望对你有用 function [Fx,Fy]=FDM_circular(X,Y,X_dot,Y_dot) Fx=0; Fy=0; epsilon= ...

  7. c语言调用dll实例 结构体传参,Matlab调用Dll,对于自定义结构体数组如何传参? - 程序语言 - 小木虫 - 学术 科研 互动社区...

    %[例子1] % Matlab调用DLL,整形数组传参 v = [1,22,33;44,-56,67]; pv = libpointer('int16Ptr',v); m = get(pv,'Valu ...

  8. matlab如何取等距点,matlab编程等距取点 - 程序语言 - 小木虫 - 学术 科研 互动社区...

    各位大神,我想用matlab编一个程序,目的是对于一组数据等距取点,但不会写,请各位赐教! 首先,我有一个21×4的矩阵A[21,4],希望对每一列数据每间隔3个点取一个数,然后将这些取出来的点保存在 ...

  9. matlab imfill holes,到底怎么实现,bwfill和imfill - 程序语言 - 小木虫 - 学术 科研 互动社区...

    认真研读一下MATLAB的help文档吧,解释最权威: BWFILL Fill background regions in binary image. BWFILL is a grandfathere ...

  10. matlab简单的程序,一段简单的matlab程序 - 程序语言 - 小木虫 - 学术 科研 互动社区...

    原程序是可以运行的,为全面理解程序内容,我将分以下几个部分进行分析: 1."for x=varx"怎么理解? 请参看Matlab关于for函数的帮助文件: Syntax:for i ...

最新文章

  1. android 7.0 mi5s,官方回应 小米5s升级安卓7.0还要再等等
  2. python 项目中的requirements.txt文件
  3. 老股民经验之谈 这些股票买入必死无疑
  4. 数据库限制内存使用方法
  5. 特定场景下取代if-else和switch的方式
  6. 一张纸折多少次可以变成珠穆朗玛峰那么高?
  7. linux8安装bbr_CentOS 7安装4.9内核开启BBR
  8. 秦始皇:请叫我数据治理专家
  9. JDK8 Stream 操作
  10. 华为年底将推多款平板产品:顶配版将搭载5nm麒麟9000
  11. 4399小游戏flash插件怎么下载_Flash即将关闭,但这个小游戏平台,或许可以帮你找回4399的回忆...
  12. python是一种代表简单主义思想的语言_setuptools-36.6.0
  13. js ajax 传输list,jQuery ajax请求返回list数据动态生成input标签,并把list数据赋值到input标签...
  14. 公司打卡少几秒分析--学会拒绝,学会选择
  15. Android反编译工具dex2jar的使用
  16. 卡西欧函数计算机广告词,卡西欧函数计算器,终身学习如有神助!
  17. 威金蠕虫(网吧杀手)肆虐互联网 九千用户十余企业遭攻击
  18. ClickHouse遇见RoaringBitmap
  19. 【JAVA长虹键法】第 十二式 外观模式(23种设计模式)
  20. 北京PMP考试考点在哪里?

热门文章

  1. Web浏览器测试,怎么提取测试点 - web测试方法总结
  2. C语言程序——math库的函数应用
  3. android 编译 libjpeg-turbo,编译Android环境的libjpeg-turbo
  4. vlfeat python
  5. QT网络编程——UDP
  6. 使用HBuilderX将H5网页打包成APP
  7. glide 4.0 之后 和圆角裁剪CenterCrop冲突
  8. ffmpeg之图片串流h264视频
  9. 反编译工具Reflector下载
  10. 正则 负数 python_如何系统的自学 Python?“从新手到大师”的百天之路(上)