Please help in plotting the graphs
Show older comments
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 Comments
darova
on 22 Apr 2020
It appears that your function has 3 independent variable

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

Answers (0)
Categories
Find more on Calculus in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!