北京科技大学matlab综合实验,北京科技大学《matlab和数学实验》第三次

1、矩阵的基本运算 解线性方程组 用数值方法计算定积分 矩阵特征值、特征向量 矩阵的基本运算 注意 k是一个数， A是一个矩阵 k*A AB AXB, XA-1B, A必须是方阵 数 乘 矩阵的左除 矩阵的右除 A/B XBA,XAB-1, B必须是方阵 矩阵的行列式 detA A必须为方阵 矩阵的逆 invA A必须为方阵， |A| 0 矩阵的乘幂 An A必须为方阵， n是正整数 矩阵行变换化简 rrefA 求 A阶梯形的行最简形式 矩阵的特征值、特征向量、特征多项式 V,DeigA 例 1 A1,-1;2,4; V,DeigA ans V -985/1393 1292/2889 985/13。

2、93 -2584/2889 方阵 A的特 征向量矩阵 D 2 0 0 3 方阵 A的特 征值矩阵 矩阵的特征值、特征向量、特征多项式 ppolyA 若 A为矩阵，则 p为 A的特征多项式系数； 若 A为行向量，则 p为以 A为根的特征多项式系数。 例 1 A1,-1;2,4; ppolyA poly2strp,x poly2strp,x 得到多项式的习惯形式 ans p1 -5 6 x2-5x6 解 线 性 方 程 组 1、逆矩阵法(求逆法) X 1.4000 0.4000 解 例 1 求方程组的解 2 3 4 1 xy xy A2,3;1,-1; b4;1 XinvA*b 相当于 2 3 。

3、4 1 1 1 x y ans 方程的解是 x1.4, y0.4 P94 A2,3;1,-1; b4;1 XAb 逆矩阵法(左除与 右 除法) 例 1 求方程组的解 2 3 4 1 xy xy 解 线 性 方 程 组 解 ans X 1.4000 0.4000 方程的解是 x1.4, y0.4 相当于 AXb,XAb 2、初等变换法 解 线 性 方 程 组 在线性代数中用消元法求线性方程组的通解的过程为 1、 用初等变换化线性方程组为阶梯形方程组，把最 后的恒等式“ 00” 去掉； 2、如果剩下的方程当中最后的一个等式是零等于非 零的数，那么方程无解。否则有解； 3、在有解的情况下 如果阶梯。

4、形方程组中方程的个数 r等于未知量 的个数，那么方程组有唯一的解； 如果阶梯形方程组中方程的个数 r小于未知量 的个数，那么方程组有无穷多个解。 例 8 求齐次线性方程组的通解 1 2 3 4 1 2 3 4 1 2 3 4 8 1 0 2 0 2 4 5 0 3 8 6 2 0 x x x x x x x x x x x x 解 Matlab命令为 1 0 4 0 0 1 -3/4 -1/4 0 0 0 0 ans A1 -8 10 2;2 4 5 -1;3 8 6 -2; 系数矩阵 rrefA 行的最简形式 解 线 性 方 程 组 分析 将 00的一行去掉，则原方程组等价于 13 2 3 。

5、4 4 31 44 xx x x x 方程的个数 未知量个数 有无穷多个解 取 3 4 1 3 x x 3 4 0 4 x x 得 1 2 4 0 x x 取 得 1 2 0 1 x x 基础解系为 1 4 0 1 3 ， 2 0 1 0 4 所以方程的通解为 1 2 12 3 4 40 01 10 34 x x kk x x 其中 k1, k2 是任意实数 解 线 性 方 程 组 例 9 求非齐次方程的解 1 2 3 1 2 3 12 4 2 2 3 1 2 1 0 1 1 3 8 x x x x x x xx 解 Matlab命令为 A4 2 -1;3 -1 2;11 3 0; b2;10。

6、;8; BA,b 增广矩阵 系数矩阵 常数项 rrefB 解 线 性 方 程 组 ans 1 0 3/10 0 0 1 -11/10 0 0 0 0 1 结果分析 行最简形式中最后一行出现了零等于 非零的情况，故方程组无解。 解 线 性 方 程 组 用数值方法计算定积分 y x a b yfx 的几何意义 b a f x d x 有三种方法 1、 矩形法 sumx 2、复合梯形公式 3、复合辛普生公式 1 20 4 , 1 dxx 例 1 计算定积分 与精确值 比较。 h0.01;x0h1; y4./1x.2; at long z1sumy1100*h 左矩形公式 z2sumy2101*h 右。

7、矩形公式 解 MATLAB命令为 1、使用矩形法求定积分 输出结果 z1 3.151575986923129 z2 3.131575986923129 u1 0.0100 u2 -0.0100 at short u1z1-pi,u2z2-pi 2、复合梯形公式 用小梯形面积代替小曲边梯形的面积，然后求和 以获得定积分的近似值，比矩形法精度高。 命令 trapzx,y 相当于 求 1 1 2 n i i y i y i x 3、复合辛普生公式 用抛物线代替小曲边梯形的曲边计算小面积，然后 求和以获得定积分的近似值，精度比前两种方法高。 命令 quadfun,a,b,tol,tr。

8、ace 1、式中 fun是被积函数表达式字符串或者是 M函数 文件； 2、 a,b是积分的下限与上限； 3、 tol代表精度，可以缺省( tol0.001)； 4、 trace1时用图形展示积分过程，省略时无图形。 例 3 用三种方法计算定积分 21 0 sin 1 x dx x 的值。 解 编程如下 x00.011; ysinx.2./x1; s1sumy1100*0.01 s2sumy2101*0.01 s3trapzx,y ffinlinesinx.2./x1,x s4quadff,0,1 s4 0.1808 运行结果为 s1 0.1787 s2 0.1829 s3 0.1808 0.1808. 按左矩形公式计算结果是 0.1787， 按右矩形公式计 算结果是 0.1829， 按梯形法和辛普生法计算结果都是 建立函数文件 jifen.m function yjifenx ysinx.2./x1; squadjifen,0,1 编程如下 作业 第七章习题 P167 17. 2 18 P115 21. 1,2 P114 14 名称线性代数相关运算及数值方法计算定积分 目的掌握矩阵的基本运算、特征值、特征向量和线 性方程组的求解；能熟练运用数值方法求定积分。