24小时热门版块排行榜

返回列表

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

chengxl

木虫 (小有名气)

应助: 5 (幼儿园)
金币: 1520.7
红花: 2
帖子: 247
在线: 64.4小时
虫号: 594608
注册: 2008-09-06
专业: 化工系统工程

[求助] 1stopt 语句修改已有4人参与

y'=k*y
已知不同初始条件下（t=0）初始值及其对应的t=10时函数值，如何用1stopt 语句拟合k值。下列语句运行不了，请高手指导，感激不尽。假定函数不能积分得到解析解。
Parameters k;
Variables a,y;
Function y=int(k*y,t=0,10)+a;
RowData;
0 1 2 3;
1 3 4 5;

回复此楼

» 猜你喜欢

散金币祈福已经有187人回复
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有190人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
澳门大学IAPME李海峰教授常年招收博士研究生已经有0人回复
深圳大学应用技术学院招聘凝聚态物理博士后已经有114人回复
深圳大学2027级光电信息工程/物理学研究生推免已经有70人回复
求助中医院临床中医中标的青年标书学习，谢谢已经有4人回复
浙大物院表面物理课题组长期招收普博及直博学生已经有257人回复
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry 已经有39人回复

1楼 2014-12-05 21:01:21

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

feixiaolin

荣誉版主 (文坛精英)

专家经验: +518
应助: 942 (博后)
贵宾: 1.275
金币: 3901
散金: 58785
红花: 532
沙发: 11
帖子: 24219
在线: 2601.8小时
虫号: 2139575
注册: 2012-11-21
专业: 光学信息获取与处理
管辖: 数学

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
chengxl: 金币+2 2014-12-09 08:33:43

引用回帖:

11楼: Originally posted by xxz903 at 2014-12-08 09:15:10
你这个可以直接用微分方程拟合啊！不过1.5版本的应该不能运行！代码如下：
Parameters k;
Variables a,y;
odeFunction y'=k*y;
RowData;
0, 1, 2, 3;
1, 3, 4, 5;

Parameters c, k;
Variable t, y;
Function  y=c*(exp(k*t)-1)+1;
data;
0  1
1  3
2  4
3  5
// ----------------------------------------
"Type your title here"

迭代数: 66
计算用时(时:分:秒:微秒): 00:00:00:94
优化算法: 粒子群算法PSO
计算结束原因: 达到收敛判断标准
均方差(RMSE): 0.0895640117970001
残差平方和(SSE): 0.0320868488366927
相关系数(R): 0.998194643831721
相关系数之平方(R^2): 0.996392546974336
决定系数(DC): 0.996332931561521
卡方系数(Chi-Square): 0.00435149314601671
F统计(F-Statistic): 552.407773509935

参数最佳估算
---------- -------------
c       -5.45183726787717
k       -0.428004961121383
====== 结果输出 =====
No 实测值y 计算值y
1 1 1.0000000
2 3 2.8982851
3 4 4.1356029
4 5 4.9420968