数值计算方法I试验手册.DOC

计算方法 实　验　手　册 江西财经大学信息管理学院 数学与决策科学系 2003.1  目 录 预备实验　MATLAB使用练习…………………………………………2 实验一　数值计算中误差的传播规律…………………………………6 实验二　数值计算中的算法稳定性……………………………………8 实验三 Gauss消去法主元素的选取与算法的稳定性………………10 实验四　观察Runge现象和对非光滑函数进行插值的可能性…13 实验五　观察及改善最小二乘拟合的数值不稳定现象…………15 实验六 曲线逼近方法的比较………………………………………17 实验七 利用数值微分方法求解偏微分方程………………………18 实验八　比较不同数值积分公式对积分方程解的影响……………20 实验九　高维积分数值计算的Monte-Carlo方法…………………23 实验十　观察前进Euler公式的收敛性和数值不稳定性…………26 实验十一　利用实际例子筛选有效算法……………………………28 实验十二　估计线性方程组的性态与条件数………………………29 实验十三　迭代法的算法设计与比较………………………………31 实验十四　迭代法收敛阶的估计……………………………………33 实验十五　迭代法过程中的分支与混沌……………………………34 综合实验一……………………………………………………………35 综合实验二……………………………………………………………36 预备实验　MATLAB使用练习 一、实验目的 MATLAB是本实验课的主要软件，尽快熟悉其基本命令和语法，以适应后续实验对MATLAB的使用要求． １．掌握MATLAB中常用的基本命令和函数； ２．掌握MATLAB的简单作图方法； ３．掌握MATLAB的简单编程技巧． 二、实验内容 １.矩阵和向量的输入和运算 在MATLAB的Commond窗口中输入以下命令，并仔细观察执行后的结果． （1）矩阵的输入： A=[1 2 3;4 5 6] （分号表示矩阵换行，表示回车） A=[1,2,3;4,5,6] （空格或逗号表示同一行元素的分隔） A=[1 2;3 4]; （在一行的末尾加分号，运行结果将不显示） Z=zeros(2,3) （阶零矩阵） E=ones(2) （阶全1矩阵） I=eye(2,4) （阶对角线元素为1的矩阵） R=rand(2,3) （阶(0,1)均匀分布随机矩阵） H=hilb(4) （4阶hilbert矩阵） 另外还有标准正态分布矩阵randn(m,n),n阶幻方矩阵magic(n)等. （2）矩阵元素的调用和赋值： a=A(1,3) （矩阵A的第１行第３列元素赋值给a） A(1,3) （查看矩阵A的第１行第３列元素） A(1,3)=0 （将A的第１行第３列元素赋值为０） A(3,4)=1 （注意矩阵A原本是阶矩阵） B(2,3)=2 （注意前面没有给出矩阵B） A （显示矩阵A） A(2,: （A的第２行） A(:,3) （A的第3列） C1=A(1:2,:) （将A的第１,２行赋值给矩阵C1） C２=A(:,２:４) （将A的第２～４列赋值给矩阵C2） C３=A(２:３,２:４) （将A的第２,３行,第２～４列赋值给C3） C4= A(1:２:３,4:-1:2) （将A的第1、３行,第４、３、２列赋值给C3） A （显示矩阵A为） B1=[A,C2] A和C2为B1） B2=[A；C1] （上下拼接矩阵A和C1为B2） （４）矩阵的运算 MATLAB为矩阵的运算提供了下列矩阵运算符： ＋（加法）；－（减法）； ＇（转置）；＊（乘法）；＾（乘幂）；＼（左除）；/（右除） .＊（点乘法）；.＾（点乘幂）；.＼（点左除）；./（点右除） 矩阵的“＋”、“－”、“＇”、“＊”、“＾”必须遵循矩阵的运算规则，点运算 “.＊”、“.＾”、“.＼”、“./”是指矩阵对应元素之间进行相应的运算． 在MATLAB的Commond窗口中继续输入以下命令，并仔细观察结果： C=zeros(3) D=ones(3) 比较C+1和C+D． A=[1 2 ２;4