24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
帖子: 2936
在线: 1329.4小时
虫号: 664177

[交流] 【求助】使用MATLAB拟合出按照公式的系数

我想拟合出y=a(1)*exp(-b(1)*B(i,1))+a(2)*exp(-b(2)*B(i,1))+a(3)*exp(-b(3)*B(i,1))+y(0),其中x,y是已知数，其他的都是需要拟合出来的系数，x,y 来自于这个文件：CT.xls,我编写了一个m文件，不知道为什么不行：

CODE:

function plottu

clc;

B=xlsread('E:\CT.xls','sheet1')

fun=inline('a(1)*exp(-b(1)*B(i,1))+a(2)*exp(-b(2)*B(i,1))+a(3)*exp(-b(3)*B(i,1))+y(0)','a','b','y(0)','x');

a(0)=[10,10,10];

b(0)=[30,30,30];

y(0)=36;

[x,resnorm]=lsqcurvefit(fun,a(1),b(1),y(1),B(i,1),B(i,2));

end

回复此楼

» 猜你喜欢

售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有9人回复
广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有6人回复
两块石头已经有5人回复
面上没中，邀请各位路过的虫友分析一下分数已经有8人回复
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生已经有3人回复
找导师已经有8人回复
国社科又开始会评了，不知道这次命运如何已经有14人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有14人回复
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够已经有17人回复
学科评审组评审是指会评吗？已经有5人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab拟合方程参数时初值的选择已经有15人回复
matlab拟合模型参数已经有5人回复
matlab拟合求值已经有17人回复
在matlab中拟合曲面的表达式已经有29人回复
需要用matlab拟合数据的外行人急需帮助已经有10人回复
matlab拟合公式中未知系数已经有5人回复
采用matlab中nlinfit拟合，出现问题了，不知道是哪错了，该怎么解决，求高手指点一下已经有9人回复
【求助】Matlab非线性最小二乘拟合活度系数模型（Willson、NRTL、UNIQUAC）已经有17人回复
如何用MATLAB拟合指数函数已经有3人回复
【求助】如何使得matlab拟合的曲线强制经过一个点已经有8人回复
【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】已经有15人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-01-11 18:03:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

信彼南山

木虫 (著名写手)

应助: 33 (小学生)
金币: 4142.9
帖子: 1178
在线: 233.5小时
虫号: 1133529

★ ★
robert2020(金币+2):多谢应助！辛苦了！鼓励继续解答！ 2011-01-12 11:06:29

% zyj_main.m
clear;clc;
B=xlsread('CT.xls','sheet1');
xdata=B(:,2);
ydata=B(:,3);
a=[10,10,10];
b=[30,30,30];
y0=36;
x=lsqcurvefit(@zyj_fun1,[a,b,y0],xdata,ydata);

-----------------------------------------------------------------------
% zyj_fun1.m
function zas=zyj_fun1(x,xdata)
zas=x(1)*exp(-x(4)*xdata)+x(2)*exp(-x(5)*xdata)+x(3)*exp(-x(6)*xdata)+x(7)';

----------------------------------------------------------------------
结果似乎还不对，你的这些数据本身就比较接近y=a*exp(-b*x)+c形式