MATLAB: How can i solve this system

To determine the damping ratio (psij and psik ), i have to solve the attached system where all other name

are constant.

How can i solve this system ?

Best Answer

A = 2*omegaj*omegak/(omegaj^2-omegak^2)*[Y*I/omegak -Y*I/omegaj ; -m*omegak m*omegaj];
b = [cs*I ; ca];
zeta = A\b;
zetaj = zeta(1);
zetak = zeta(2);

Best wishes

Torsten.

Related Solutions

MATLAB: I want to find the determinant for 5.5 matrix

>> syms s1 s2 s3 s4 s5 s6 s7 s8 s9 s10
>> M = [0 -s1 -s2 -s3 -s4; s1 0 -s5 -s6 -s4; s2 s5 0 -s8 -s9; s3 s6 s8 0 -s10; s4 s7 s9 s10 0];
>> det(M)
ans =
s1*s4^2*s8^2 - s1*s4*s7*s8^2 - s2*s4^2*s6*s8 + s3*s4^2*s5*s8 - s3^2*s4*s5*s9 - s2^2*s4*s6*s10 + s3^2*s5*s7*s9 + s2^2*s6*s7*s10 + s2*s3*s4*s5*s10 + s2*s3*s4*s6*s9 + s1*s2*s4*s8*s10 - s1*s3*s4*s8*s9 - s2*s3*s5*s7*s10 - s2*s3*s6*s7*s9 + s2*s4*s6*s7*s8 - s3*s4*s5*s7*s8 - s1*s2*s7*s8*s10 + s1*s3*s7*s8*s9
>> simplify(ans)
ans =
(s4 - s7)*(s3*s5 - s2*s6 + s1*s8)*(s2*s10 - s3*s9 + s4*s8)

MATLAB: SUbstitions

M = @(csi) [1, csi(1), csi(2), csi(3), csi(4), csi(1)^2 - 1,...
            csi(1)*csi(2), csi(1)*csi(3), csi(1)*csi(4),...
            csi(2)^2 - 1, csi(2)*csi(3), csi(2)*csi(4),...
            csi(3)^2 - 1, csi(3)*csi(4), csi(4)^2 - 1];
s=[1 2 3 4];
M(s)  % You can use this directly or assign the results to an array.

% EDIT

If you want to be able to do it all at once for n-by-4, then:

M = @(csi) [ones(size(csi,1),1), csi(:,1), csi(:,2), csi(:,3), csi(:,4),...
            csi(:,1).^2 - 1,csi(:,1).*csi(:,2), csi(:,1).*csi(:,3),...
            csi(:,1).*csi(:,4), csi(:,2).^2 - 1, csi(:,2).*csi(:,3),...
            csi(:,2).*csi(:,4), csi(:,3).^2 - 1, csi(:,3).*csi(:,4),...
            csi(:,4).^2 - 1];
s = rand(15,4);
M(s)

Best Answer

Related Solutions

MATLAB: I want to find the determinant for 5.5 matrix

MATLAB: SUbstitions

Related Question