Matlab最小二乘法曲线拟合的应用实例_第1页
Matlab最小二乘法曲线拟合的应用实例_第2页
Matlab最小二乘法曲线拟合的应用实例_第3页
Matlab最小二乘法曲线拟合的应用实例_第4页
Matlab最小二乘法曲线拟合的应用实例_第5页
已阅读5页，还剩1页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、MATLAB机械工程最小二乘法曲线拟合的应用实例班 级:姓 名:学 号:指导教师: 一，实验目的 通过Matlab上机编程，掌握利用Matlab软件进行数据拟合分析及数据可视化方法 二，实验内容 1.有一组风机叶片的耐磨实验数据，如下表所示，其中X为使用时间，单位为小时h，Y为磨失质量，单位为克g。要求：对该数据进行合理的最小二乘法数据拟合得下列数据。 x=10000 11000 12000 13000 14000 15000 16000 17000 18000 1

2、9000 20000 21000 22000 23000; y=24.0 26.5 29.8 32.4 34.7 37.7 41.1 42.8 44.6 47.3 65.8 87.5 137.8 174.2 三，程序如下X=10000:1000:23000;Y=24.0,26.5,29.8,32.4,34.7,37.7,41.1,42.8,44.6,47.3,65.8,87.5,137.8,174.2dy=

3、1.5; %拟合数据y的步长for n=1:6 a,S=polyfit(x,y,n); An=a; da=dy*sqrt(diag(inv(S.R´*S.R);Dan=da´freedom(n)=S.df;ye,delta=polyval(a,x,S);YEn=ye;Dn=delta;chi2(n)=sum(y-ye).2)/dy/dy;end Q=1-chi2cdf(chi2,freedom); %判断拟合良好度clf,shgsubplot(1,2,1),plot(1：6,abs(chi2-freedom),b)xlabel(阶次)，title(chi2与自由度)subp

4、lot(1,2,2),plot(1：6,Q,r,1：6,ones(1,6)*0.5)xlabel(阶次),title(Q与0.5线)nod=input(根据图形选择适当的阶次（请输入数值）)；elf,shg,plot(x,y,kx)；xlabel(x),ylabel(y)；axis(8000,23000,20.0,174.2)；hold onerrorbar(x,YEnod,Dnod,r)；hold offtitle(较适当阶次的拟合)text(10000,150.0,chi2=num2str(chi2(nod)int2str(freedom(nod)text(10000,140.0,freedom=int2str(freedom(nod)text(20000,40.0,Q=num2str(Q(nod)0.5)disp()disp(拟合多项式系数)，disp(Anod)disp(拟合系数的离差)，disp(DAnod)运行结果分为两个阶段，第一阶段先判断拟合度，第二阶段根据拟合度，选择合适的拟合阶次，再绘出拟合结果。运行结果如下：>>根据图形选择适当的阶次（请输入数值）4拟合多项式系数00000 -0.000

人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论