计算方法实验指导书

1、计算方法实验指导书河北大学电子信息工程学院2007年9月目录计算方法实验教学计划 2实验一 MATLAB基本操作 .3实验二插值法 .5实验三线性方程组的直接解法 .7实验四线性方程组的迭代方法 10实验五函数逼近与曲线拟合 .12计算方法实验教学计划指导教师：李昆仑学时数：10学时 周4学时2次实验，共3周6次实验，第79教学周， 每次实验2学时所用仪器设备：MATLAB7.0实验软件系统实验指导书：Matlab实验指导书自编实验参考书：计算方法，徐萃薇等编著，高等教育出版社，2003.7数值分析与实验，薛毅编著，北京工业大学出版社，2005.3实验项目：A, MATLAB基本操作内容：矩阵

2、操作，基本数学函数，逻辑函数操作等；要求：完成一些基本练习题B, 插值法内容：Lagra nge插值、分段线性插值、三次样条在计算机上用MATLAB软件实现要求：完成一些基本练习题C, 线性方程组的求解一一直接解法内容：高斯消元法，列主元素法及其误差分析等在计算机上用MATLAB软件实现。要求：用实例在计算机上计算和作图。D, 线性方程组的求解迭代法内容：向量和矩阵的范数，雅可比迭代法，高斯-赛德尔迭代法及其收敛性等在计算机上用MATLAB软件实现。要求：用实例在计算机上计算。E, 函数逼近与曲线拟合内容：曲线拟和与最小二乘方法在计算机上用MATLAB软件实现，并用实例在计算机上计算和作图。要

3、求：用实例在计算机上计算。实验一 MATLAB基本操作实验目的熟悉MATLAB的使用方法及特点；学会建立MATLAB搜索路径；熟悉 MATLAB工作空间、MATLAB集成环境、命令窗口；掌握 MATLAB的通用命令、管理命令和函数、管理变量和工 作空间的使用方法；掌握MATLAB基本操作及矩阵基础知识，包括：输入矩阵、矩阵的转置、矩阵元素求和、矩 阵下际、矩阵连接、矩阵行列删除、矩阵产生和操作、逻辑和关系运算、操作符和特殊字符、 基本矩阵和矩阵操作、基本矩阵和阵列；掌握特殊变量和常数；掌握基本数学函数。二、实验内容利用基本矩阵产生 3x3和15x8的单位阵，全间)，正态分布随机阵(方差 4,均

4、值1 )2 利用diag()函数和rot90()产生下列矩阵：0-731阵，全0阵，均匀分布的随机阵(卜1, 1之0a= 02850b=然后求解a阵的逆矩阵行向量。产生一均匀分布在(-5,5)随机阵(50x2)，精确到小数点后一位。4 .编程实现当a - n,n ，间隔为1o正弦，余弦曲线。5 .禾U用rand函数产生(0,1)间均匀分布的素个数。6.利用randn函数产生均值为0,方差为1小于0.5的元素个数。7 .编程实现下表功能：a baa及b阵的特征值和对应特征向量，并利用reshape将aa阵变换成时，求解正弦和余弦的值，并利用 plot()函数绘制10x10随机矩阵a,然后统计a中

5、大于等于0.6的元的10x10正态分布随机阵，然后统计其中大于-0.5,b0.5a=1显示case显示errora,找出矩阵中其值大于1的元素，并将他们重新排列成列向量temp_Baoding_sep 中(48x30)，存在有奇异值&有一矩阵9.在一保定市区9月份平均气温变化测量矩阵(大于42o C,小于0o C),编程实现删除奇异值所在的行。10.在给定的100x100矩阵中，删除整行内容全为0的行，删除整列内容全为0的列。三、实验要求要求在实验前必须预习，将实验内容事先准备好，否则不允许上机。上机过程中由指导老师检查结果后方可做其他内容。每次实验结束后完成实验报告并在下次实验之前由学委统一

6、交给指导教师。实验二插值法一、实验目的和要求(1) 学会Langrange插值、Newton插值和Hermite插值等基本插值方法.(2) 学会Matlab提供的插值函数的使用方法，会用这些函数解决实际问题按照题目要求完成实验内容、写出相应的Matlab程序给出实验结果.(4)对实验结果进行分析讨论.写出相应的实验报告.实验内容1. Lagrange插值公式.Matlab 程序：按照Lagrange插值公式编写 Matlab程序(函数名：Lagrange.m)function yi = agrange(x, y, xi)% Lagra nge插值多项式，其中，%x为向量，全部的插值节点；%y为

7、向量，插值节点处的函数值；% xi为标量，被估计函数的自变量：% yi为xi处的函数估计值.n = length(x); m=1ength(y);%输入的插值点与它的函数值应有相同的个数ifn = merror(The length of X must be equal!);return;endp=zeros(1, n);for k = 1: nt=ones(1 , n);for j=1 : nif j=k%输入的插值节点必须互异if abs(x(k)x(j)a_max a_max=abs(A(i,k); r=i;endendif a_maxkfor j=k: nz=A(k, j);A(k,j

8、)=A(r,j);A(r,j)=z;end%消元过程for i=k+11 : Hm = A(i, k)/A(k, L);for i = k 十 1:nm=A(i,k)/A(k,k);for j=k+1: nA(i,j)=A(i,j)-m*A(k,j);endb(i)=b(i)-m*b(k);enddet=det*A(k,k);enddet=det*A( n,n);%回代过程if abs(A( n,n) v 1e-10in dex=O;retur n;endfor k=n :-1:1for j=k+1: nb(k)=b(k)-A(k,j)*x(j)endx(k)=b(k)/A(k,k)end练习

9、题1：用列主元法求解线性方程组注意：在Matlab函数中，求解线性方程组的命令非常简单。只需使用矩阵除法。如X=Ab,或利用矩阵的逆 X=i nv(A)*b。而求行列式的函数为det(A).请同学自行比较计算结果。练习题2 :编写Matlab程序，求列矩阵的LU分解。三、实验要求要求在实验前必须预习，将实验内容事先准备好，否则不允许上机。上机过程中由指导老师检查结果后方可做其他内容。每次实验结束后完成实验报告并在下次实验之前由学委统一交给指导教师。实验四线性方程组的迭代解法实验目的(1) 学会用Jacobi迭代法、Gauss- Seidel迭代法和超松弛迭代法求线性方程组解(2) 学会对各种迭

10、代法作收敛性分析，研究求方程组解的最优迭代方法.(3) 按照题目要求完成实验内容，写出相应的Matlab程序，给出实验结果.(4) 对实验结果进行分析讨论.(5) 写出相应的实验报告.二、实验内容1. 熟悉Jacobi迭代法，并编写 Matlab程序matlab程序按照算法(Jacobi迭代法)编写Matlab程序(Jacobi.m)function x,k,i ndex=Jacobi(A,b,ep,it_max)% 求解线性方程组的Jacobi迭代法，其中，% A为方程组的系数矩阵；% b为方程组的右端项;% ep为精度要求，缺省值 1e-5;% it_max为最大选代次数，缺省值 100;

11、% x为方程组的解；% k为迭代次数；% index为指标变量index=1表示迭代收敛到指定要求%,index=0表示迭代失败.if nargin4 it_max=100;e ndif nargin3 ep=1e-5;e ndn=le ngth(A);k=0;x=zeros( n,1);y=x;i ndex=1;while lfor i=1: ny(i)=b(i)for j=1: nif j = iy(i)=y(i)-A(i,j)*x(j);endendif abs(A(i,i)1e-10|k= it_maxin dex=0; retur n;endy(i)=y(i)/A(i,i);endi

12、f no rm(y-x,i nf)epbreak;ende=y;k=k+1;end2. 熟悉Gauss-Seidel迭代法，并编写 Matlab程序3. 练习练习题1.用Jacobi迭代法求方程组的解。练习题2.用Gauss-Seidel迭代法求解上题的线性方程组，取。三、实验要求要求在实验前必须预习，将实验内容事先准备好，否则不允许上机。上机过程中由指导老师检查结果后方可做其他内容。每次实验结束后完成实验报告并在下次实验之前由学委统一交 给指导教师。实验五函数逼近与数据拟合实验目的(1) 学会函数的最佳平方逼近，理解最佳平方逼近的意义.(2) 学会数据的最小二乘拟合.(3) 按照题目要求完成

13、实验内容，写出相应的Matlab程序，给出实验结果.(4) 对实验结果进行分析讨论.写出相应的实验报告.二、实验内容1 .最小二乘算法Matlab程序算法：(1) 输入数据和权系数选择函数逼近系。(2) 求解法方程组。(3) 求出最小二乘逼近。Matlab 程序：编写成Matlab程序(函数名：square_least.m)计算数据最小二乘拟合系数.fun ctio n S=square_least(x, y,n,w)%数据的最小二乘拟合，其中% x,y为数据的(x,y)坐标；%n为数据拟合的次数，缺省时 n=1;% w为权值，缺省值 w=1;%S为数据拟合的系数；%需要另写的外部函数.%phi(x)为基函数，通常为多项式1,x,xA2,.global;global;if nargin 4w = 1; end if nargin 3n = 1; end phi2=zeros (n+1);for i=0: nfor i=0: nphi2(i+1, j+1)=sum(w.* phi_k(x,i).*phi_k(x,j);endendphiF= zeros (n+1,1);endS=phiphiF1 数据的最小二乘拟和假定某天的气温变化记录如表所示 时间123