matlab拟合4参数的公式,matlab自定义公式拟合

clc,clear

%定义一个符号变量t

syms t

%写入要拟合数据

x=[0;1;3;10;30;100;300;1000;3000;5000;7000;10000];

y=[0;0.003915;-0.01044;-0.01184;-0.03737;-0.14699;-0.41946;-1.00311;-1.94593;-2.56821;-3.07218;-3.65212];

%编辑要拟合的公式,设置变量,设置系数

ft=fittype(‘a(1-exp(-(t/tau)^beta))‘,‘independent‘,‘t‘,‘coefficients‘,{‘a‘,‘tau‘,‘beta‘});

options = fitoptions(ft);

options.StartPoint=[1 10 1];

%设定拟合所用系数的上下限

options.Lower = [-999999 -99999 0];

options.Upper=[100 99999 1];

%进行拟合

cfun=fit(x,y,ft,options)

%显示拟合函数,数据必须为列向量形式

xi=0:10:10000;

yi=cfun(xi);

figure

plot(x‘,y‘,‘r‘,xi,yi,‘b-‘);

title(‘拟合函数图形‘);

得到的拟合参数:

cfun =

General model:

cfun(t) = a*(1-exp(-(t/tau)^beta))

Coefficients (with 95% confidence bounds):

a = -5.807 (-6.909, -4.705)

tau = 1.021e+04 (5945, 1.448e+04)

beta = 0.7357 (0.6822, 0.7892)

拟合效果图:

2f0626b3f64459dec68a95fff1353433.png

原文:https://www.cnblogs.com/lucaschan/p/11690383.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值