how to plot this equation

1 visualización (últimos 30 días)
shiv gaur
shiv gaur el 15 de Sept. de 2022
Comentada: shiv gaur el 15 de Sept. de 2022
equation is
x^6+k^2*X4+k*x^2+3=0;
by varying the value of k=1:10
plot the graph between real(x) vs k
  2 comentarios
Torsten
Torsten el 15 de Sept. de 2022
For each value of k, you get 6 values for x. So what do you want to plot ? Trace one x-value over the k-range ?
shiv gaur
shiv gaur el 15 de Sept. de 2022
equation is this type in 123and plot is this type 1234
prameter are Ct =0.5 Cr=1.5 wj =0

Iniciar sesión para comentar.

Respuestas (2)

Torsten
Torsten el 15 de Sept. de 2022
Editada: Torsten el 15 de Sept. de 2022
syms x y k
f = y^3 + k^2*y^2 + k*y + 3;
s = solve(f==0,y,'Maxdegree',3);
s1 = sqrt(s);
s2 = -sqrt(s);
k = 0:0.1:10;
s11 = matlabFunction(s1(1));
s12 = matlabFunction(s1(2));
s13 = matlabFunction(s1(3));
s21 = matlabFunction(s2(1));
s22 = matlabFunction(s2(2));
s23 = matlabFunction(s2(3));
plot(k,[real(s11(k));real(s12(k));real(s13(k));real(s21(k));real(s22(k));real(s23(k))])

Sam Chak
Sam Chak el 15 de Sept. de 2022
This only shows the propagation of the complex-valued roots of the 6th-order polynomial equation for the integer .
So, you want to extract the real part of the roots?
k = 1:10;
hold on
for j = 1:length(k)
num = 3;
den = [1 0 k(j)^2 0 k(j) 0 3];
r = rlocus(num, den, 1);
plot(r, 'o')
end
hold off
title('Root Locus Plot')
grid on, xlabel('Real axis'), ylabel('Imaginary axis')
  2 comentarios
Sam Chak
Sam Chak el 15 de Sept. de 2022
Continuation.
Not sure if this is what you want...
k = 1:10;
hold on
for j = 1:length(k)
num = 3;
den = [1 0 k(j)^2 0 k(j) 0 3];
r = rlocus(num, den, 1);
plot(k(j), real(r), 'o')
end
hold off
grid on, xlabel('\it{k}'), ylabel('Re(\it{x})')
shiv gaur
shiv gaur el 15 de Sept. de 2022
sir I have sent fig file as well as equation file pl help this

Iniciar sesión para comentar.

Categorías

Más información sobre Stability Analysis en Help Center y File Exchange.

Etiquetas

Productos


Versión

R2022a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by