1、在直角坐标系下绘制(同一个窗口):
笛卡尔叶形线、星形线、摆线;
%在直角坐标系下绘制(同一个窗口):笛卡尔叶形线、星形线、摆线;clc
figure
%Descartes folium
theta_1=-2*pi:0.01:2*pi;%角度
t=tan(theta_1);
a=1;
x1=3*a*t./(1+t.^3);%参数方程
y1=3*a*t.^2./(1+t.^3);%参数方程
subplot(1,3,1);plot(x1,y1);
legend('笛卡尔叶形线');
axis([-4,4,-4,4]);%只显示局部
grid on;
%星形线
a=2;
theta=-2*pi:0.01:2*pi;
x2=a*cos(theta).^3;
y2=a*sin(theta).^3;
subplot(1,3,2);plot(x2,y2);
legend('星形线');
axis([-4,4,-4,4]);%只显示局部
grid on;
%摆线
a=2;
theta=-2*pi:0.001:2*pi;
x3=a.*(theta-sin(theta));
y3=a.*(1-cos(theta));
subplot(1,3,3);plot(x3,y3);
legend('摆线');
axis([-8,8,-8,8]);%只显示局部
grid on;
2、在极坐标系下绘制(加注释):
心形线,对数螺线、四叶玫瑰线
%在极坐标系下绘制(加注释):心形线
clc
figure
%心形线
a=2;
t=-2*pi:0.01:2*pi;
r=a.*(1+cos(t));
r=a.*(1+sin(t));
polar(t,r);
legend('心形线');
%在极坐标系下绘制(加注释):对数螺线clc
figure
%对数螺线
a=0.1;
t=-2*pi:0.001:2*pi;
r=exp(a*t);
polar(t,r);
legend('对数螺线');
%在极坐标系下绘制(加注释):四叶玫瑰线clc
figure
%四叶玫瑰线
a=4;
t=-2*pi:0.001:2*pi;
r=a*sin(2*t);
polar(t,r);
legend('四叶玫瑰线');
3、绘制双曲抛物面、单叶双曲面。
%双曲抛物面
clc
sqra=1;
sqrb=4;
[x,y]=meshgrid(-2:0.01:2);
z=(x.^2/sqra-y.^2/sqrb)/2;
plot3(x,y,z);
legend('双曲抛物面');
%单叶双曲面
clc
xi=-10:0.1:10;
yi=xi;
[x,y]=meshgrid(xi,yi);
a=1;b=1;c=1;
z=c^2*(x.^2/a^2+y.^2/b^2+1);
mesh(x,y,z);
hold on;
mesh(x,y,-z); legend('单叶双曲面');