The Student Room Group
Uni students - Complete our survey >>
Uni students - Complete our survey >>
Uni students - Complete our survey >>

This discussion is now closed.

Check out other Related discussions

Quick Matlab Question relating matrices

Let L and M be square matrices of equal size with known values.

syms x lamda

Z = det (L -lamda*M)

This would give us a polynomial in lamda

c0 + c1*lamda + c2* lamda ^2 ...

Is there a way to find the roots of lamda?

thanks
Reply 1
Avatar for Thusitha
Thusitha
OP
This is what I came up with.

The reason I have to do this is, I can't do an fzero using the function Z for some reason

xmin = 0;
xmax = 5000;
meshSize = 0.01

Result = subs(Z, {lamda}, {xmin:meshSize:xmax});
sizeOfVector = length(Result)
for i=1:sizeOfVector-1
current = Result(i);
next = Result(i+1);
if (current < 0 & next > 0)
i * meshSize
[current next]
end

if (current > 0 & next < 0)
i * meshSize
[current next]
end

if (current == 0)
i * meshSize
[current next]
end
end;
Reply 2
Avatar for Thusitha
Thusitha
OP
Ignore the question. I found what I wanted to do.

There is a thing called generalize eigen function. Which will give me the roots.

[P Q] = eig (L, M);
diag(Q)

Related discussions

Latest

Trending

Trending

Articles for you