Please help in plotting the graphs

Neelanjan Pal

2020 4 22

0 个回答

13 次查看（30 天）

0 个投票

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 个评论
显示 1更早的评论隐藏 1更早的评论

Neelanjan Pal 2020-5-3

在 MATLAB Online 中打开

The above graph I got by executing the following function containing 2 variables which is not proper. Please give any suggestions.

fun1=@(l,ki) ((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)

darova 2020-5-4

在 MATLAB Online 中打开

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