用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('悬链线')