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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab求解方程中的参数
222/3返回列表上一页123下一页
查看: 7164  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工 求调剂 +12 enenenhui 2026-03-13 13/650 2026-03-16 08:30 by Linda Hu
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 0703化学调剂，求各位老师收留 +7 秋有木北 2026-03-14 7/350 2026-03-15 17:30 by 小物理化学
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 求调剂 +3 清风问长安 2026-03-09 3/150 2026-03-14 02:15 by JourneyLucky
[考研] 288求调剂 +14 王晓阳- 2026-03-09 19/950 2026-03-14 02:05 by JourneyLucky
[考研] 332分材料工程调剂 +3 莓好时光海苔 2026-03-09 3/150 2026-03-14 02:03 by JourneyLucky
[考研] 333求调剂 +3 球球古力 2026-03-09 3/150 2026-03-14 01:57 by JourneyLucky
[考研] 环境调剂 +6 晓看天暮看云 2026-03-09 6/300 2026-03-14 01:16 by JourneyLucky
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] 材料与化工304求B区调剂 +5 邱gl 2026-03-11 6/300 2026-03-13 22:37 by JourneyLucky
[考研] 336求调剂 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 一志愿西南交大，材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
[考研] 一志愿：武汉理工，材料工程，英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
关闭
请填处理意见
关闭