用MatLab制作的几个数学函数图像

文字加注：
x=-1.5:0.001:1.5;
y=(x.^2-1).^3+1;
plot(x,y)
title('\fontsize{14}\fontname{宋体}函数图像：y=(x^2-1)^3+1')
xlabel('\fontsize{14}x'),ylabel('\fontsize{14}y')
text(-1,1.1,'\fontsize{8}点（1，1）处倒数为零，但无极值')
x=-10:1:10;
y=-(x-5).^2+2;
[y_max,x_max]=max(y);
num2str(y_max);
num2str(x_max);
plot(x,y)
hold on
plot(y_max,t_max,'r.')
hold off
字符串的应用:
a=2;
w=3;
t=0:0.01:10;
y=exp(-a*t).*sin(w*t);
[y_max,t_max]=max(y);
t_text=['t=',num2str(t (t_max))];
y_text=['y=',num2str(y_max)];
max_text=char('maxinum',t_text,y_text);
tit=['字符串的应用:y=exp(-',num2str(a),'t)*sin(',num2str(w),'t)']; hold on
plot(t,y,'b')
plot(t(t_max),y_max,'r.')%最大值处以红点标示
text(t(t_max)+0.3,y_max+0.05,max_text)
title(tit),xlabel('t'),ylabel('y')
hold off
求近似极限，修补图形缺口：
t=-2*pi:pi/10:2*pi;
y=sin(t)./t;
tt=t+(t==0)*eps;%逻辑数组参与运算，用“机械零”代替零元素
yy=sin(tt)./tt;%用数值可算的sin(eps)/eps近似替代sin(0)/0
subplot(1,2,1),plot(t,y),title('残缺图形
'),xlabel('t'),ylabel('y'),axis([-7,7,-0.5,1.2])
subplot(1,2,2),plot(tt,yy),title('正确图形
'),xlabel('tt'),ylabel('yy'),axis([-7,7,-0.5,1.2])
三维图像：
z=peaks(50);
[x,y]=meshgrid([-4:0.1:4]);
m=peaks(x,y);
subplot(2,2,1);h=mesh(z);
subplot(2,2,2);surf(z);
subplot(2,2,3);plot3(x,y,m)
摆线方程图形：
a=1;
t=0:0.1:4*pi;
x=a*(t-sin(t));
y=a*(1-cos(t));
plot(x,y)
axis image
title('\fontsize{14}\fontname{宋体}摆线方程图形')
xlabel('x')
ylabel('y')
箕形线图形:
a=1;
x=-8:0.1:8;
y=8*a.^3./(x.^2+3*a.^2);
plot(x,y)
axis image
title('\fontsize{14}\fontname{宋体}箕形线图形');
xlabel('\fontsize{14}x'),ylabel('\fontsize{14}y')
阿基米德螺旋线
a=2;
t=0:0.1:5*pi;
m=a*t;
x=m.*cos(t);
y=m.*sin(t);
plot(x,y)
hold on
plot(-x,y,'b:')
axis image
title('\fontsize{14}\fontname{宋体}阿基米德螺旋线:m=at') xlabel('\fontsize{14}x')
ylabel('\fontsize{14}y')
四叶玫瑰线：
a=2;
t=0:0.01:2*pi;
b=a*sin(2*t);
x=b.*cos(t);
y=b.*sin(t);
plot(x,y)
axis image
title('\fontsize{14}\fontname{仿宋}四叶玫瑰线：b=a*sin(2*t)') xlabel('x')
ylabel('y')
对数螺线：
a=2;
t=-pi:0.01:pi;
b=a.*exp(t);
x=b.*cos(t);
y=b.*sin(t);
plot(x,y)
axis image
title('\fontsize{14}\fontname{仿宋}对数螺线：b=a*e^t')
xlabel('x')
ylabel('y')
渐伸线：
a=5;
t=0:0.01:1.5*pi;
x=a*(cos(t)+t.*sin(t));
y=a*(sin(t)-t.*cos(t));
plot(x,y)
hold on
t1=0:0.01:2*pi;
plot(a.*cos(t1),a.*sin(t1))
hold off
title('渐开线')
xlabel('x'),ylabel('y')
悬链线方程：
x=-10:0.1:20;
a=5;
y=a.*(exp(x./a)+exp(-x./a))./2;
plot(x,y)
title('悬链线')。