24小时热门版块排行榜    

查看: 625  |  回复: 0

candyzhang

新虫 (初入文坛)

[求助] 求助准二级吸附动力学matlab代码

clc;
clear;
X=[10 20 30 40 50 60 80 100 120 150 180 210 270 330 390];
Y1=[5.74359 6.48718 7.53846 7.89744 7.94872 8.25641 8.46154 8.46154 8.66667 8.61538 8.76923 8.82051 9.4359 9.41026 9.35897];
Y2=[4.92308 5.17949 5.58974 6.46154 6.5641  6.51282 7.12821 7.48718 7.74359 8.87179 9.58974 10.25641 10.87179 10.66667 10.89744];
Y3=[6.97436 7.64103 8.51282 9.4359 9.53846 10.10256 10.82051 11.17949 11.79487 12.66667 12.97436 13.14103 13.28205 13.38462 13.17949];


%% function 2
myfun=inline('(x.*A(1)*A(2)^2)./(1+x.*A(1)*A(2))','A','x');
A=nlinfit(X,Y1,myfun,[-3 12.4]);
xx=X;
yy=(xx.*A(1)*A(2)^2)./(1+xx.*A(1)*A(2));
plot(X,Y1,'o',xx,yy,'r');
corrcoef(Y1,yy)
准二级动力学拟合公式t/qt=1/K2qe2+t/qe,我输入进去后算出的相关系数是负值,而且k值绝对大,最后的拟合结果是Qt等于Qe,为什么,到底是哪里错了呢,求大侠给出正确的matlab代码
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 candyzhang 的主题更新
信息提示
请填处理意见