MATLAB机械工程最小二乘法曲线拟合的应用实例
班姓学
级名号
指导教师
f一，实验目的通过Matlab上机编程，掌握利用Matlab软件进行数据拟合分析及数据可视化方法
二，实验内容1有一组风机叶片的耐磨实验数据，如下表所示，其中X为使用时间，单位为小时h，Y为磨失质量，单位为克g。要求：对该数据进行合理的最小二乘法数据拟合得下列数据。x1000011000120001300014000150001600017000180001900020000210002200023000
y24026529832434737741142844647365887513781742
三，程序如下X10000100023000Y24026529832434737741142844647365887513781742dy15for
16aSpolyfitxy
A
a拟合数据y的步长
fdadysqrtdiagi
vSRSRDa
dafreedom
SdfyedeltapolyvalaxSYE
yeD
deltachi2
sumyye2dydye
dQ1chi2cdfchi2freedomclfshgsubplot121plot1：6abschi2freedom‘b’xlabel‘阶次’，title‘chi2与自由度’subplot122plot1：6Q‘r’1：6o
es1605xlabel‘阶次’title‘Q与05线’
odi
put‘根据图形选择适当的阶次（请输入数值）’；elfshgplotxy‘kx’；xlabel‘x’ylabel‘y’；axis8000230002001742；holdo
errorbarxYE
odD
od‘r’；holdofftitle‘较适当阶次的拟合’text100001500‘chi2’
um2strchi2
od‘’i
t2strfreedom
od判断拟合良好度
ftext100001400‘freedom’i
t2strfreedom
odtext20000400‘Q’
um2strQ
od‘05’disp‘’disp‘拟合多项式系数’，dispA
oddisp‘拟合系数的离差’，dispDA
od运行结果分为两个阶段，第一阶段先判断拟合度，第二阶段根据拟合度，选择合适的拟合阶次，再绘出拟合结果。运行结果如下：根据图形选择适当的阶次（请输入数值）4拟合多项式系数0．000000000拟合系数的离差00000000000000000399153313000000016095042295
f四．输出图像
ffr