四阶含参数矩阵求逆怎么求不出来?。

我在用matlab算矩阵迭代的时候,发现老是各种错误,后来看了一下,发现四阶矩阵在求逆的时候老是失败,这啥意思,三阶的就可以算??这是我算矩阵迭代的程序
clear;clc
format long;
A1=[1;1;1;1];  
eps = 0.000001;   
for i = 1:500  
    syms m k A
    M=[m 0 0 0;0 m 0 0;0 0 m 0;0 0 0 m];
    K=[-k k 0 0;-k 0 k 0;0 -k 0 k;0 0 -k k];
    Q=inv(K);
    D=Q*M;
    B=D*A1;
    A=B./B(4,1);
    C=B./A;
    C=C(1,1);
    if(abs(A-A1) < eps)  
        break;  
    end  
    A1 = A;
    p1=sqrt(1./C);
end  
disp('第一主振型:');  
vpa(A1)  
disp('第一阶固有频率:');  
vpa(p1)
disp('迭代次数:');  
i
然后结果是:
错误使用 mupadmex
Error in MuPAD command: The operand is invalid. [_index]
出错 sym/privBinaryOp (line 820)
            Csym = mupadmex(op,args{1}.s, args{2}.s, varargin{:});
出错  ./  (line 248)
        X = privBinaryOp(A, B, 'symobj::zip', 'symobj::divide');
然后我看了一下求逆的情况:
syms m k A
    M=[m 0 0 0;0 m 0 0;0 0 m 0;0 0 0 m];
    K=[-k k 0 0;-k 0 k 0;0 -k 0 k;0 0 -k k];
    Q=inv(K)
结果是失败的:
Q =
FAIL

 Accepted Answer

0 votes

syms k
K=[-k k 0 0;-k 0 k 0;0 -k 0 k;0 0 -k k];
det(K)
ans =0
奇异矩阵,存在逆矩阵?

More Answers (0)

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!