| 查看: 3059 | 回复: 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 |
回复此楼» 猜你喜欢
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有254人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有20人回复
-
考博自荐
已经有1人回复
-
山东大学第二批博士研究生招生
已经有0人回复
16楼2017-10-08 13:53:46
2楼2016-03-03 19:41:21
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
3楼2016-03-03 21:12:00
4楼2016-03-03 21:50:59
5