2016春计算方法实验指导-作业资料.docVIP

计算方法实验指导材料 目 录 第1次 解线性方程组的直接解法 - 2 - 1.1 例题 - 2 - 1.2 Matlab解线性方程组常用命令介绍 - 6 - 1.3 实验作业 - 7 - 第2次 解线性方程组的迭代解法 - 8 - 2.1 例题 - 8 - 2.2 Matlab迭代解法常用函数介绍 - 12 - 2.3 实验作业 - 12 - 第3次 非线性方程求根 - 13 - 3.1 例题 - 13 - 3.2 Matlab非线性方程求根的命令 - 22 - 3.3 实验作业 - 22 - 第4次 插值法 23 4.1 例题 23 4.2 Matlab插值函数介绍 29 4.3 实验作业 30 第5次 利用最小2乘法进行曲线拟合 31 5.1 例题 31 5.2 Matlab数据拟合命令介绍 34 5.3 实验作业 35 第6次 数值积分与数值微分 36 6.1 例题 36 6.2 Matlab数值积分函数介绍 42 6.3 实验作业 43 第 解线性方程组的直接解法 1.1 例题 例.1 用Gauss消元法解方程组： 解：直接建立求解该方程组的M文件Gauss.m如下 % 求解例题.1 % 高斯法求解线性方程组Ax=b % A为输入矩阵系数,b为方程组右端系数 % 方程组的解保存在x变量中 % 先输入方程系数 A=[1 2 3;2 7 5;1 4 9]; b=[1 6 -3]; [m,n]=size(A); %检查系数正确性 if m~=n error(矩阵A的行数和列数必须相同); return; end if m~=size(b) error(b的大小必须和A的行数或A的列数相同); return; end %再检查方程是否存在唯一解 if rank(A)~=rank([A,b]) error(A矩阵的秩和增广矩阵的秩不相同,方程不存在唯一解); return; end %这里采用增广矩阵行变换的方式求解 c=n+1; A(:,c)=b; %%消元过程 for k=1:n-1 A(k+1:n, k:c)=A(k+1:n, k:c)-(A(k+1:n,k)/ A(k,k))*A(k, k:c); End %%回代结果 x=zeros(length(b),1); x(n)=A(n,c)/A(n,n); for k=n-1:-1:1 x(k)=(A(k,c)-A(k,k+1:n)*x(k+1:n))/A(k,k); end %显示计算结果 disp(x=); disp(x); 直接运行上面的M文件或在Matlab命令窗口中直接输入Gauss即可得出结果. 在Matlab命令窗口中输入Gauss得出结果如下: Gauss x= 2.0000 1.0000 -1.0000 扩展： Matlab求解线性方程的几种命令如下（方程组的一般形式可用矩阵和向量表示成，但运用下列方法的前提必须保证所求解的方程为恰定方程，即方程组存在唯一解）． 运用求逆思想：: 或 ； 左除法：原理上是运用高斯消元法求解，但Matlab在实际执行过程中是通过分解法进行的（即先将矩阵A作分解,再回代计算）：； 符号矩阵法：这种计算方法最接近精确值，但计算速度最慢：； 将矩阵施行初等行变换化成行简化阶梯形的办法：可以这样实现之: ; ． 上面四种常用的办法示例如下： A=[1 2 3;2 7 5;1 4 9] % 上面示例方程组系数 A = 1 2 3 2 7 5 1 4 9 b=[1 6 -3] % 方程组右端的系数 b = 1 6 -3 x1_1=inv(A)*b,x1_2=A^(-1)*b % 方法一，求逆思想 x1_1 = 2.0000 1.0000 -1.0000 x1_2 = 2.0000 1.0000 -1.0000 x2=A\b % 方法二，左除思想 x2 = 2 1 -1 x3=sym(A)\sym(b) % 方法三，符号法 x3 = 2 1 -1 C=[A,b],rref(C) % 方法四，行简化阶梯形思想,最后输出结果的一列为解 C = 1 2 3 1 2 7 5 6 1 4 9 -3 ans = 1 0 0