Я пытаюсь извлечь мнимые термины из некоторых собственных значений.Я установил, что символика является реальной, и затем я использовал imag () для извлечения результата, но в выводе есть функции imag ().
k=1;
syms a mu xt sigma X r real;
C = (.01365+.2808*k*1i-k^2/2)/(.01365+.3455*k*1i-k^2);
lh = 1 - 2*1i*C/k;
lt = -a - 1i/k - 2*C/k^2 - 2*1i*(.5-a)*C/k;
mh = -a + 2*1i*(.5+a)*C/k;
mt = 1/8 + a^2 - 1i*(.5-a)/k + 2*(.5+a)*C/k^2 + 2*1i*(.25-a^2)*C/k;
trash = [mu*(1-sigma^2*X)+lh, mu*xt+lt; mu*xt+mh, mu*r^2*(1-X)+mt];
detrash = det(trash) == 0;
eigtrash = solve(detrash, X);
imag(soltrash(1))
Я ожидаю символическое выражение без imag ()функции.
Вот что я получаю:
ans =
- (36028797018963968 * imag ((- ^ 4 * сигма ...... этодлинное выражение, но вы можете увидеть изображение там.