汕头大学海洋科学、生物学、生物与医药等3个专业接受调剂
查看: 6753  |  回复: 21
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yupeifeng

金虫 (小有名气)

[交流] matlab求解方程中的参数已有4人参与

搞了这么久还是没有解决这个问题。现在求助于高手给予指点。谢谢了!!
实验数据为:(t,c)=(0,0.69)(2,0.645)(4,0.635)(8,0.62)(24,0.61)(48,0.61).其中t为时间,c为某离子的浓度。
动力学方程模型为:-dc/dt=k*(c0-c)^(1/3)*(c-c~).
其中c0为初始浓度可以取0.7,c~为平衡浓度取0.61.
怎么样才能拟合出参数k的值呢?谢谢大家给出程序代码,再次感谢
回复此楼
憧憬未来
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
fegg7502: 金币+2, 多谢交流 2012-07-09 08:07:55
1stOpt求解:
CODE:
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
残差平方和(SSE): 0.000143143811043369
相关系数(R): 0.955733024504782
相关系数之平方(R^2): 0.913425614129058
决定系数(DC): 0.852429060780033

参数                  最佳估算
--------------------        -------------
k        1.02159845408008
5楼2012-07-03 08:52:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 22 个回答

dbb627

荣誉版主 (著名写手)

★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
fegg7502: 金币+4, 多谢交流 2012-07-09 08:07:15
这个其实论坛里有很多例子,参考就能写出来。
给你写了个
CODE:
function parafit
clear all;
t=[0 2 4 8 24 48];
y=[0.69 0.645 0.635 0.62 0.61 0.61];
y0=0.69;

% Nonlinear least square estimate using lsqnonlin()
beta0=0.5;
lb=[0];ub=[inf];
[beta,resnorm,residual,exitflag,output,lambda,jacobian] = ...
    lsqnonlin(@Func,beta0,lb,ub,[],t,y,y0);         
ci = nlparci(beta,residual,jacobian);
beta;
% result
fprintf('\n Estimated Parameters by Lsqnonlin():\n')
fprintf('\t k = %.4f ± %.4f\n',beta(1),ci(1,2)-beta(1))
fprintf('  The sum of the residual squares is: %.1e\n\n',sum(residual.^2))

% plot of fit results
tspan = [0  max(t)];
[tt yc] = ode45(@ModelEqs,tspan,y0,[],beta);
tc=linspace(0,max(t),200);
yca = spline(tt,yc,tc);
plot(t,y,'ro',tc,yca,'r-');
hold on
xlabel('Time');
ylabel('Concentration');
hold off
% =======================================
function f1 = Func(beta,t,y,y0)        % Define objective function
tspan =t;
[tt yy] = ode45(@ModelEqs,tspan,y0,[],beta);
yc= spline(tt,yy,t);
f1=y-yc;
% ==================================
function dydt = ModelEqs(t,y,beta)          % Model equations
dydt = -beta*(0.7-y).^(1/3)*(y-0.61);

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
2楼2012-07-02 17:38:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

★ ★ ★ ★
fegg7502: 金币+4, 应助指数+1, 专家考核, 多谢交流 2012-07-09 08:07:38
结果如下
Local minimum possible.

lsqnonlin stopped because the final change in the sum of squares relative to
its initial value is less than the default value of the function tolerance.




Estimated Parameters by Lsqnonlin():
         k = 1.0138 ± 0.3000
  The sum of the residual squares is: 1.5e-004


The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
3楼2012-07-02 17:39:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yupeifeng

金虫 (小有名气)

引用回帖:
3楼: Originally posted by dbb627 at 2012-07-02 17:39:27
结果如下
Local minimum possible.

lsqnonlin stopped because the final change in the sum of squares relative to
its initial value is less than the default value of the function tolerance.

< ...

首先感谢这位 学长啊,我有点看不懂啊,你帮我 分析下思路,步骤是什么?就是先做设么计算,
憧憬未来
4楼2012-07-02 22:01:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿哈工大材料专硕 085600 初试总分318分 求调剂 +5 aaaa33! 2025-04-04 5/250 2025-04-05 07:22 by 我是小康
[考研] 苏州大学 化学 317求调剂 +11 小宋必定上岸 2025-04-04 13/650 2025-04-05 07:11 by 河大薇薇
[考研] 轻工类 082200 考了数学 +8 Pdz666 2025-03-31 8/400 2025-04-04 22:41 by 诗与自由
[考研] 化学0703调剂293分 +5 太阳花向日葵 2025-04-04 5/250 2025-04-04 22:34 by 有机合成小白100
[考研] 312求调剂 +5 一只大嵩鼠 2025-03-31 5/250 2025-04-04 21:31 by Ade.
[考研] 085602化学工程,273分 +18 青阳21st 2025-04-01 19/950 2025-04-04 19:01 by DNDYR
[考研] 273求调剂 +3 攀登T1 2025-04-04 3/150 2025-04-04 15:43 by lin-da
[考研] 0856(材料工程)368求调剂 +5 Bovenz 2025-04-03 5/250 2025-04-04 07:48 by liangzixue
[考研] 273求调剂 +6 鱼与鱼欲语 2025-04-03 6/300 2025-04-03 21:18 by xxcao001
[考研] 264材料与化工求调剂 +6 泳圈Swim_RING 2025-03-30 7/350 2025-04-03 20:46 by yujianx
[考研] 341求调剂 +4 moonlight! 2025-04-03 4/200 2025-04-03 16:48 by 阳光花
[考研] 273求调剂 +10 攀登T1 2025-04-03 10/500 2025-04-03 16:48 by mini小九
[考研] 344求调剂 +3 萌新吃席 2025-04-02 3/150 2025-04-03 14:55 by fmesaito
[考研] 一志愿哈工大材料与化工,总分390 +4 夜洱yer 2025-04-03 4/200 2025-04-03 14:50 by Ade.
[考研] 333求调剂 +13 无问西东1231213 2025-03-30 13/650 2025-04-03 14:10 by 17313308663
[考研] 来看看我吧材料工程277寻求调剂 +12 游夜天使 2025-04-01 13/650 2025-04-03 10:51 by jean5056
[考研] 材料与化工310求调剂院校 +18 五毛冰淇淋a 2025-04-01 20/1000 2025-04-02 23:19 by 星空星月
[考研] 材料与化工专硕,085600考研求调剂 +21 c一定要上岸 2025-03-31 30/1500 2025-04-02 14:45 by yanenwang
[考研] 284求考研调剂 +4 禾鹤^_^ 2025-04-02 4/200 2025-04-02 14:03 by 李昀睿
[考研] 085602化学工程-299分求调剂 +3 李司煵 2025-03-29 3/150 2025-04-01 10:07 by houyaoxu
信息提示
请填处理意见
澶嶈瘯璋冨墏澶嶈瘯璋冨墏
澶嶈瘯璋冨墏浜岀淮鐮�鑾峰彇澶嶈瘯 & 璋冨墏鎸囧崡锛�寰俊鎵竴鎵�
鍟嗗姟鍚堜綔鍟嗗姟鍚堜綔
鍟嗗姟鍚堜綔浜岀淮鐮�鍟嗗姟鍚堜綔鍜ㄨ寰俊鎵竴鎵�(澶囨敞锛氬晢鍔″悎浣�)
瀛︽湳骞茶揣瀛︽湳骞茶揣
瀛︽湳骞茶揣浜岀淮鐮�寰俊鎵竴鎵�浜嗚В鏈€鏂板鏈共璐�
璁烘枃杈呭璁烘枃杈呭
璁烘枃杈呭浜岀淮鐮�寰俊鎵竴鎵�娣诲姞鑰佸笀锛屽厤璐瑰挩璇�
甯堣祫鎷涘嫙甯堣祫鎷涘嫙