资源描述：
《用MATLAB作线性和非线性最小二乘法拟合.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、新乡学院数学与信息科学系实验报告实验名称插值与拟合Ⅱ所属课程数学软件与实验实验类型综合型实验专业信息与计算科学班级2011级1班学号姓名李欢丽指导教师朱耀生老师一、实验概述【实验目的】学会用MATLAB作线性和非线性最小二乘法拟合.【实验原理】1.用polyfit作线性最小二乘法拟合:a=polyfit(x,y,m),a表示输出的拟合多现实的系数[a1,a2,…an]（数组）,x,y指输出同长度的数组x,y,m指拟合多项式的次数。2,用lsqcorvefit作非线性最小二乘法拟合x=lsqcurvefit(‘fun’,x0,xda
2、td,ydata)fun是一个事先建立的定义函数f(x,xdata)的M文件，自变量是x,xdata.x0迭代初值。Xdata,ydata已知数据点。【实验环境】MatlabR2010b二、实验内容1．问题一某年美国旧车价格的调查资料如下表，其中xi表示轿车的使用年数，yi表示相应的平均价格。试分析用什么形式的曲线来拟合上述的数据，并预测使用4.5年后轿车的平均价格大致为多少？表1xi12345678910yi26151943149410877655384842902262041.分析问题用polyfit作线性最小二乘法拟合2.问
3、题求解x1=1:10;y=[2615194314941087765538484290226204];A=polyfit(x1,y,2);z1=polyval(A,x1);plot(x1,z1)x=4.5;z=polyval(A,x)3.结果z=955.70474.结论及分析经多次试验可知分析无误问题二增加生产、发展经济所依靠的主要因素有增加投资、增加劳动力以及技术革新等，在研究国民经济产值与这些因素的数量关系时，由于技术水平不像资金、劳动力那样容易定量化，作为初步的模型，可认为技术水平不变，只讨论产值和资金、劳动力之间的关系。在科
4、学技术发展不快时，如资本主义经济发展的前期，这种模型是有意义的。用Q，K，L分别表示产值、资金、劳动力，要寻求的数量关系。经过简化假设与分析，在经济学中，推导出一个著名的Cobb-Douglas生产函数：（*）式中要由经济统计数据确定。现有美国马萨诸塞州1900—1926年上述三个经济指数的统计数据，如下表，试用数据拟合的方法，求出式（*）中的参数。表2tQKLtQKL19001.051.041.0519011.181.061.0819021.291.161.1819031.301.221.2219041.301.271.1719
5、051.421.371.3019061.501.441.3919071.521.531.4719081.461.571.3119091.602.051.4319101.692.511.5819142.013.241.6519152.003.241.6219162.093.611.8619171.964.101.9319182.204.361.9619192.124.771.9519202.164.751.9019212.084.541.5819222.244.541.6719232.564.581.8219242.344.581.
6、6019111.812.631.5919121.932.741.6619131.952.821.6819252.454.581.6119262.584.541.64提示：由于（*）式对参数α,β，a是非线性的，因此，可以有两种方式进行拟合，一是直接使用MATLAB软件中的曲线或曲面拟合命令。另一个是将非线性函数转化成线性函数的形式，使用线性函数拟合1.分析问题用lsqcorvefit作非线性最小二乘法拟合2.问题求解a=[1.041.061.161.221.271.371.441.531.572.052.512.63...2.74
7、2.823.243.243.614.14.364.774.754.544.544.58...4.584.584.54;1.051.081.181.221.171.31.391.471.31...1.431.581.591.661.681.651.621.861.931.961.951.9...1.581.671.821.61.611.64];z=[1.051.181.291.31.31.421.51.521.461.61.691.81...1.931.952.012.002.091.962.22.122.162.082.242.5
8、6...2.342.452.58];x0=[0.10.10.2];x=lsqcurvefit('fun3',x0,a,z)m=linspace(0,2.7,27);n=linspace(0,2.7,27);[M,N]=meshgrid(m,n);z
