matlab基础教程实验报告,matlab基础性实验报告5.doc

matlab基础性实验报告5.doc

一、实验要求了解插值与拟合的基本原理和方法；掌握用MATLAB计算插值与作最小二乘多式拟合和曲线拟合的方法；通过范例展现求解实际问题的初步建模过程；通过动手作实验学习如何用插值与拟合方法解决实际问题，提高探索和解决问题的能力。这对于学生深入理解数学概念，掌握数学的思维方法，熟悉处理大量的工程计算问题的方法具有十分重要的意义。二、实验仪器、设备或软件电脑,MATLAB软件三、实验内容1．编写插值方法的函数M文件；2．用MATLAB中的函数作函数的拟合图形；3．针对实际问题，试建立数学模型，并求解。四、实验步骤1．开启软件平台MATLAB，开启MATLAB编辑窗口；2．根据各种数值解法步骤编写M文件；3．保存文件并运行；4．观察运行结果数值或图形；5．写出实验报告，并浅谈学习心得体会。五、实验要求与任务根据实验内容和步骤，完成以下具体实验，要求写出实验报告(实验目的→问题→数学模型→算法与编程→计算结果→分析、检验和结论→心得体会)。1．天文学家在1914年8月的7次观测中，测得地球与金星之间距离(单位米)，并取得常用对数值，与日期的一组历史数据如下表日期(号)18202224262830距离对数996177995436994681993910993122992319991499由此推断何时金星与地球的距离(米)的对数值为993518解输入命令DAYS18202224262830DISTANCELOGS996177995436994681993910993122992319991499T1INTERP1DISTANCELOGS,DAYS,993518线性插值T2INTERP1DISTANCELOGS,DAYS,993518, NEAREST 最近邻点插值T3INTERP1DISTANCELOGS,DAYS,993518, SPLINE 三次样条插值T4INTERP1DISTANCELOGS,DAYS,993518, CUBIC 三次插值计算结果T1249949T224T3250000T4250000综上所得，可推断25日金星与地球的距离(米)的对数值为993518。2．在某海域测得一些点X,Y处的水深Z由下表给出，船的吃水深度为5英尺，在矩形区域(75，200)(50，150)里的哪些地方船要避免进入。XYZ12914010358818551951057514152314722513758554868688XYZ157510757781162162117565813565665843359988949(1)输入插值基点数据；(2)在矩形区域(75，200)(50，150)作二维插值；(3)作海底曲面图；(4)作出水深小于5的海域范围，即Z5的等高线。解程序输入插值基点数据X12914010358818551951051575107577811621621175Y7514152314722513758556581356566584335Z48686889988949ZZ在矩形区域(75，200)(50，150)作二维插值CX7505200CY5005150CZGRIDDATAX,Y,Z,CX,CY , CUBIC 作海底曲面图SUBPLOT1,2,1,MESHZCX,CY,CZXLABEL X ,YLABEL Y ,ZLABEL Z 作出水深小于5的海域范围，即Z5的等高线SUBPLOT1,2,2,C,HCONTOURCX,CY,CZCLABELC,H,5插值后作出的海底曲面图及等高线图如下若船的吃水深度为5英尺，在矩形区域(75，200)(50，150)里如上图等高线5M内的地方船要避免进入。3．用电压V10伏的电池给电容器充电，电容器上T时刻的电压为，其中V0是电容器的初始电压，是充电常数。试由下面一TETV0组(T，V)数据确定V0和。T秒051234579V伏636648726822866899943963解一(1)用命令LSQCURVEFIT。1编写M文件CURVEFUN1MFUNCTIONFCURVEFUN1X,TDATAF1010X1EXPTDATA/X22)主程序XITITHREE1M如下TDATA051234579CDATA636648726822866899943963X004316,1XLSQCURVEFIT CURVEFUN1 ,X0,TDATA,CDATAFCURVEFUN1X,TDATA3)运行主程序，得结果为X5557735002F6149066616749138114785832893539398796604即拟合得V055577，35002。(2)用命令LSQNONLIN。1编写M文件CURVEFUN2MFUNCTIONFCURVEFUN2XTDATA051234579CDATA636648726822866899943963FCDATA1010X1EXPTDATA/X22)主程序XITITHREE2M如下X004316,1XLSQNONLIN CURVEFUN2 ,X0FCURVEFUN2X3)运行主程序，得结果为X5557735002F0211001816023130105300768005470031300304结果同上，即拟合得V055577，35002。解二(1)对将要拟合的非线性模型，建立M文件VOLUMMTEVTV0如下FUNCTIONYHATVOLUMBETA,TYHAT1010BETA1EXPT/BETA2(2)输入数据T051234579Y636648726822866899943963BETA053 (3)求回归系数BETA,R,JNLINFITT ,Y , VOLUM ,BETA0BETA得结果BETA5557735002即得回归模型为5023710TETV(5)预测及作图YY,DELTANLPREDCI VOLUM ,T ,BETA,R,JPLOTT,Y, K ,T,YY, R 012345678966577588599510