Matlab 代码
a = 2.09;
b = 1.31;
c = 1.16;
d = 0.885;
m = 0.036;
n = 0.615;
p = 35.237;
syms alpha;
alpha = -pi/2:0.01:pi/2;
x = a*exp(c*alpha).*cos(c*alpha) + exp(-(alpha-n).*(alpha-n)/(2*m*m))./(sqrt(2*pi*m*m)*p);
y = b*exp(b*alpha).*sin(d*alpha);
plot(x,y);hold on;
axis([0 4 -2 8]);
效果图:
matlab 代码
syms alpha;
alpha = 0:0.01:2*pi;
x = 13*power(sin(alpha),4);
y = -17*cos(alpha) + 9*cos(2*alpha) -0.1*cos(3*alpha) + cos(4*alpha) +0.1*cos(5*alpha);
plot(x,y);hold on;
axis([0 14 -16 28]);
效果图: