资源描述：
《重庆大学数学实验报告04》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、实验项目类型课程实验项目数学实验MATLAB软件入门名称名称验证演示综合设计其他指导成绩教师一、实验目的及意义[1]了解最小二乘拟合的基本原理和方法；[2]掌握用MATLAB作最小二乘多项式拟合和曲线拟合的方法；[3]通过实例学习如何用拟合方法解决实际问题，注意与插值方法的区别。[4]了解各种参数辨识的原理和方法；[5]通过范例展现由机理分析确定模型结构，拟合方法辨识参数，误差分析等求解实际问题的过程；通过该实验的学习，掌握几种基本的参数辨识方法，了解拟合的几种典型应用，观察不同方法得出的模型的准确程度，学习参数的误差分析，进一步了解数学建模过程。这对于学生深入理
2、解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。二、实验内容1．用MATLAB中的函数作一元函数的多项式拟合与曲线拟合，作出误差图；2．用MATLAB中的函数作二元函数的最小二乘拟合，作出误差图；3．针对预测和确定参数的实际问题，建立数学模型，并求解。三、实验步骤1．开启软件平台——MATLAB，开启MATLAB编辑窗口；2．根据各种数值解法步骤编写M文件3．保存文件并运行；4．观察运行结果(数值或图形)；5．根据观察到的结果写出实验报告，并浅谈学习心得体会。四、实验要求与任务根据实验内容和步骤，完成以下具体实验，要求写出实验报
3、告（实验目的→问题→数学模型→算法与编程→计算结果→分析、检验和结论→心得体会）应用实验1．旧车价格预测某年美国旧车价格的调查资料如下表，其中xi表示轿车的使用年数，yi表示相应的平均价格。试分析用什么形式的曲线来拟合上述的数据，并预测使用4.5年后轿车的平均价格大致为多少？表1xi12345678910yi2615194314941087765538484290226204解：问题分析：已知xi表示的是轿车的使用年数，yi表示的是相应的平均价格，a表示拟合多项式的系数，dl表示多项式在在x=4.5处的拟合值。可运用拟合将各个系数以及dl求出，再比较理论值与实际值
4、偏离的程度（误差）可以不断精确4.5年后轿车的平均价格。模型假设与建立：运用两组数据进行比较，即4次多项式和9次多项式的比较。已知数据可在matlab中输入以下程序：xi=1:10;yi=[2615194314941087765538484290226204];a=polyfit(xi,yi,4)dl=polyval(a,4.5)e=yi-dl;el=sum(e.*e)或：xi=1:10;yi=[2615194314941087765538484290226204];a=polyfit(xi,yi,9)dl=polyval(a,4.5)e=yi-dl;el=sum
5、(e.*e)模型求解：运行程序后得到如下数据：adlel4次多项式1.0e+003*923.08276.0255e+0060.0001-0.00540.1002-0.90823.42329次多项式1.0e+004*6.0254e+006-0.00000.0001-0.00230.0260-0.1838923.16960.8231-2.32183.9496-3.69051.6611比较数据后可知9次多项式的误差较小，所以应优先选择9次多项式的数据。即4.5年后的轿车平均价格为：923.1696元。模型检验：由以上表格可以看出，取平均价格为：923.1696的数据误差
6、小，基本复合要求。2．机器人识别定形工具柄问题机器人在不同层次上应用于工业生产、水下探测、核点开发、军事研究等领域和部门。当一个机器人工作时，经常需要识别那些从外形上看来是圆形或椭圆形的仪器或工具柄等基本设备，以便执行进一步的操作。通常在所需操纵的工具柄上放置适当数量的传感器，这些传感器不断向四周发射电信号，机器人身上安置有接收电信号的硬件装置，根据这些信号，机器人将估算出各个传感器当时所在的位置，然后，再利用这些数据获得工具柄的位置。由于硬件设备的限制和测量的随机偏差，所获得的传感器位置数据是有误差的。因此，为了增强识别的准确性和可靠性，工具柄上放置的传感器应多
7、于确定该定形曲线所需的最少点数。（能否获得比较准确的工具柄位置，对机器人能否有效抓握、操作该工具柄起着关键的作用。）现有一个圆形工具柄，其边缘上放置了6个传感器，一机器人在某一个时刻测得这些传感器的位置坐标为：(1,7),(2,6),(5,8),(7,7),(9,5),(3,7),如何确定该圆形工具柄的圆心坐标和半径。解：用拟合的方法确定该圆形工具柄的圆心坐标（x(1),x(2）)和半径x(3)M文件程序：functionf=yuan(x,tdata,cdata)tdata=[125793];cdata=[768757];f=(tdata-x(1)).^2+(cd
8、ata-x
