sumarize8_5_2.m

来自「matlab7.0光盘内容。都是一些具体的例子。」· M 代码 · 共 13 行

H = hilb(3)                                    %生成三阶希尔伯特数值矩阵
H = sym(H)                                   %将数值矩阵转换成为符号矩阵
inv(H)                                        %求符号矩阵的逆矩阵
det(H)                                        %求符号矩阵的行列式
V = vpa(H,16)                                          %用16位精度转换成数值矩阵
inv(V)                                                %对带有误差的矩阵进行求逆
det(V)                                                %求带有误差的矩阵的行列式
syms s;
H(1,1) = s                                            %定义基本符号变量并替换H(1,1)
Z = det(H)
sol = solve(Z)
H = subs(H,s,sol)
inv(H)