| 查看: 4216 | 回复: 10 | ||
| 本帖产生 1 个 程序强帖 ，点击这里进行查看 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
急切求助用matlab曲线拟合色散方程
|
||
|
现在有如下色散方程： y^2 = A + B/(x^2-C) - D x^2 ----------------Eq. 1 其中A，B, C, D为待定系数，现在有(x,y)的一组实验值； 如何得到最佳的A、B、C、D 我的思路如下： step1 通过四对(x_i,y_i) 得到 A，B，C，D 值， step2 代入Eq.1 输入x实验值，得到y拟合值， step3 判断y拟合值和y实验值差值， step4 如果差值大于10^-4则调节A,B,C,D; step5 用新的ABCD，返回到step2 现在遇到的问题是 1.不收敛 2. 调节过的A，B，C，D有时候使得Eq.1方程右边小于零 求助：要是用matlab如何处理， [ Last edited by yangger on 2011-6-3 at 12:53 ] |
回复此楼» 本帖已获得的红花（最新10朵）
liuweihong1980» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有23人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有22人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有14人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高金求助matlab解微分方程组
已经有12人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab 解矩阵方程 重金悬赏解决方案
已经有10人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【原创】多项式和非多项式曲线拟合的MATLAB实现视频教程
已经有10人回复
- 求助曲线拟合
已经有6人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
- 求助：用Origin做微分方程的曲线拟合
已经有7人回复
3楼2011-06-04 13:36:21
smithjohn
金虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1334.9
- 散金: 372
- 红花: 4
- 帖子: 673
- 在线: 148.8小时
- 虫号: 1030430
- 注册: 2010-05-28
2楼2011-06-03 23:28:10
4楼2011-06-04 14:39:09
【答案】应助回帖
★ ★
余泽成(金币+2, 程序强帖+1): 谢谢参与应助！ 2011-06-04 19:30:05
yangger(金币+3): 非常感谢~~@ 2011-06-17 19:02:33
余泽成(金币+2, 程序强帖+1): 谢谢参与应助！ 2011-06-04 19:30:05
yangger(金币+3): 非常感谢~~@ 2011-06-17 19:02:33
|
或者你仿照这个代码改改 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 |
5楼2011-06-04 14:39:52
