计算方法作业2讲述.docxVIP

《计算方法》上机指导书  实验1 MATLAB基本命令 1．掌握MATLAB的程序设计 实验内容：对以下问题，编写M文件。 生成一个5×5矩阵，编程求其最大值及其所处的位置。 编程求 。 一球从100米高度自由落下，每次落地后反跳回原高度的一半，再落下。求它在第10次落地时，共经过多少米？第10次反弹有多高？ 2．掌握MATLAB的绘图命令 实验内容：对于自变量的取值属于[0，]，在同一图形窗口画出如下图形。 （1）； （2）；  实验2 插值方法与数值积分 1. 研究人口数据的插值与预测 实验内容：下表给出了从1940年到1990年的美国人口，用插值方法推测1930年、1965年、2010年人口的近似值。 美国人口数据 年194019501960197019801990人口(千人)132,165151,326179,323203,302226,542249,6331930年美国的人口大约是123,203千人，你认为你得到的1965年和2010年的人口数字精确度如何？ 2．最小二乘法拟合经验公式 实验内容：某类疾病发病率为‰和年龄段(每五年为一段，例如0~5岁为第一段，6~10岁为第二段……)之间有形如的经验关系，观测得到的数据表如下 1234567890.8982.383.071.842.021.942.222.774.02101112131415161718194.765.466.5310.916.522.535.750.661.681.8（1）用最小二乘法确定模型中的参数和。 （2）利用MATLAB画出离散数据及拟合函数图形。 3. 复化求积公式 实验内容：对于定积分。 （1）分别取利用复化梯形公式计算，并与真值比较。再画出计算误差与之间的曲线。 （2）取[0，1]上的9个点，分别用复化梯形公式和复化辛普森公式计算，并比较精度。  实验3 非线性方程与线性方程组 1．矩阵的范数与条件数 实验内容：已知矩阵 求，，和。 2．研究高斯消去法的数值稳定性 实验内容：设方程组，其中 （1）， （2）， 分别对以上两个方程组 （1）计算矩阵的条件数，判断系数矩阵是良态的还是病态的？ （2）用列主元消去法求得L和U及解向量； （3）用不选主元的高斯消去法求得L和U及解向量； （4）观察小主元并分析对计算结果的影响。 3. 求解非线性方程，比较不同方法的计算量 实验内容：比较求的根到三位小数所需的计算量： （1）在区间［0,1］内用二分法； （2）用迭代法，初值； （3）用牛顿迭代法，取初值。  《计算方法》上机实验报告 姓名： 陶成川 学号： U201410820 班级： 机械09 问题 1. 研究人口数据的插值与预测 实验内容：下表给出了从1940年到1990年的美国人口，用插值方法推测1930年、1965年、2010年人口的近似值。 美国人口数据 年194019501960197019801990人口(千人)132,165151,326179,323203,302226,542249,6331930年美国的人口大约是123,203千人，你认为你得到的1965年和2010年的人口数字精确度如何？ 2．最小二乘法拟合经验公式 实验内容：某类疾病发病率为‰和年龄段(每五年为一段，例如0~5岁为第一段，6~10岁为第二段……)之间有形如的经验关系，观测得到的数据表如下 1234567890.8982.383.071.842.021.942.222.774.02101112131415161718194.765.466.5310.916.522.535.750.661.681.8（1）用最小二乘法确定模型中的参数和。 （2）利用MATLAB画出离散数据及拟合函数图形。 3. 复化求积公式 实验内容：对于定积分。 （1）分别取利用复化梯形公式计算，并与真值比较。再画出计算误差与之间的曲线。 （2）取[0，1]上的9个点，分别用复化梯形公式和复化辛普森公式计算，并比较精度。  Matlab程序 %构造lagrange插值函数 function y1=lagrange(x,y,x1) m=length(x); n=length(y); p=length(x1); if m~=n error; end s=0; for k=1:n t=ones(1,p); for j=1:n if