Please help in plotting the graphs
Mostra commenti meno recenti
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 Commenti
darova
il 22 Apr 2020
It appears that your function has 3 independent variable

How to plot this madness?
Neelanjan Pal
il 3 Mag 2020
darova
il 4 Mag 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

Risposte (0)
Categorie
Scopri di più su Calculus in Centro assistenza e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
