| 查看: 1377 | 回复: 2 | |||
[交流]
微分方程组参数拟合 已有2人参与
|
|
编一个关于微分方程组拟合的问题，能运行的了，也能出结果，但是得到的三个参数最佳值是错误的，且随着上下限取的不同而不同，能帮忙看看怎么办吗？ 这个是我科研立项最后也是最关键一步，但是Matlab用的不熟，希望有大神帮我解决一下问题，真心感谢。 具体的微分方程是 然后求 上面L和下面L不一样，在程序里有区分。最后是得到M矩阵里M、第一个L和tsd的最佳参数值，拟合数据是t和第二个L的值（程序中已经给出了）。 代码贴出来不知道怎么回事会乱码，所以就这么直接贴了 function HH clear all clc T=[502272 762325 1114354 1364122 1706100 1964980 2140916 2313059 2489034 2742752 3001044]; L=[5.65432E+42 1.30331E+43 9.00958E+42 5.3601E+42 3.56514E+42 2.5611E+42 ... 2.27497E+42 1.94825E+42 1.50971E+42 1.22059E+42 1.04568E+42]; x0=[0.1;0.1]; M0=[0,0,0]; %M0=[0.001*1.989*10^33,1,10^42]; lb=[0.1,0.1,0.1]; ub=[0.5*1.989*10^33,1728000,1*10^46]; Lexp=... [502272 5.65432E+42 762325 1.30331E+43 1114354 9.00958E+42 1364122 5.3601E+42 1706100 3.56514E+42 1964980 2.5611E+42 2140916 2.27497E+42 2313059 1.94825E+42 2489034 1.50971E+42 2742752 1.22059E+42 3001044 1.04568E+42]; %t=0.1:1:402800; [M,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@g,M0,lb,ub,[],T,x0,Lexp); A=M(1) B=M(2) C=M(3) end function LL=g(M,T,x0,Lexp) c=3*10^10; k=0.5; [t,Eint]=ode45(@f,T,x0,[],M); LL=Eint(:,2).*t./(3.*k.*M(:,1)).*(4.*pi.*c.*Eint(:,1))-Lexp(:,2); end function dEint=f(t,Eint,M) %目标函数,M(1)=M,M(2)=tsd,M(3)=L,Eint(1)代表v，Eint(2)代表Eint c=3*10^10; k=0.5; dEint(1)=Eint(2)./(M(1).*Eint(1).*t); dEint(2)=(M(3)./((1+t./M(2)).^2)-4.*pi.*c.*Eint(2).*t.*Eint(1)./(3.*k.*M(1))-Eint(2)./t); dEint=dEint(  ;end 麻烦哪位大神能帮忙认真看一下吗。这是我科研立项最重要一步，已经卡这里很久了，真心没什么办法了。 |
回复此楼» 猜你喜欢
070200求调剂，一志愿某211，288分
已经有18人回复
-
法国博士后职位
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有298人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有9人回复
2楼2016-06-07 23:19:02
3楼2019-05-03 10:11:27
10
