clear;
clc;
%%
x=1;
n=0;
for phi=0:3:360
n=n+1;
phi1=phi/180*pi;
y1(n)=((sin(pi*cos(phi1)*x)/(pi*cos(phi1)*x))^2*(cos(phi1))^2+(sin(pi*sin(phi1)*x)/(pi*sin(phi1)*x))^2*(sin(phi1))^2)^0.5;
end
%%
x=1/2;
n=0;
for phi=0:3:360
n=n+1;
phi1=phi/180*pi;
y2(n)=((sin(pi*cos(phi1)*x)/(pi*cos(phi1)*x))^2*(cos(phi1))^2+(sin(pi*sin(phi1)*x)/(pi*sin(phi1)*x))^2*(sin(phi1))^2)^0.5;
end
x=1/3;
n=0;
for phi=0:3:360
n=n+1;
phi1=phi/180*pi;
y3(n)=((sin(pi*cos(phi1)*x)/(pi*cos(phi1)*x))^2*(cos(phi1))^2+(sin(pi*sin(phi1)*x)/(pi*sin(phi1)*x))^2*(sin(phi1))^2)^0.5;
end
x=1/4;
n=0;
for phi=0:3:360
n=n+1;
phi1=phi/180*pi;
y4(n)=((sin(pi*cos(phi1)*x)/(pi*cos(phi1)*x))^2*(cos(phi1))^2+(sin(pi*sin(phi1)*x)/(pi*sin(phi1)*x))^2*(sin(phi1))^2)^0.5;
end
x=1/5;
n=0;
for phi=0:3:360
n=n+1;
phi1=phi/180*pi;
y5(n)=((sin(pi*cos(phi1)*x)/(pi*cos(phi1)*x))^2*(cos(phi1))^2+(sin(pi*sin(phi1)*x)/(pi*sin(phi1)*x))^2*(sin(phi1))^2)^0.5;
end
x=1/8;
n=0;
for phi=0:3:360
n=n+1;
phi1=phi/180*pi;
y8(n)=((sin(pi*cos(phi1)*x)/(pi*cos(phi1)*x))^2*(cos(phi1))^2+(sin(pi*sin(phi1)*x)/(pi*sin(phi1)*x))^2*(sin(phi1))^2)^0.5;
end
x=1/14;
n=0;
for phi=0:3:360
n=n+1;
phi1=phi/180*pi;
y14(n)=((sin(pi*cos(phi1)*x)/(pi*cos(phi1)*x))^2*(cos(phi1))^2+(sin(pi*sin(phi1)*x)/(pi*sin(phi1)*x))^2*(sin(phi1))^2)^0.5;
end
x=1/20;
n=0;
for phi=0:3:360
n=n+1;
phi1=phi/180*pi;
y20(n)=((sin(pi*cos(phi1)*x)/(pi*cos(phi1)*x))^2*(cos(phi1))^2+(sin(pi*sin(phi1)*x)/(pi*sin(phi1)*x))^2*(sin(phi1))^2)^0.5;
end
figure(100)
plot(0:3:360,y2,'-*k','linewidth',2)
hold on
plot(0:3:360,y3,'-ok','linewidth',2)
plot(0:3:360,y4,'-^k','linewidth',2)
plot(0:3:360,y5,'.-k','linewidth',2)
plot(0:3:360,y8,'--k','linewidth',2)
plot(0:3:360,y14,'-.k','linewidth',2)
plot(0:3:360,y20,'.k','linewidth',2)
legend('2','3','4','5','8','14','20')
hold off
grid on
phi=0:3:360;
phi1=phi/180*pi;
% figure(101)
% polar(phi1,y2(phi),'-*k','linewidth',2)
% hold on
% polar(phi1,y3(phi),'-ok','linewidth',2)
% polar(phi1,y4(phi),'-^k','linewidth',2)
% polar(phi1,y5(phi),'.-k','linewidth',2)
% polar(phi1,y8(phi),'--k','linewidth',2)
% polar(phi1,y14(phi),'-.k','linewidth',2)
% polar(phi1,y20(phi),'.k','linewidth',2)
% legend('2','3','4','5','8','14','20')
% hold off
% grid on
评论0