计算方法实验指导书

实验参考书：
《计算方法》,徐萃薇等编著，高等教育，2003.7
《数值分析与实验》，薛毅编著，工业大学，2005.3
实验项目：
A, MATLAB基本操作
容：矩阵操作，基本数学函数，逻辑函数操作等；
要求：完成一些基本练习题
B,插值法
容：Lagrange插值、分段线性插值、三次样条在计算机上用MATLAB软件实现
掌握MATLAB基本操作及矩阵基础知识，包括：输入矩阵、矩阵的转置、矩阵元素求和、矩阵下际、矩阵连接、矩阵行列删除、矩阵产生和操作、逻辑和关系运算、操作符和特殊字符、基本矩阵和矩阵操作、基本矩阵和阵列；掌握特殊变量和常数；掌握基本数学函数。
二、实验容
利用基本矩阵产生3x3和15x8的单位阵，全1阵，全0阵，均匀分布的随机阵（[-1，1]之间），正态分布随机阵（方差4，均值1）
计算方法实验指导书
大学电子信息工程学院
2007年9月
计算方法实验教学计划…………………………………………………………2
实验一MATLAB基本操作……………………………………………………..3
实验二插值法…………………………………………………………………..5
实验三线性方程组的直接解法………………………………………………..7
％Lagrange插值多项式，其中，
％x为向量，全部的插值节点；
％y为向量，插值节点处的函数值；
％xi为标量，被估计函数的自变量：
％yi为xi处的函数估计值．
n＝length(x)；m=1ength(y);
％输入的插值点与它的函数值应有相同的个数
ifn~＝m
error('The length of X must be equal!')；
实验四线性方程组的迭代方法………………………………………………10
实验五函数逼近与曲线拟合…………………………………………….…. 12
计算方法实验教学计划
指导教师：昆仑
学时数：10学时周4学时2次实验，共3周6次实验，第7—9教学周，
每次实验2学时
所用仪器设备：MATLAB7.0实验软件系统
实验指导书：Matlab实验指导书自编
要求：完成一些基本练习题
C,线性方程组的求解——直接解法
容：高斯消元法，列主元素法及其误差分析等在计算机上用MATLAB软件实现。
要求：用实例在计算机上计算和作图。
D,线性方程组的求解——迭代法
容：向量和矩阵的数，雅可比迭代法，高斯-赛德尔迭代法及其收敛性等在计算机上用MATLAB软件实现。
要求：用实例在计算机上计算。
E,函数逼近与曲线拟合
容：曲线拟和与最小二乘方法在计算机上用MATLAB软件实现，并用实例在计算机上计算和作图。
要求：用实例在计算机上计算。
实验一来自百度文库ATLAB基本操作
实验目的
熟悉MATLAB的使用方法及特点；学会建立MATLAB搜索路径；熟悉MATLAB工作空间、MATLAB集成环境、命令窗口；掌握MATLAB的通用命令、管理命令和函数、管理变量和工作空间的使用方法；
5．利用rand函数产生(0,1)间均匀分布的10x10随机矩阵a,然后统计a于等于0.6的元素个数。
6．利用randn函数产生均值为0，方差为1的10x10正态分布随机阵，然后统计其于-0.5，小于0.5的元素个数。
7．编程实现下表功能：
a b
b<=0.5
b>0.5
a<1
显示‘case 1’
显示‘case 2’
return;
end
p=zeros(1,n)；
for k＝1：n
t=ones(1，n)；
for j=1：n
if j~=k
％输入的插值节点必须互异
if abs(x(k)—x(j))< eps
error('the DATA is erroe!');
return;
end
end
p(k)=prod(t);
end
三、实验要求
要求在实验前必须预习，将实验容事先准备好，否则不允许上机。
上机过程中由指导老师检查结果后方可做其他容。每次实验结束后完成实验报告并在下次实验之前由学委统一交给指导教师。
实验二插值法
一、实验目的和要求
(1)学会Langrange插值、Newton插值和Hermite插值等基本插值方法．
(2)学会Matlab提供的插值函数的使用方法，会用这些函数解决实际问题
(3)按照题目要求完成实验容、写出相应的Matlab程序给出实验结果．
(4)对实验结果进行分析讨论．
(5)写出相应的实验报告．
实验容
1. Lagrange插值公式.
Matlab程序：
按照Lagrange插值公式编写Matlab程序(函数名：Lagrange.m)
function yi＝agrange(x，y，xi)
yi=sum(y.*p);
练习1：已知，，，利用Lagrange插值法求的值。
练习2：已知，，，利用Newton插值法求的值，并与上题的结果进行比较。
三、实验要求
要求在实验前必须预习，将实验容事先准备好，否则不允许上机。
上机过程中由指导老师检查结果后方可做其他容。每次实验结束后完成实验报告并在下次实验之前由学委统一交给指导教师。
实验三线性方程组的直接解法
实验目的
(1)学会用Gauss消去法和Matab函数求线性方程组解。
(2)学会对矩阵进行各种分解，如LU分解、Cholskey分解等，并研究相应的Matlab函数。
(3)按照题日要求完成实验容，写出相应的Matlab程序，给出实验结果．对实验结果进行分析讨论．写出相应的实验报告．
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阵变换成行向量。
产生一均匀分布在(-5,5)随机阵(50x2)，精确到小数点后一位。
4．编程实现当α∈[-π，π]，间隔为1o时，求解正弦和余弦的值，并利用plot()函数绘制正弦，余弦曲线。
a>=1
显示‘case 3’
显示‘error’
8．有一矩阵a,找出矩阵中其值大于1的元素，并将他们重新排列成列向量b。
9．在一市区9月份平均气温变化测量矩阵temp_Baoding_sep中(48x30)，存在有奇异值（大于42o C，小于0o C）,编程实现删除奇异值所在的行。
10．在给定的100x100矩阵中，删除整行容全为0的行，删除整列容全为0的列。