MATLAB实验报告
班级：10自动化 姓名：王万飞 学号：18
1.在0≤x≤2区间内 y=2e-0.5xcos(4πx) 程序如下： >> x=0:pi/100:2*pi;
y=2*exp(-0.5*x).*cos(4*pi*x); plot(x,y)
2 绘制曲线。 程序如下： >> t=0:0.1:2*pi; x=t.*sin(3*t); y=t.*sin(t).*sin(t); plot(x,y);
3 分析下列程序绘制的曲线 >> x1=linspace(0,2*pi,100); x2=linspace(0,3*pi,100); x3=linspace(0,4*pi,100); y1=sin(x1); y2=1+sin(x2); y3=2+sin(x3);
x=[x1;x2;x3]'; y=[y1;y2;y3]'; plot(x,y,x1,y1-1)
4 用不同标度在同一坐标内绘制曲线y1=0.2e-0.5xcos(4πx) 和y2=2e-0.5xcos(πx)。 程序如下： >> x=0:pi/100:2*pi;
y1=0.2*exp(-0.5*x).*cos(4*pi*x); y2=2*exp(-0.5*x).*cos(pi*x); plotyy(x,y1,x,y2);
5 采用图形保持，在同一坐标内绘制曲线y1=0.2e-0.5xcos(4πx) 和y2=2e-0.5xcos(πx)。 程序如下： >> x=0:pi/100:2*pi;
y1=0.2*exp(-0.5*x).*cos(4*pi*x); plot(x,y1) hold on
y2=2*exp(-0.5*x).*cos(pi*x); plot(x,y2); hold off
6 在同一坐标内，分别用不同线型和颜色绘制曲线y1=0.2e-0.5xcos(4πx) 和y2=2e-0.5xcos(πx)，标记两曲线交叉点。 程序如下：
>> x=linspace(0,2*pi,1000); y1=0.2*exp(-0.5*x).*cos(4*pi*x); y2=2*exp(-0.5*x).*cos(pi*x); k=find(abs(y1-y2)<1e-2);
x1=x(k); y3=0.2*exp(-0.5*x1).*cos(4*pi*x1); plot(x,y1,x,y2,'k:',x1,y3,'bp');
7 在0≤x≤2区间内，绘制曲线y1=2e-0.5x和y2=cos(4πx)，并给图形添加图形标注。 程序如下： >> x=0:pi/100:2*pi; y1=2*exp(-0.5*x); y2=cos(4*pi*x); plot(x,y1,x,y2)
title('x from 0 to 2{\\pi}'); %加图形标题 xlabel('Variable X'); %加X轴说明 ylabel('Variable Y'); %加Y轴说明
text(0.8,1.5,'曲线y1=2e^{-0.5x}'); %在指定位置添加图形说明 text(2.5,1.1,'曲线y2=cos(4{\\pi}x)');
legend(‘y1’,‘ y2’) %加图例
??? legend(‘y1’,‘ y2’) %加图例 |
Error: The input character is not valid in MATLAB statements or expressions.
8 在同一坐标中，可以绘制3个同心圆，并加坐标控制。 程序如下： >> t=0:0.01:2*pi; x=exp(i*t); y=[x;2*x;3*x]'; plot(y)
grid on; %加网格线 box on; %加坐标边框
axis equal ；%坐标轴采用等刻度
??? Error using ==> get Invalid handle object.
Error in ==> plotyy>localUpdatePosition at 378 newDestPos =
hgconvertunits(hFig,newPos,get(axSource,'Units'),get(axDest,'Units'),get(axSource,'Parent'));
Warning: Error occurred while evaluating listener callback. > In axis>LocSetEqual at 361 In axis at 135 ??? Error using ==> get Invalid handle object.
Error in ==> plotyy>localUpdatePosition at 378 newDestPos =
hgconvertunits(hFig,newPos,get(axSource,'Units'),get(axDest,'Units'),get(axSource,'Parent'));
Warning: Error occurred while evaluating listener callback. > In axis>LocSetEqual at 362 In axis at 135
??? Error using ==> axis at 183 Unknown command option ；
9 用fplot函数绘制f(x)=cos(tan(πx))的曲线。 命令如下：
>> fplot('cos(tan(pi*x))',[ 0,1],1e-4)
10.分别以条形图、阶梯图、杆图和填充图形式绘制曲线y=2sin(x)。
程序如下： >> x=0:pi/10:2*pi; y=2*sin(x);
subplot(2,2,1);bar(x,y,'g'); title('bar(x,y,''g'')');axis([0,7,-2,2]); subplot(2,2,2);stairs(x,y,'b'); title('stairs(x,y,''b'')');axis([0,7,-2,2]); subplot(2,2,3);stem(x,y,'k'); title('stem(x,y,''k'')');axis([0,7,-2,2]); subplot(2,2,4);fill(x,y,'y'); title('fill(x,y,''y'')');axis([0,7,-2,2]);
11.
隐函数绘图应用举例。
程序如下：
>> subplot(2,2,1);
ezplot('x^2+y^2-9');axis equal subplot(2,2,2);
ezplot('x^3+y^3-5*x*y+1/5') subplot(2,2,3);
ezplot('cos(tan(pi*x))',[ 0,1]) subplot(2,2,4);
ezplot('8*cos(t)','4*sqrt(2)*sin(t)',[0,2*pi])
12绘制了peaks函数曲面并且将它绕z轴旋转。 程序如下
>> [X,Y,Z]=peaks(30); surf(X,Y,Z)
axis([-3,3,-3,3,-10,10])
axis off; shading interp; colormap(hot);
m=moviein(20); for i=1:20
view(-37.5+24*(i-1),30) m(:,i)=getframe; end
movie(m,2);
13.光照处理后的球面。 程序如下：
>> [x,y,z]=sphere(20); subplot(1,2,1); surf(x,y,z);axis equal; light('Posi',[0,1,1]);
shading interp; hold on;
plot3(0,1,1,'p');text(0,1,1,' light'); subplot(1,2,2); surf(x,y,z);axis equal; light('Posi',[1,0,1]); shading interp; hold on;
plot3(1,0,1,'p');text(1,0,1,' light');
