试验6数据拟合及参数辨识方法

1、6 数据拟合及参数辨识方法一、实验目的及意义1 了解最小二乘拟合的基本原理和方法；2掌握用MATLAB乍最小二乘多项式拟合和曲线拟合的方法；3 通过实例学习如何用拟合方法解决实际问题，注意与插值方法的区别。4 了解各种参数辨识的原理和方法；5 通过范例展现由机理分析确定模型结构，拟合方法辨识参数，误差分析等求解实际问题的过程；通过该实验的学习，掌握几种基本的参数辨识方法，了解拟合的几种典型应用，观察不同方法得出的模型的准确程度，学习参数的误差分析，进一步了解数学建模过程。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。二、实验内容1 .用MA

2、TLA井的函数作一元函数的多项式拟合与曲线拟合，作出误差图；2 .用MATLA井的函数作二元函数的最小二乘拟合，作出误差图；3针对预测和确定参数的实际问题，建立数学模型，并求解。三、实验步骤3 开启软件平台一一 MATLAB开启 MATLA编辑窗口；2 .根据各种数值解法步骤编写M文件3保存文件并运行；4观察运行结果（ 数值或图形） ；5根据观察到的结果写出实验报告，并浅谈学习心得体会。四、实验要求与任务根据实验内容和步骤，完成以下具体实验，要求写出实验报告（实验目的-问题-数学 模型-算法与编程-计算结果-分析、检验和结论-心得体会）应用实验1 旧车价格预测某年美国旧车价格的调查资料如下表，

3、其中xi 表示轿车的使用年数，yi 表示相应的平均价格。试分析用什么形式的曲线来拟合上述的数据，并预测使用4.5 年后轿车的平均价格大致为多少?表1Xi12345678910V、2615194314941087765538484290226204由题意知用matlab编程：t=1:1:10;r=2615 1943 1494 1087 765 538 484 290 226 204;aa=polyfit(t,r,2);a=aa(1)b=aa(2)c=aa(3)y=polyval(aa,t);plot(t,r,'k+',t,y,'r')z=a*4.5*4.5+b*4

4、.5+c图像如下：z = 955.7047从而可知第4.5年的预测结果为 956辆2 .机器人识别定形工具柄问题机器人在不同层次上应用于工业生产、水下探测、核点开发、军事研究等领域和部门。当一个机器人工作时，经常需要识别那些从外形上看来是圆形或椭圆形的仪器或工具柄等基本设备，以便执行进一步的操作。 通常在所需操纵的工具柄上放置适当数量的传感器，这些传感器不断向四周发射电信号，机器人身上安置有接收电信号的硬件装置，根据这些信号，机器人将估算出各个传感器当时所在的位置，然后，再利用这些数据获得工具柄的位置。由于硬件设备的限制和测量的随机偏差，所获得的传感器位置数据是有误差的。因此，为了增 强识别的

5、准确性和可靠性，工具柄上放置的传感器应多于确定该定形曲线所需的最少点数。(能否获得比较准确的工具柄位置，对机器人能否有效抓握、操作该工具柄起着关键的作用。)现有一个圆形工具柄， 其边缘上放置了 6个传感器，一机器人在某一个时刻测得这些传 感器的位置坐标为：(1,7),(2,6),(5,8),(7,7),(9,5),(3,7),如何确定该圆形工具柄的圆心坐标和半径。用matlab编程如下：xdata=1,2,5,7,9,3'ydata=9,6,8,7,5,7'B=-(xdata.A2+ydata.A2);A=xdata,ydata,ones(6,1);X=AB;a=-1/2*X(

6、1)b=-1/2*X(2)r=sqrt(-X(3)+1/4*(X(1)A2+X(2)A2)alpha=0:pi/20:2*pi;R=r;x=a+R*cos(alpha);y=b+R*sin(alpha);plot(xdata,ydata,'o' ,x,y,'r')axis equal 运行结果：a =4.8007b =6.2908r =3.1901图像：3 .经济增长模型增加生产、发展经济所依靠的主要因素有增加投资、增加劳动力以及技术革新等，在研究国民经济产值与这些因素的数量关系时，由于技术水平不像资金、劳动力那样容易定量化，作为初步的模型，可认为技术水平不变，

7、 只讨论产值和资金、 劳动力之间的关系。 在科学技 术发展不快时，如资本主义经济发展的前期，这种模型是有意义的。用Q, K, L分别表示产值、资金、劳动力，要寻求的数量关系Q(K,L)。经过简化假设与分析，在经济学中，推导出一个著名的Cobb-Douglas生产函数：Q(K,L) = aK “L° , 0< a , 3 <1(*)式中a , 3 , a要由经济统计数据确定。现有美国马萨诸塞州19001926年上述三个经济指数的统计数据，如下表，试用数据拟合的方法，求出式(*)中的参数“，3 , a。表2tQKLtQKL19001.051.041.0519142.013.2

8、41.6519011.181.061.0819152.003.241.6219021.291.161.1819162.093.611.8619031.301.221.2219171.964.101.9319041.301.271.1719182.204.361.9619051.421.371.3019192.124.771.9519061.501.441.3919202.164.751.9019071.521.531.4719212.084.541.5819081.461.571.3119222.244.541.6719091.602.051.4319232.564.581.8219101.692.511.5819242.344.581.6019111.812.631.5919252.454.581.6119121.932.741.661