| 查看: 3286 | 回复: 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 |
回复此楼» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有193人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
|
您好 您能帮我运行一下这个程序么 我用的1stopt没有注册 以前可以运行出来 现在类似的程序就运行不出来了。谢谢您呀 ，  或者您看我的程序错误的地方我自己改也可以 谢谢您呀Title "ir"; Parameters a,b,c,d,e,f,g,h,m; Variable x1,x2,x3,x4; Function y=a+b*x1+c*x1^2+d*x1^3+e^x2+f*x3+g*x3^2+h*x3^3+m^x4; date; x1 x2 x3 x4 y 0.25 1 18.1 138 29.4 0.25 2 9.47 330 10.3 0.25 4 3.39 138 32.1 0.25 8 2.01 330 7.3 1 1 9.47 138 34.9 1 2 18.1 330 30.1 1 4 2.01 138 24.4 1 8 3.39 330 10.4 4 1 3.39 330 30.1 4 2 2.01 138 32.4 4 4 18.1 330 32.1 4 8 9.47 138 27.2 8.75 1 2.01 330 42.8 8.75 2 3.39 138 64.7 8.75 4 9.47 330 40.9 8.75 8 18.1 138 84.7 |
14楼2016-03-17 16:37:13
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
投票: