| 查看: 3464 | 回复: 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朵）
» 猜你喜欢
0856调剂，是学校就去
已经有9人回复
-
298-一志愿中国农业大学-求调剂
已经有9人回复
-
265求调剂
已经有7人回复
-
求调剂
已经有4人回复
-
学校已经提交到NSFC，还能修改吗？
已经有7人回复
-
材料学硕297已过四六级求调剂推荐
已经有6人回复
-
295复试调剂
已经有7人回复
-
本人考085602 化学工程 专硕
已经有20人回复
-
一志愿中南化学337求调剂
已经有5人回复
-
281求调剂（0805）
已经有23人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 高金求助matlab解微分方程组
已经有12人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab 解矩阵方程 重金悬赏解决方案
已经有10人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【原创】多项式和非多项式曲线拟合的MATLAB实现视频教程
已经有10人回复
- 求助曲线拟合
已经有6人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
- 求助：用Origin做微分方程的曲线拟合
已经有7人回复
【答案】应助回帖
★ ★
余泽成(金币+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
smithjohn
金虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 1334.9
- 散金: 372
- 红花: 4
- 帖子: 673
- 在线: 148.8小时
- 虫号: 1030430
- 注册: 2010-05-28
2楼2011-06-03 23:28:10
3楼2011-06-04 13:36:21
4楼2011-06-04 14:39:09
6楼2011-06-17 19:15:06
7楼2011-06-17 19:16:02
送鲜花一朵 | 支持 |
8楼2011-06-21 10:10:08
9楼2011-06-21 11:32:29
liang201101
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 169.3
- 帖子: 120
- 在线: 25.2小时
- 虫号: 1444020
- 注册: 2011-10-15
- 性别: GG
- 专业: 制造系统与自动化
10楼2013-05-23 21:44:09
