4 符号推演与计算.ppt

第四章 符号推演与计算 4.3 矩阵代数中的符号计算 矩阵代数中的符号计算包括对矩阵进行各种四则运算、求矩阵的行列式、矩阵的逆、矩阵的谱分解以及求线性方程组的解等。部分命令和含义见表8.3.1 注意：有些命令如求矩阵模norm，求矩阵的条件数cond不能用在符号矩阵中 【例4.3.1】求希尔伯特符号矩阵H的特征值、行列式、逆矩阵等矩阵操作。 clear all;clc H = sym(hilb(5)) % 定义一个5阶希尔伯特矩阵 H =[ 1, 1/2, 1/3, 1/4, 1/5] [ 1/2, 1/3, 1/4, 1/5, 1/6] [ 1/3, 1/4, 1/5, 1/6, 1/7] [ 1/4, 1/5, 1/6, 1/7, 1/8] [ 1/5, 1/6, 1/7, 1/8, 1/9] d = det(H) % 求希尔伯特矩阵的行列式的值 d =1/266716800000 X = inv(H) % 求希尔伯特矩阵的行列式的逆矩阵 X = [ 25, -300, 1050, -1400, 630] [ -300, 4800, -18900, 26880, -12600] [ 1050, -18900, 79380, -117600, 56700] [ -1400, 26880, -117600, 179200, -88200] [ 630, -12600, 56700, -88200, 44100] p = poly(H) % 求希尔伯特矩阵的特征多项式 p = x^5-563/315*x^4+735781/2116800*x^3-852401/222264000*x^2+61501/53343360000*x-1/266716800000 e = eig(vpa(H)) % 将符号矩阵转为10位有效数的数字矩阵并求特征根 e = [ .3287939270e-5] [ .3058980252e-3] [ .1140749157e-1] [ .2085342187] [ 1.567050691] t = sym('t'); [I,J] = meshgrid(1:5); H = 1./(I+J-t) % 产生带参数t的希尔伯特矩阵 H=subs(H,t,1) % 将参数替换为1 H = [ 1/(2-t), 1/(3-t), 1/(4-t), 1/(5-t), 1/(6-t)] [ 1/(3-t), 1/(4-t), 1/(5-t), 1/(6-t), 1/(7-t)] [ 1/(4-t), 1/(5-t), 1/(6-t), 1/(7-t), 1/(8-t)] [ 1/(5-t), 1/(6-t), 1/(7-t), 1/(8-t), 1/(9-t)] [ 1/(6-t), 1/(7-t), 1/(8-t), 1/(9-t), 1/(10-t)] H=subs(H,t,1) % 将参数替换为1 H = 1.0000 0.5000 0.3333 0.2500 0.2000 0.5000 0.3333 0.2500 0.2000 0.1667 0.3333 0.2500 0.2000 0.1667 0.1429 0.2500 0.2000 0.1667 0.1429 0.1250 0.2000 0.1667 0.1429 0.1250 0.1111 【例4.3.2】 计算矩阵的秩 syms a11 a12 a13 a21 a22 a23 x y A=[a11, a12; a21, a22]；% 这是一个满秩矩阵 rank(A) ans = 2 B=[a11, a11; a22, a22]；% 第一行的元素相同，第二行的元素也相同 rank(B) ans = 1 C=magic(6) % 产生6阶魔方矩阵 rank(sym(C)) C = 35 1 6 26 19 24 3 32