应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 隐式的微分方程组求解问题。
81/1返回列表
查看: 1978  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

cyl19881108

铁虫 (初入文坛)

[求助] 隐式的微分方程组求解问题。

方程组如下：
dx/dt=a*x*(1-x/b)*(1/(1+s/c));
dp/dt=d*dx/dt+ex;
ds/dt=f*dx/dt+g*dp/dt+h*x;
其中a-h都是未知参数，解这个方程目的是得到x(t)、p（t）、s（t），进而使用实验数据拟合这3个方程并得到a-h参数的具体数值。
我就是想知道怎么通过上面的微分方程组求得x(t)、p（t）、s（t）这三个方程啊？
哪位大神可以教教我啊。编程这块我确实是个小白，要是哪里写的不清楚各位大神可以问我啊。
使用什么软件的什么程序能把这个微分方程组解出我想要的这三个方程呢？求各位大神指教啊。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-09-04 21:56:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
就是微分方程拟合，有数据吗？放上来看看。
赞一下 回复此楼
2楼2013-09-05 08:57:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyl19881108

铁虫 (初入文坛)
引用回帖:
2楼: Originally posted by dingd at 2013-09-05 08:57:06
就是微分方程拟合，有数据吗？放上来看看。

t[0,6,12,18,24,30,36,42,48,54,60]
s[24,22,19.3,16.0,13.4,9.9,6.5,3.7,1.92,1.3,0.96]
x[0.12,0.27,1.321,1.264,1.22,1.33,1.321,1.245,1.264,1.151,1.151]
p[0,1.02,1.98,3.37,5.05,6.66,8.65,10.95,14.00,15.50,16.40.
赞一下 回复此楼
3楼2013-09-06 13:47:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
★ ★
csgt0: 金币+2, 谢谢 2013-09-06 15:41:28
1stOpt求解代码，x项的效果较差：
CODE:
ConstStr f1=a*x*(1-x/b)*(1/(1+s/c)),
         f2=d*f1+e*x,
         f3=f*f1+g*f2+h*x;
Variable t,x,p,s;
ODEFunction x'=f1;
            p'=f2;
            s'=f3;
Data;
t=[0,6,12,18,24,30,36,42,48,54,60];
x=[0.12,0.27,1.321,1.264,1.22,1.33,1.321,1.245,1.264,1.151,1.151];
p=[0,1.02,1.98,3.37,5.05,6.66,8.65,10.95,14.00,15.50,16.40];
s=[24,22,19.3,16.0,13.4,9.9,6.5,3.7,1.92,1.3,0.96];

均方差(RMSE):0.499058750578204
残差平方和(SSE):7.47178909586035
相关系数(R): 0.999326893708328
相关系数之平方(R^2): 0.998654240488736
确定系数(DC): 0.998768314212006
F统计(F-Statistic): 0.0522044372000502

参数                  最佳估算
--------------------        -------------
a        0.0446360761798174
b        -4203.33960099553
c        18491.7942131467
d        -14454.704777378
e        645.620163822298
f        738.371259987294
g        -2.60499188266666
h        -32.3637736613547
隐式的微分方程组求解问题。
c1.jpg


隐式的微分方程组求解问题。-1
c2.jpg


隐式的微分方程组求解问题。-2
c3.jpg
赞一下(2人) 回复此楼
4楼2013-09-06 14:45:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyl19881108

铁虫 (初入文坛)
引用回帖:
4楼: Originally posted by dingd at 2013-09-06 14:45:43
1stOpt求解代码，x项的效果较差：

ConstStr f1=a*x*(1-x/b)*(1/(1+s/c)),
         f2=d*f1+e*x,
         f3=f*f1+g*f2+h*x;
Variable t,x,p,s;
ODEFunction x'=f1;
            p'=f2;
            s'=f ...

你用的是1stOpt 哪个版本的啊 ？ 第一个拟合的好像不好啊。。
赞一下 回复此楼
5楼2013-09-06 20:21:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyl19881108

铁虫 (初入文坛)
引用回帖:
4楼: Originally posted by dingd at 2013-09-06 14:45:43
1stOpt求解代码，x项的效果较差：

ConstStr f1=a*x*(1-x/b)*(1/(1+s/c)),
         f2=d*f1+e*x,
         f3=f*f1+g*f2+h*x;
Variable t,x,p,s;
ODEFunction x'=f1;
            p'=f2;
            s'=f ...

我只有istOpt1.5哦，好像不能带着微分拟合啊。。所以想把微分方程组先解出来。
赞一下 回复此楼
6楼2013-09-06 20:26:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

怀疑你x的数据是否正确（看看图）。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
7楼2013-09-06 22:11:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyl19881108

铁虫 (初入文坛)
引用回帖:
7楼: Originally posted by dingd at 2013-09-06 22:11:08
怀疑你x的数据是否正确（看看图）。

X数据没错啊，我源数据的图也是那样的
赞一下 回复此楼
8楼2013-09-06 22:26:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cyl19881108 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭