版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

quuq432

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 11
帖子: 13
在线: 9.5小时
虫号: 4110317
注册: 2015-09-29
专业: 恒星与星际物质

[交流] 微分方程组拟合问题已有1人参与

本人MATLAB新手，编一个关于微分方程组拟合的问题，怎么就是运行不了，能帮忙看看怎么办吗？
具体的微分方程是
$dv/dt=Eint/(M*v*t)$
$dEint/t=(L/(1+t/tsd)^2)-4*pi*c*Eint*t*v/(3*k*M)-Eint/t$
然后求得
$L=Eint*t*(4*pi*c*v)/(3*k*M)$
拟合的数据对象是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的值

回复此楼

» 猜你喜欢

1楼 2016-02-26 20:04:52

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

参考：
http://muchong.com/bbs/viewthread.php?tid=10036072&page=2&fpage=1

赞一下