24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

quuq432

新虫 (初入文坛)

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

[求助] 1stopt微分方程组参数拟合已有2人参与

真心希望有人能帮我一下，这个问题已经困扰我很多天了，这10个金币是我所有的金币了。
具体的内容是这样 $dv/dt=Eint/(M*v*t)$
$dEint/t=(L/(1+t/tsd)^2)-4*pi*c*Eint*t*v/(3*k*M)-Eint/t$
然后再得
$LL=Eint*t*(4*pi*c*v)/(3*k*M)$
c是光速为3*10^10，k=0.5，M、tsd和L是未知参数，需要拟合求解，现在只知道能得到t和LL的实验值。具体的计算应该是v和Eint都是不知道的，若知道M、tsd和L的值，我用的是ode45，参数v和Eint根据t来回算（初始值是0,0），然后带入到LL里面。现在反过来了不知道M、tsd和L的值我就不知道怎么写了，用MATLAB的ode45和lsqnonlin写了个程序，最后求出来不是很理想，而且跟初始值那些影响蛮大。用1stopt写了一个，运行不了（我觉得肯定是错的），不知道怎么写，能帮我改改，然后算出来吗？真的麻烦了，在线等啊。。。
Parameters M,tsd,L;
Constant c=3*10^10,k=0.5;
Variable t,v,E,LL;
Function v'=E/(M*v*t);
         E'=(L/((1+t/tsd)^2))-4*pi*c*E*t*v/(3*k*M)-E/t;
         LL=E*t*(4*pi*c*v)/(3*k*M);
Data;
//t          LL
502272       5.65432*10^42
762325       1.30331*10^43
1114354       9.00958*10^42
1364122       5.3601*10^42
1706100       3.56514*10^42
1964980       2.5611*10^42
2140916       2.27497*10^42
2313059       1.94825*10^42
2489034       1.50971*10^42
2742752       1.22059*10^42
3001044       1.04568*10^42

回复此楼

» 猜你喜欢

求助VISSIM破解版软件已经有0人回复
求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有110人回复
求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有19人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助！！！1stopt进行常微分方程组参数模拟已经有5人回复
1stOpt软件微分方程参数拟合，运行后没有任何反应。已经有5人回复
1stopt 微分方程数据拟合求助已经有2人回复
1stopt微分方程组的参数估计的请教已经有3人回复
用1stOpt求解微分方程已经有2人回复
[求助] 1stopt 微分方程代码求助已经有5人回复
求帮忙运行一下下面的1stopt的程序，微分方程-参数估计已经有6人回复
用1stopt软件做微分方程组的参数估计已经有13人回复
求助，matlab/ 1stOpt 自定义微分方程拟合实验数据来求方程中的参数已经有9人回复
求助1stopt或者matlab 微分方程拟合已经有0人回复
求助，matlab/ 1stOpt 自定义微分方程拟合实验数据来求方程中的参数已经有1人回复
求助，matlab/ 1stOpt 自定义微分方程拟合实验数据来求方程中的参数已经有4人回复
帮忙看看能否用1stopt计算出微分方程组已经有0人回复
1stOpt求解常微分方程边值问题已经有5人回复

1楼2016-03-03 17:11:35

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

quuq432

新虫 (初入文坛)

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

引用回帖:

5楼: Originally posted by dingd at 2016-03-04 11:37:12
建议将一般公式：
LL=E*t*(4*pi*c*v)/(3*k*M);

写成微分方程形式：LL'=。。。。

这样才能进行微分方程计算

Parameters M,tsd,L;
Constant c=3*10^10,k=0.5;
InitialODEValue t=0.1,E=0.1,v=0.1,LL=0;
Variable t,LL；
ODEFunction v'=E/(M*v*t);
         E'=(L/((1+t/tsd)^2))-4*pi*c*E*t*v/(3*k*M)-E/t;
         LL'=E*(4*pi*c*v)/(3*k*M)
Data;
//t          LL
0.1          0.1
502272       1.12575E+37
762325          1.70965E+37
1114354          8.08503E+36
1364122          3.92934E+36
1706100          2.08964E+36
1964980       1.30337E+36
2140916          1.06261E+36
2313059          8.42281E+35
2489034          6.06545E+35
2742752          4.45025E+35
3001044          3.4844E+35

改成这样了，还是不能运行