四阶含参数矩阵求逆怎么求不出来?。
Mostra commenti meno recenti
我在用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
Risposta accettata
Più risposte (0)
Categorie
Scopri di più su Common Operations in Centro assistenza e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!