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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab求解方程中的参数
222/3返回列表上一页123下一页
查看: 7055  |  回复: 21
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

yupeifeng

金虫 (小有名气)
★ ★
fegg7502: 金币+2, 多谢交流 2012-07-09 08:08:39
引用回帖:
5楼: Originally posted by dingd at 2012-07-03 08:52:49
1stOpt求解：

Constant c0=0.7,cp=0.61;
Variable t,c;
ODEFunction c'=-k*(c0-c)^(1/3)*(c-cp)
Data;
t,c
0,0.69
2,0.645
4,0.635
8,0.62
24,0.61
48,0.61


均方差(RMSE): 0.00535058522113925
...

Parameter k;
Constant c0=0.7,cp=0.61;
Variable t,c;
ODEFunction c'=-k*(c0-c)^(1/3)*(c-cp);
Data;
t,c
0,0.69
2,0.645
4,0.635
8,0.62
24,0.61
48,0.61
我这样定义的，但运行后没有结果 出现，我哪里出现错误了？我把你写得复制过去，提示没有定义定义参数。不知道啥原因啊
赞一下 回复此楼
憧憬未来
11楼2012-07-03 10:39:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
版本老了。
回复此楼
12楼2012-07-03 11:01:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yupeifeng

金虫 (小有名气)
引用回帖:
10楼: Originally posted by dingd at 2012-07-03 10:38:35
你用的是网上的1.5PJ版吧？老了点，好像也不支持微分方程计算。用新版本没问题。

你可以帮我发个新版本吗？谢谢啦，麻烦你了664363288@qq.com
赞一下 回复此楼
憧憬未来
13楼2012-07-03 16:15:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yupeifeng

金虫 (小有名气)
引用回帖:
10楼: Originally posted by dingd at 2012-07-03 10:38:35
你用的是网上的1.5PJ版吧？老了点，好像也不支持微分方程计算。用新版本没问题。

我在网上找不到新版本的？麻烦你给我发下啊，
赞一下 回复此楼
憧憬未来
14楼2012-07-03 16:30:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fegg7502: 金币+3, 多谢交流 2012-07-09 08:08:51
与硬件关联，没法给你。

难得有个像样的国产软件，有能力该支持下啊。
赞一下 回复此楼
15楼2012-07-03 20:12:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yupeifeng

金虫 (小有名气)
引用回帖:
15楼: Originally posted by dingd at 2012-07-03 20:12:29
与硬件关联，没法给你。

难得有个像样的国产软件，有能力该支持下啊。

怎么解决这个问题呢
赞一下 回复此楼
憧憬未来
16楼2012-07-03 20:49:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yupeifeng

金虫 (小有名气)
引用回帖:
10楼: Originally posted by dingd at 2012-07-03 10:38:35
你用的是网上的1.5PJ版吧？老了点，好像也不支持微分方程计算。用新版本没问题。

请问，新版本去哪里搞呢？纠结啊。
赞一下 回复此楼
憧憬未来
17楼2012-07-04 08:53:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
用新版本的都是自己或单位买的正版。
赞一下 回复此楼
18楼2012-07-04 13:56:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yupeifeng

金虫 (小有名气)
引用回帖:
18楼: Originally posted by dingd at 2012-07-04 13:56:45
用新版本的都是自己或单位买的正版。

不知道该怎么称呼你，我有三组数据，分别想套用7个动力学方程，找出最匹配的那个（R最接近1），你可以 帮我用这个软件做一下吗？非常感谢！
三组数据分别为：（1）
（t,c)
0              6.94
1        6.39
2        6.25
4        6.1
12        5.95
22        5.8
48        5.79
(2)组：
0        6.94
1        6.45
2        6.35
4        6.2
12        6.1
24        5.99
48        5.97
第三组：
0        6.94
2        6.39
4        6.3
8        6.2
12        6.1
48        6.05

动力学方程分别为：
(1)  -dc/dt=k(c0-c)^(1/3)(c-c~);
(2) -dc/dt=k(c0-c)^(2/3)(c-c~)^p。p分别取1、2、3、4。也就是相当于有四个方程。
（3） -dc/dt=k(c0-c)^(4/3)(c-c~)^p.p分别取1、2、3、4。也就是相当于有四个方程。
(4) -dc/dt=k(c-c~)^n,其中n分别取1、2、3、4。也就是相当于有四个方程。
我想三组数据分别用软件拟合这些方程，由k的误差判断应该取哪一个。
赞一下 回复此楼
憧憬未来
19楼2012-07-04 15:12:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gangbs

新虫 (小有名气)
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fegg7502: 金币+3, 多谢交流 2012-07-09 08:09:09
两种思路，合适与否大家一起探讨
1.由微分方程分析得出带k的解析解，再通过最小二乘法拟合得到参数k;
2.由数据点，用一个3次或4次的多项式来拟合c=f(t), 再将得到的f(t)对t求导，得到c'=g(t),将以上数据点中的t代入，则得到数据点(t,c'),再对上面的微分方程进行拟合得到k值。
赞一下 回复此楼
20楼2012-07-08 16:44:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yupeifeng 的主题更新
222/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭