| 查看: 663 | 回复: 1 | |||
[交流]
微分方程组拟合问题已有1人参与
|
|
本人MATLAB新手，编一个关于微分方程组拟合的问题，怎么就是运行不了，能帮忙看看怎么办吗？ 具体的微分方程是 然后求得 拟合的数据对象是t与L 其中t和L的数据是 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 写的程序比较乱，也运行不了，哪位大神能看看吗？ 我的渣渣程序是 function H clear all clc T=xlsread('数据.xlsx','sheet1','A1:A11'); L=xlsread('数据.xlsx','sheet1','B1:B11'); x0=[1 1]; t=0.1:1:402800; M=lsqcurvefit(@g,x0,T,L) end function LL=g(M,t,Eint) c=3*10^10; k=0.5; [t,Eint]=ode45(@f,t,[1,1]); size(Eint(:,1)); size(Eint(:,2)); LL=Eint(:,2).*t./(3.*k.*M(:,1)).*(4.*pi.*c.*Eint(:,1)); end function dEint=f(t,Eint,M); %目标函数,M(1)=M,M(2)=tsd,M(3)=L 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 想拟合出曲线，并且知道微分方程M tsd L的值 |
回复此楼» 猜你喜欢
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2016-03-06 21:45:54