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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 用Matlab求反应动力学参数
11/1返回列表
查看: 1804  |  回复: 0
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

amenboss

铁虫 (初入文坛)

[求助] 用Matlab求反应动力学参数

中间先拟合为一段一级和一段二级，并求出各点的导数，最后拟合出需要的参数。拟合过程中出现导数为负数，求各位大神指导！
代码如下：
clc;
clear all;
format short g
%%%分段求解  一级
expE1=[0        2054.794521
300        903.6835616
600        449.3047945];

t1=expE1(:,1);
C0=expE1(1,2);
C1=expE1(:,2);
LCO=log(C0);
LC=log(C1);
d=LCO-LC;

K0=[0.02];
K=lsqcurvefit(@FunS,K0,t1,d);

mdcdt1=K.*C1;

%%%%分段求解 二级

expE2=[900        200.3150685
1200        107.0773973
1500        78.67123288
2100        57.00136986
2700        47.67328767
3300        38.26164384
3900        28.4
4500        21.78219178];
C02=expE2(1,2);
C2=expE2(:,2);
t2=expE2(:,1);
RC2=1./C2;
d2=RC2-1/C02;

K20=[0.02];
K2=lsqcurvefit(@FunS2,K20,t2,d2);

mdcdt2=K2.*(C2.^2);
%%%%合并求参数
mdcdt=[mdcdt1;mdcdt2];
C=[C1;C2];

%%%% T=250C下，mnts初值P0
P0=[89 0.0328 110 2284.62585];
P=lsqcurvefit(@FunP,P0,C,mdcdt)

mdcdt3=(P(1)+P(2).*C)./(P(3)+1/P(4)./C./C)
figure
plot(C,mdcdt,'r-');
hold on;
plot(C,mdcdt3,'b*')


function f=FunS(K,t1)
f=K(1).*t1;
function f=FunS2(K2,t2)
f=K2(1).*t2;

function f=FunP(P,C)
f=(P(1)+P(2).*C)./(P(3)+1/P(4)./C./C);
回复此楼

» 猜你喜欢
1楼 2014-12-13 15:06:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

科研从小木虫开始，人人为我，我为人人
相关版块跳转 我要订阅楼主 amenboss 的主题更新
11/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭