Графический матричный график для петли

Question

Я очень сомневаюсь, что мой код верен для моего задания. Ниже приведено краткое описание задачи с некоторым синтаксисом matlab, но не интерпретируйте ее как код. Мой код можно найти в конце этого сообщения.

M = [100 0;0 10] 

K = [k+oa^2m -oa^2m;-oa^2m oa^2m] 

, где k=8950000, m=10 и oa изменяется в течение петли.

A = inv([-(100pi)^2*M+K]) 

матричные элементы A11 и A12 должна быть построена в интервале 30pi<=oa<=150pi. Тогда я должен найти, где A11=0.

Причина, по которой я скептически отношусь к моему коду (см. Ниже), заключается в том, что график не пересекает ось x, но имеет больше асимптотического характера (y = 0). Нет «чистого» A11=0. Что-то интересное также происходит в oa=~69.46pi. Может ли кто-нибудь найти ошибку (ы)? Благодаря!

M = [100 0;0 10]; 
A11 = []; 
A12 = []; 
oav = linspace(0.3*100*pi,1.5*100*pi,10000); 
for oa=oav; 
A = inv(-(100*pi)^2*M+[8950000+oa^2*10 -oa^2*10; -oa^2*10 oa^2*10]); 
A11 = [A11 A(1,1)]; 
A12 = [A12 A(1,2)]; 
end 
plot(oav,A11,oav,A12) 

Answer 1

A11 делает крест ноль, если вы увеличите достаточно:

>> min(A11) 
ans = -0.0079732 

69.46 < 30*pi так oa=~69.46 даже не в диапазоне вы определили ...