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

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

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

Fly_LP

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 4
在线: 4.5小时
虫号: 3030737
注册: 2014-03-09

引用回帖:

13楼: Originally posted by dingd at 2016-03-06 15:37:16
光看效果的话应该还可以吧：

均方差(RMSE):1.10232446404849E42
残差平方和(SSE):1.45814306884775E85
相关系数(R): 0.957906895473774
相关系数之平方(R^2): 0.917585620396204
修正R平方(Adj. R^2): 0.899 ...

您好您能帮我运行一下这个程序么我用的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