用Matlab做线性拟合方法
简介
(01)非线性最小二乘优化在曲线拟合、参数估计等问题中有着广泛的应用。例如,我们要拟合一系列观测数据(t,y),拟合函数为F(t,x),他是x的非线性函数。对于这种最小二乘曲线拟合问题,可以通过Matalb优化工具箱中的lsqcurvefit命令求解,可以根据实际问题进行曲线拟合。
例题
(01)在工程实验中,测得下面一组数据。求系数a、b、c、d,使得函数为表中数据的最佳拟合函数。f(t)=a+b·sin(t)+c·cos(t)+dt3
观测数据表
(01)——————————————————————————————————t | 0 0.5 1 1.5 2 2.5 3 3.5 4——————————————————————————————————y | 0 3.4 4.1 4.6 5.9 6.9 8.1 9.8 11——————————————————————————————————
操作方法
(01)首先建立拟合函数M文件如下:
(02)function f=example8_15(x,ti)n=length(ti);for i=1:nf(i)=x(1)+x(2)*sin(ti(i))+x(3)*cos(ti(i))+x(4)*ti(i)^3;end
(03)从命令窗口输入
(04)>> ti=[0 0.5 1 1.5 2 2.5 3 3.5 4];>> yi=[0 3.4 4.1 4.6 5.9 6.9 8.1 9.8 11];>> x0=[1 1 1 1]x27;; %初始点选为全1向量>> [x,resnorm,residual,exitflag,output,lambda,J]=lsqcurvefit(@example8_15,x0,ti,yi)
(05)输出结果为
(06)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Optimization completed because the size of the gradient is less thanthe default value of the function tolerance.<stopping criteria details>x =1.87062.7714-1.04770.1708resnorm =2.9080residual =Columns 1 through 70.8228 -1.0989 -0.2927 0.5373 0.2929 0.1372 -0.1897Columns 8 through 9-0.5977 0.3887exitflag =1output =firstorderopt: 6.6428e-08iterations: 2funcCount: 15cgiterations: 0algorithm: 'trust-region-reflective'message: [1x425 char]lambda =lower: [4x1 double]upper: [4x1 double]J =(1,1) 1.0000(2,1) 1.0000(3,1) 1.0000(4,1) 1.0000(5,1) 1.0000(6,1) 1.0000(7,1) 1.0000(8,1) 1.0000(9,1) 1.0000(2,2) 0.4794(3,2) 0.8415(4,2) 0.9975(5,2) 0.9093(6,2) 0.5985(7,2) 0.1411(8,2) -0.3508(9,2) -0.7568(1,3) 1.0000(2,3) 0.8776(3,3) 0.5403(4,3) 0.0707(5,3) -0.4161(6,3) -0.8011(7,3) -0.9900(8,3) -0.9365(9,3) -0.6536(2,4) 0.1250(3,4) 1.0000(4,4) 3.3750(5,4) 8.0000(6,4) 15.6250(7,4) 27.0000(8,4) 42.8750(9,4) 64.0000%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
(07)再在命令窗口中输入:
(08)>> xi=0:0.1:4;>> y=example8_15(x,xi);>> plot(ti,yi,'r*')>> grid on>> hold on>> plot(xi,y)>> legend('观测数据点','拟合曲线')>> title('最小二乘曲线拟合')
(09)输出结果如下图所示:
-
怎么样校验Win10 ISO的SHA1 如何验证MD5值
如果你想了解自己下载的Windows10ISO文件是否为原版,这时该怎么样操作呢?我们可以通过验证ISO文件的SHA1与MD5值,下面就来看看如何验证SHA1与MD5的方法吧。操作方法(01)首先我们在下载Windows10ISO文件的页面,可以看到MD5与SHA1的数值,把这些数值复制保存起来。(02)接下...
-
如何使用PPT输入化学式
今天我给大家尤其式学生们介绍一个怎么使用PPT输入化学式的传统方法,希望能帮到你们,一起来看看吧!操作方法(01)首先你需要做的就是打开一个新建的文本,在上方的菜单栏里面选择插入选项卡,看到右边选择绘制横排文本框。(02)线使用鼠标选中文本框,右键单击鼠标,选择跳出选...
-
QQ通讯录怎么添加联系人?
QQ通讯录怎么添加联系人?操作方法(01)我们首先登录自己的QQ号。(02)然后点击软件右下角的应用管理器按钮,如图:(03)点击通讯录按钮,如图:(04)在通讯录中点击添加按钮,如图所示:(05)输入联系人的姓名等信息后点击确定按钮,如图所示:...
-
如何将CorelDraw文件转为jpg图片
CorelDraw保存的源文件是CDR格式,将文件转为jpg格式的图片。操作方法(01)准备好CorelDraw文件(02)然后打开CorelDraw软件进入界面,选择菜单栏中的-文件-打开,(03)找到准备好的cdr文件,将文件置入到文档中(04)当导出时,点击菜单栏中的-文件-导出(05)在保存类型中找到jpg格式,点...