matlab求四次方程的解,matlab解四次方程的问题 - 程序语言 - 小木虫 - 学术 科研 互动社区...

CODE:

function k = res2

syms x

w = 0.9 : 0.01 : 1.1; %21个

gamma = 1.1 : 0.5 : 10;   %18个

c = 3 * 10^8;

R = 2.3 * 10^-2;

n1 = length(w);

n2 = length(gamma);

v = zeros(n2);

y = zeros(n1);

wpb = zeros(n2);

f = sym(zeros(n1, n2));

k = zeros(n1 * n2, 1);

for j = 1 : 18        %gamma

v(j) = c * sqrt(1 - 1 ./ gamma(j).^2);

for i = 1 : 21        %w

y(i) = 1 + 0.2^2 ./ (1 - w(i).^2 + 2 * i * 0.015 * w(i));

wpb(j) = sqrt((1.602 * 10^-19)^2 * 10^9 ./ (9.109 * 10^-31 * 8.85 * 10^-12 * gamma(j)^3));

f(i, j) = (y(i) * (w(i) * 2 * pi * 24 * 10^9)^2 / c^2 - x^2) * (1 - (wpb(j) * 2 * pi * 24 * 10^9)^2 / (((w(i) * 2 * pi * 24 * 10^9) - x * v(j))^2 * y(i))) - 2.4048^2 / R^2;

k((i-1) * n2 + j) = max(sqrt(eval(solve(f(i, j)))));

end

end

semilogy(gamma, k(1:18), 'b--');

figure;

plot(gamma, k(1:18), 'r');