计算方法实验指导书.doc

计算方法实验指导书河北大学电子信息工程学院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基本操作内容：矩阵操作，基本数学函数，逻辑函数操作等；要求：完成一些基本练习题B, 插值法 内容：Lagrange插值、分段线性插值、三次样条在计算机上用MATLAB软件实现要求：完成一些基本练习题C, 线性方程组的求解直接解法内容：高斯消元法，列主元素法及其误差分析等在计算机上用MATLAB软件实现。要求：用实例在计算机上计算和作图。D, 线性方程组的求解迭代法内容：向量和矩阵的范数，雅可比迭代法，高斯-赛德尔迭代法及其收敛性等在计算机上用MATLAB软件实现。要求：用实例在计算机上计算。E, 函数逼近与曲线拟合内容：曲线拟和与最小二乘方法在计算机上用MATLAB软件实现，并用实例在计算机上计算和作图。 要求：用实例在计算机上计算。实验一 MATLAB基本操作一、 实验目的1 熟悉MATLAB的使用方法及特点；学会建立MATLAB搜索路径；熟悉MATLAB工作空间、MATLAB集成环境、命令窗口；掌握MATLAB的通用命令、管理命令和函数、管理变量和工作空间的使用方法；2 掌握MATLAB基本操作及矩阵基础知识，包括：输入矩阵、矩阵的转置、矩阵元素求和、矩阵下际、矩阵连接、矩阵行列删除、矩阵产生和操作、逻辑和关系运算、操作符和特殊字符、基本矩阵和矩阵操作、基本矩阵和阵列；掌握特殊变量和常数；掌握基本数学函数。二、实验内容1 利用基本矩阵产生3x3和15x8的单位阵，全1阵，全0阵，均匀分布的随机阵（-1，1之间），正态分布随机阵（方差4，均值1）2利用diag()函数和rot90()产生下列矩阵： 0 0 8 2 0 4 a= 0 -7 5 b= 0 5 0 2 3 0 7 0 8然后求解a阵的逆矩阵aa及b阵的特征值和对应特征向量，并利用reshape将aa阵变换成行向量。3 产生一均匀分布在(-5,5)随机阵(50x2)，精确到小数点后一位。 4编程实现当-，间隔为1o时，求解正弦和余弦的值，并利用plot()函数绘制正弦，余弦曲线。5利用rand函数产生(0,1)间均匀分布的10x10随机矩阵a,然后统计a中大于等于0.6的元素个数。 6利用randn函数产生均值为0，方差为1的10x10正态分布随机阵，然后统计其中大于-0.5，小于0.5的元素个数。7编程实现下表功能：a b b0.5a=1显示case 3显示error8有一矩阵a,找出矩阵中其值大于1的元素，并将他们重新排列成列向量b。9在一保定市区9月份平均气温变化测量矩阵temp_Baoding_sep中(48x30)，存在有奇异值（大于42o C，小于0o C）,编程实现删除奇异值所在的行。10在给定的100x100矩阵中，删除整行内容全为0的行，删除整列内容全为0的列。三、实验要求要求在实验前必须预习，将实验内容事先准备好，否则不允许上机。上机过程中由指导老师检查结果后方可做其他内容。每次实验结束后完成实验报告并在下次实验之前由学委统一交给指导教师。实验二 插值法一、实验目的和要求(1)学会Langrange插值、Newton插值和Hermite插值等基本插值方法(2)学会Matlab提供的插值函数的使用方法，会用这些函数解决实际问题(3) 按照题目要求完成实验内容、写出相应的Matlab程序给出实验结果(4)对实验结果进行分析讨论(5)写出相应的实验报告二、 实验内容1. Lagrange插值公式.Matlab 程序：按照Lagrange插值公式编写Matlab程序(函数名：Lagrange.m)function yiagrange(x，y，xi)Lagrange插值多项式，其中，x为向量，全部的插值节点；y为向量，插值节点处的函数值；xi为标量，被估计函数的自变量：yi为xi处的函数估计值nlength(x)；m=1ength(y);输入的插值点与它的函数值应有相同的个数ifnm error(The length of X must be equal!)；return;endp=zeros(1,n)；for k1：n t=ones(1，n)； for j=1：n if j=k 输入的插值节点必须互异 if abs(x(k)x(j)a_max a_max=abs(A(i,k); r=i; end end if a_maxk for j=k:n z=A(k，j);A(k,j)=A(r,j);A(r,j)=z; end 消元过程 for i=k+11：H mA(i，k)A(k，L); for ik十1:n m=A(i,k)/A(k,k); for j=k+1:n A(i,j)=A(i,j)-m*A(k,j); end b(i)=b(i)-m*b(k); end det=det*A(k,k); end det=det*A(n,n); 回代过程if abs(A(n,n)1e-10 index=0;return;endfor k=n:-1:1 for j=k+1:n b(k)=b(k)-A(k,j)*x(j) end x(k)=b(k)/A(k,k)end练习题1：用列主元法求解线性方程组注意：在Matlab函数中，求解线性方程组的命令非常简单。只需使用矩阵除法。如X=Ab,或利用矩阵的逆X=inv(A)*b。而求行列式的函数为det(A).请同学自行比较计算结果。练习题2：编写Matlab程序，求列矩阵的LU分解。三、实验要求要求在实验前必须预习，将实验内容事先准备好，否则不允许上机。上机过程中由指导老师检查结果后方可做其他内容。每次实验结束后完成实验报告并在下次实验之前由学委统一交给指导教师。实验四 线性方程组的迭代解法一、 实验目的(1) 学会用Jacobi迭代法、Gauss- Seidel迭代法和超松弛迭代法求线性方程组解(2) 学会对各种迭代法作收敛性分析，研究求方程组解的最优迭代方法(3) 按照题目要求完成实验内容，写出相应的Matlab程序，给出实验结果(4) 对实验结果进行分析讨论(5) 写出相应的实验报告二、实验内容1.熟悉Jacobi迭代法，并编写Matlab程序matlab程序按照算法(Jacobi迭代法)编写Matlab程序(Jacobi.m)functionx,k,index=Jacobi(A,b,ep,it_max) 求解线性方程组的Jacobi迭代法,其中, A为方程组的系数矩阵; b为方程组的右端项; ep为精度要求，缺省值1e-5; it_max为最大选代次数，缺省值100; x为方程组的解; k为迭代次数; index为指标变量 index=1表示迭代收敛到指定要求, , index=0表示迭代失败.if nargin4 it_max=100;endif nargin3 ep=1e-5;endn=length(A);k=0;x=zeros(n,1);y=x;index=1;while l for i=1:ny(i)=b(i) for j=1:n if j = i y(i)=y(i)-A(i,j)*x(j); end end if abs(A(i,i)1e-10|k= it_max index=0; return; end y(i)=y(i)/A(i,i); end if norm(y-x,inf)ep break; end e=y;k=k+1;end2. 熟悉Gauss-Seidel迭代法，并编写Matlab程序3.练习练习题1. 用Jacobi迭代法求方程组的解。练习题2. 用Gauss-Seidel迭代法求解上题的线性方程组，取。三、实验要求要求在实验前必须预习，将实验内容事先准备好，否则不允许上机。上机过程中由指导老师检查结果后方可做其他内容。每次实验结束后完成实验报告并在下次实验之前由学委统一交给指导教师。实验五 函数逼近与数据拟合一、 实验目的(1) 学会函数的最佳平方逼近，理解最佳平方逼近的意义(2) 学会数据的最小二乘拟合(3) 按照题目要求完成实验内容，写出相应的Matlab程序，给出实验结果(4)对实验结果进行分析讨论(5)写出相应的实验报告二、实验内容1最小二乘算法Matlab程序算法：(1) 输入数据和权系数选择函数逼近系。(2) 求解法方程组。(3)求出最小二乘逼近。Matlab程序：编写成Matlab程序(函数名：square_least.m)计算数据最小二乘拟合系数.function S=square_least(x,y,n,w)数据的最小二乘拟合,其中x,y为数据的(x,y)坐标;n为数据拟合的次数，缺省时n=1;w为权值，缺省值w=1;S为数据拟合的系数;需要另写的外部函数.phi(x)为基函数,通常为多项式1,x,x2,.global;global;if nargin 4w = 1; endif nargin 3n = 1; endphi2=zeros(n+1);for i=0:n for i=0:n phi2(i+1, j+1)=sum(w.* phi_k(x,i).*phi_k(x,j); endendphiF= zeros(n+1,1);endS=phiphiF 1数据的最小二乘拟和假定某天的气温变化记录如表所示时间(t)0123456