9楼: Originally posted by 山鹰HPU at 2014-11-07 18:34:46
求教您个问题，根据一组数据，下面这个关系式可以拟合出系数R和k吗？可以给予金币，谢谢。
y=5.2/(5.2+u)*exp(57.2-11*u)*erf((45*R-u*x)/(2*(2.1*R*x)^0.5))+ 5.2/(5.2+u)*exp(57.2-11*u)*erf((45*R+u*x)/(2*(2. ...

function feixianxingnihe_3
clear all;clc
format long


data=[298,1.51941000000000;473,1.47000000000000;673,1.38000000000000;873,1.32000000000000;1073,1.28000000000000];
xspan=data(:,1);    %x的数据，在此输入
Texp=data(:,2);    %T的数据，在此输入

k0=[10 10];
lb=-[1 1]*1e9;
ub=[1 1]*1e9;


%-------------------------------------------------------------------------

% 使用函数lsqnonlin()进行参数估计

OPTIONS=optimset('MaxFunEvals',1000);
[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...
    lsqnonlin(@ObjFunc,k0,lb,ub,OPTIONS,xspan,Texp);


fprintf('\n\n拟合结果:\n')
fprintf('\n\t参数 A = %.16f',k(1))
fprintf('\n\t参数 a = %.16f',k(2))

y=KineticsEqs(xspan,k);
R2=1-sum((Texp-y).^2)./sum((Texp-mean(y)).^2);
fprintf('\n\t相关系数之平方R^2 = %.16f',R2);
figure(1)
mm=max(xspan)-min(xspan);
xspan1=min(xspan)-0.05*mm:0.01:max(xspan)+0.05*mm;
plot(xspan1,KineticsEqs(xspan1,k),'b',xspan,Texp,'or'),legend('计算值','实验值','Location','Best')


%-------------------------------------------------------------------------

function f = ObjFunc(k,xspan,Texp)
f=KineticsEqs(xspan,k)-Texp;

%------------------------------------------------------------------------
function xt = KineticsEqs(x,k)
global b
A=k(1);a=k(2);
xt=A/(3*137.8)+A./x.^0.5.*2/a^0.5.*atan((-0.5+0.5*(1+a*137.8./x.^2).^0.5).^0.5)-2*2^0.5/3*A*137.8^0.5./x.^1.5./(1+(1+a*137.8./x.^2).^0.5).^1.5;

拟合结果:

        参数 A = 502.1696785220468100
        参数 a = 23729.5421378699030000
        相关系数之平方R^2 = 0.9274532902868037>>

均方差(RMSE): 0.0241425766168092
残差平方和(SSE): 0.0029143200284925
相关系数(R): 0.963175718167507
相关系数之平方(R^2): 0.927707464067493
决定系数(DC): 0.927453222207639
卡方系数(Chi-Square): 0.00102747119728203
F统计(F-Statistic): 38.4980600874317

参数        最佳估算
----------        -------------
a                  502.170022602279
a1                 23730.0107275906

====== 结果输出 =====

No        实测值y        计算值y
1        1.51941        1.5410049
2        1.47        1.4301342
3        1.38        1.3659321
4        1.32        1.3285206
5        1.28        1.3042527

2楼: Originally posted by 月只蓝 at 2014-07-20 19:57:04
一、MATLAB代码和结果如下，图形结果见附图1。

function feixianxingnihe_3
clear all;clc
format long


data=;
xspan=data(:,1);    %x的数据，在此输入
Texp=data(:,2);    %T的数据，在此输入

k0= ...

3楼: Originally posted by qj2006 at 2014-07-20 20:04:45
谢谢你，我这还有一组数据，麻烦你再给我拟合一次。
采用关系式：
y=A/(3*137.8)+A/x^0.5*2/a^0.5*atan((-0.5+0.5*(1+a*137.8/x^2)^0.5)^0.5)-2*2^0.5/3*A*137.8^0.5/x^1.5/(1+(1+a*137.8/x^2)^0.5)^1.5
x=;
y ...