2.实验二 MATLAB绘图一 答案

实验二 MATLAB 绘图一

1. 编程绘制y=sin(t)/t 的曲线，t 的定义域是[-10Π,10Π]，绘图时加网格 解：t=[-10*pi:0.2:10*pi];y=sin(t)./t;plot(t,y),grid on

-40

-30-20-10010203040

-0.4-0.2

0.2

0.4

0.6

0.8

1

2. 在[0，10]之间用一张图画出y=sin(t)，y1=cos(t)的曲线，y 用红色实线绘制，

y1用蓝色长划线绘制，绘图时加网格，横纵坐标比例相同，横轴标明“时间”，纵轴标明“正弦、余弦”，图题“正弦和余弦曲线”，要有图例说明，且用鼠标拖动来标注“sin(t)”、“cos(t)”。 解： t=0:0.1:10;

y1=sin(t);y2=cos(t);plot(t,y1,'r',t,y2,'b--'); title('正弦和余弦曲线'); legend('正弦','余弦')

xlabel('时间t'),ylabel('正弦、余弦') grid axis square

gtext('sin(t)'),gtext('cos(t)')

2

4

6

8

10

-1-0.8-0.6-0.4-0.200.20.40.60.81正弦和余弦曲线

时间t

正弦、余弦

3. 用三种方法编程，同时在一张图上观察常用对数、自然对数函数在[0，10]

之间的曲线，其中在两种方法中，常用对数曲线用黑色实线绘制，自然对数曲线用红色“+”绘制，绘图时，MATLAB 不要提示“W arning ” 解一： t=[0.1:0.1:10]; y1=log10(t); y2=log(t);

plot(t,y1,'-k'),hold on plot(t,y2,'+r'),hold off

2

4

6

8

10

-2.5

-2-1.5-1-0.500.511.52

2.5

解二： t=[0.1:0.1:10]; y1=log10(t); y2=log(t); plot(t,[y1;y2])

2

4

6

8

10

-2.5

-2-1.5-1-0.500.511.52

2.5

解三： t=[0.1:0.1:10]; y1=log10(t); y2=log(t); plot(t,y1,'-k',t,y2,'+r')

2

4

6

8

10

-2.5

-2-1.5-1-0.500.511.52

2.5

4.曲线y=x+2x2+3x3，x的定义域为[-3，3]，在一张图上用排成一行的三幅子

图分别显示该曲线：黑色实线图、脉冲图、条形图，每幅图均有图题及横纵坐标轴说明

解：

x=[-3:0.1:3];

y=x+2*x.^2+3*x.^3;

subplot(1,3,1),plot(x,y,'k')

title('plot(x,y)')

xlabel('x'),ylabel('y')

subplot(1,3,2),stem(x,y)

title(' stem(x,y)')

xlabel('x'),ylabel('y')

subplot(1,3,3),bar(x,y)

title(' bar(x,y)')

xlabel('x'),ylabel('y')

5.通过MATLAB的help功能自学如何绘制饼图，在一张图上分上下两幅分别

绘制“通信08-1”、“通信08-2”、“电子08-1”、“电子08-2”的“MATLAB 大侠”比例为3：3：2：2的饼图和立体饼图，其中，“通信08-1”的饼被抽出。

解：

x=[3 3 2 2];y=[1 0 0 0];

subplot(2,1,1)

pie(x,y,{'通信08-1','通信08-2','电子08-1','电子08-2'})

subplot(2,1,2)

pie3(x,y,{'通信08-1','通信08-2','电子08-1','电子08-2'})

通

07-1

07-1

通信

6.尝试用几种方法画一个以原点为圆心，半径为1的圆。解：

方法一：

x=-1:0.01:1;

y=sqrt(1-x.^2);

y1=-sqrt(1-x.^2);

plot(x,y,x,y1)

axis equal

方法二：

t=0:0.1:10;

x=sin(t);

y=cos(t);

plot(x,y)

axis equal

方法三：

fplot('sqrt(1-x.^2)',[-1,1]);

hold on;

fplot('-sqrt(1-x.^2)',[-1,1])

axis equal

方法四：

fplot('[sqrt(1-x.^2),-sqrt(1-x.^2)]',[-1,1])

axis equal

方法五：

ezplot('cos(t)','sin(t)')

方法六：

ezplot ('x^2+y^2-1',[-1,1]) axis equal

方法七：

t=-2*pi:0.1:2*pi;

z=sin(t)+i*cos(t);

plot(z)

axis equal

方法八：

t=0:0.01:2*pi;

x=exp(i*t);

plot(x)

axis equal