| 查看: 2900 | 回复: 16 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
1stopt微分方程组参数拟合已有2人参与
|
||
|
真心希望有人能帮我一下，这个问题已经困扰我很多天了，这10个金币是我所有的金币了。 具体的内容是这样 然后再得 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 |
回复此楼» 猜你喜欢
求2026年在台湾举行的物理和材料领域国际学术会议信息
已经有0人回复
-
求国际会议网站
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有248人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有19人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
17楼2017-10-11 21:25:03
2楼2016-03-03 19:41:21
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
3楼2016-03-03 21:12:00
4楼2016-03-03 21:50:59