Please help in plotting the graphs
Afficher commentaires plus anciens
syms kp ki w l c alp T ki1 D2 m o L pm
w=10
L=0
b=[1]
a=[0.6 2 1]
%t=tf(b,a)
nb=[0]
na=[0.7 1.2 0]
pm=70
h=fotf(a,na,b,nb)
m=(i*w).^(na)
a1=[0.6;2;1]
%d=m*a1
n=(i*w).^(nb)
b1=[1]
d=(m*a1)/(n*b1);
A0=real(d)
B0=imag(d)
D2=tan(atan(B0/A0)+pm*pi/180+L*w)
fun= @(l) D2/((10^(-l))*sin(l*pi/2)-(10^(-l))*cos(l*pi/2)*D2)
fplot(fun)
hold on
fun1=@(l,ki) ((1/10^(l + 1)*ki*l*sin((pi*l)/2))/(kp*((ki*cos((pi*l)/2))/10^l + 1)) - (1/10^(l + 1)*ki^2*l*cos((pi*l)/2)*sin((pi*l)/2))/(10^l*kp*((ki*cos((l*pi)/2))/10^l + 1)^2))/((1/10^(2*l)*ki^2*sin((pi*l)/2)^2)/(kp^2*((ki*cos((l*pi)/2))/10^l + 1)^2) + 1) + (1/(25*((2^(1/2)*10^(1/2))/5 + 1)^2) - (2^(1/2)*10^(1/2))/(100*((2^(1/2)*10^(1/2))/5 + 1)))/(4/(5*((2^(1/2)*10^(1/2))/5 + 1)^2) + 1)-L
fimplicit(fun1)

this is the graph of fplot(fun) pic above. I know that this graph is is not right. Its not coming properly. Also I need to incorporate another graph [fimplicit(fun1)]. But on giving the fimplicit command no plot of fun1 is coming, matlab is unable to plot a large equation as that.
So help me with both the above problems please. I request you to look into it kindly.
3 commentaires
darova
le 22 Avr 2020
It appears that your function has 3 independent variable

How to plot this madness?
Neelanjan Pal
le 3 Mai 2020
darova
le 4 Mai 2020
I don't have any suggestions. MATLAB did exactly what you told it to do
syms l ki
L=0;
fun1 = ((1/10^(l + 1)*ki*l*sin((pi*l)/2))/((ki*cos((pi*l)/2))/10^l + 1) - (1/10^(l + 1)*ki^2*l*cos((pi*l)/2)*sin((pi*l)/2))/(10^l*((ki*cos((l*pi)/2))/10^l + 1)^2))/((1/10^(2*l)*ki^2*sin((pi*l)/2)^2)/((ki*cos((l*pi)/2))/10^l + 1)^2 + 1);
F = matlabFunction(fun1);
[K,L] = meshgrid(-5:.1:5);
surf(K,L,F(K,L))
axis vis3d

Réponses (0)
Catégories
En savoir plus sur Line Plots dans Centre d'aide et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
