matlab非线性参数拟合,matlab非线性拟合模型参数,大神你在哪! - 计算模拟 - 小木虫 - 学术 科研 互动社区...

CODE:

function fit_nonl

clear all;clc

format long

data=[    8.6667    2.0000

10.0000    5.0000

10.6667   10.0000

12.6667   15.0000

13.3333   30.0000

15.0000   60.0000

15.3333  120.0000

16.6667  200.0000

17.3333  300.0000

18.3333  550.0000

20.0000  700.0000

];

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

texp=data(:,1);    %t的数据,在此输入

k0=[0.1 0.1 1 100];

lb=-[1 1 1 1]*1e9;

ub=[1 1 1 1]*1e9;

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

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

options=optimset('maxfunevals',1000);

[k,resnorm,residual,exitflag,output,lambda,jacobian] = ...

lsqnonlin(@objfunc,k0,lb,ub,options,xspan,texp);

ci = nlparci(k,residual,jacobian);

%residual;

fprintf('\n\n拟合结果:\n')

fprintf('\n\t参数 e1 = %.16f',k(1))

fprintf('\n\t参数 e2 = %.16f',k(2))

fprintf('\n\t参数 e3 = %.16f',k(3))

fprintf('\n\t参数 t = %.16f',k(4))

y=kineticseqs(xspan,k);

r2=1-sum((texp-y).^2)./sum((texp-mean(y)).^2);

fprintf('\n\tr^2 = %.16f',r2);

figure

plot(xspan,kineticseqs(xspan,k),'b',xspan,texp,'or'),legend('计算值','实验值','location','best')

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

function f = objfunc(k,xspan,texp)

f=kineticseqs(xspan,k)-texp;

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

function xt = kineticseqs(x,k)

a=2;

e1=k(1);e2=k(2);e3=k(3);t=k(4);

xt=a/e1 +(a/e2)*(1-exp(-x/t))+a*x/e3;

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值