四阶含参数矩阵求逆怎么求不出来?。
Show older comments
我在用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
More Answers (0)
Categories
Find more on Common Operations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!