搜索
实验四医用薄膜渗透率的确定——数据拟合一、实验目的及意义[1]了解最小二乘拟合的基本原理和方法；[2]掌握用MATLAB作最小二乘多项式拟合和曲线拟合的方法；[3]通过实例学习如何用拟合方法解决实际问题，注意与插值方法的区别。[4]了解各种参数辨识的原理和方法；[5]通过范例展现由机理分析确定模型结构，拟合方法辨识参数，误差分析等求解实际问题的过程；通过该实验的学习，掌握几种基本的参数辨识方法，了解拟合的几种典型应用，观察不同方法得出的模型的准确程度，学习参数的误差分析，进一步了解数学建模过程。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。二、实验内容1．用MATLAB中的函数作一元函数的多项式拟合与曲线拟合，作出误差图；2．用MATLAB中的函数作二元函数的最小二乘拟合，作出误差图；3．针对预测和确定参数的实际问题，建立数学模型，并求解。三、实验步骤1．开启软件平台——MATLAB，开启MATLAB编辑窗口；2．根据各种数值解法步骤编写M文件3．保存文件并运行；4．观察运行结果(数值或图形)；5．根据观察到的结果写出实验报告，并浅谈学习心得体会。四、实验要求与任务根据实验内容和步骤，完成以下具体实验，要求写出实验报告（实验目的→问题→数学模型→算法与编程→计算结果→分析、检验和结论→心得体会）应用实验1．旧车价格预测某年美国旧车价格的调查资料如下表，其中xi表示轿车的使用年数，yi表示相应的平均价格。试分析用什么形式的曲线来拟合上述的数据，并预测使用4.5年后轿车的平均价格大致为多少？表1xi12345678910yi2615194314941087765538484290226204用matlab软件作出旧车价格随使用年数变化的散点图xdata=1:1:10;ydata=[2615,1943,1494,1087,765,538,484,290,226,204]plot(xdata,ydata,'ro')做出半对数图形，可直观的看出车价随使用年限服从负对数分布,可用对数函数模型进行拟合xdata=1:1:10;ydata=[2615,1943,1494,1087,765,538,484,290,226,204];semilogy(xdata,ydata,'+')，可直观的看出车价随使用年限服从负对数分布,可用对数函数模型进行拟合编写程序：xdata=1:1:10;ydata=[2615,1943,1494,1087,765,538,484,290,226,204]k=polyfit(xdata,log(ydata),1)y=polyval(k,xdata);z=exp(y);subplot(1,2,1),plot(xdata,y,xdata,log(ydata),'ro'),gtext('logy=-0.296*x+8.1591'),title('拟合曲线');subplot(1,2,2),plot(xdata,z,xdata,ydata,'k+'),title('拟合曲线');运行y=polyval(k,4.5);z=exp(y)结果：z=918.6543可得4.5年车龄的汽车价格大概为918.65432．机器人识别定形工具柄问题机器人在不同层次上应用于工业生产、水下探测、核点开发、军事研究等领域和部门。当一个机器人工作时，经常需要识别那些从外形上看来是圆形或椭圆形的仪器或工具柄等基本设备，以便执行进一步的操作。通常在所需操纵的工具柄上放置适当数量的传感器，这些传感器不断向四周发射电信号，机器人身上安置有接收电信号的硬件装置，根据这些信号，机器人将估算出各个传感器当时所在的位置，然后，再利用这些数据获得工具柄的位置。由于硬件设备的限制和测量的随机偏差，所获得的传感器位置数据是有误差的。因此，为了增强识别的准确性和可靠性，工具柄上放置的传感器应多于确定该定形曲线所需的最少点数。（能否获得比较准确的工具柄位置，对机器人能否有效抓握、操作该工具柄起着关键的作用。）现有一个圆形工具柄，其边缘上放置了6个传感器，一机器人在某一个时刻测得这些传感器的位置坐标为：(1,7),(2,6),(5,8),(7,7),(9,5),(3,7),如何确定该圆形工具柄的圆心坐标和半径。根据已知点，用圆方程进行拟合222()()xaybr进行拟合首先建立M文件myfun2.m：functionf=myfun2(x)xdata=[125793];ydata=[768757];f=(xdata-x(1)).^2+(ydata-x(2)).^2-x(3)^2end运行x0=[0.10.10.1];x=lsqnonlin('myfun2',x0)可得结果a=4.7423b=3.8351r=4.1088建立M文件：zuotu.mxdata=[125793];ydata=[768757];ezplot('(x-4.7423)^2+(y-3.8385)^2=4.1088^2',[-4,10]);holdon,plot(xdata,ydata,'or'),axisequal得拟合圆与已知点图形