MATLAB实验三

  • 格式:doc
  • 大小:471.50 KB
  • 文档页数:30

下载文档原格式

  / 30
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

MATLAB实验三

五、课堂演示

1. 画正弦曲线sin(),[0,2]

=∈。

y x xπ

MATLAB代码:

x = 0:pi/100:2*pi;

y = sin(x);

plot(x,y)

axis([0 2*pi -1 1])

其中axis([0 2*pi -1 1])的作用是设置x轴坐标范围是从0到2π,y轴坐标是从-1到1。

效果图:

以上画图的效果也可以用hold on实现:x = 0:pi/100:2*pi;

y1 = sin(x);

y2 = sin(2*x);

y3 = sin(3*x);

plot(x,y1);

hold on;

plot(x,y2);

plot(x,y3);

hold off

title('y = exp(x)*sin(3x)')

%画第四幅子图

y = exp(x).*sin(4*x);

subplot(2,2,4);

plot(x,y)

title('y = exp(x)*sin(4x)')

注:title的作用是给图形写上标题效果图:

4. 作隐函数3330

+-=(笛卡尔儿叶形线)的图形

x y xy

MATLAB代码:

ezplot('x^3+y^3-3*x*y',[-2 2 -2 2])

axis equal

注: axis equal的作用是把纵横坐标的刻度比例一致效果图:

5. 画出空间螺旋线cos(),sin(),

===的图形

x t y t z t MATLAB代码:

t = 0:pi/100:10*pi;

x = cos(t);

y = sin(t);

z = t;

plot3(x,y,z)

title('x = cos(t),y = sin(t) z = t')

效果图:

6. 画曲面()

2222sin 2()x y z x y +=+的网线图

MATLAB 代码:

x = -pi:pi/50:pi;

y = x;

[X,Y] = meshgrid(x,y);

Z = sin(2*(X.^2+Y.^2))./(X.^2+Y.^2+eps);

mesh(X,Y,Z)

注:

(1) meshgrid产生一个以向量x为行、向量y为列的矩阵

(2) 分母加上eps是为了保证分母不为0

效果图:

bone ——兰色调灰色图

copper ——铜色色图

pink ——粉红色图

prism ——光谱色图

jet ——饱和值色图II

flag ——红、白、蓝交替色图

六、课堂练习

1. 在同一坐标系中画出函数

1sin(2*3)

y x

=+,2sin(32)

y x

=+的图形,且对于不以同的曲线使用不同的线型;

>> x=0:pi/100:2*pi;

y1=sin(2*x+3);

y2=sin(3*x+2);

plot(x,y1,'r-',x,y2,'g+')

legend('y1=sin(2*x+3)','y2=sin(3*x+2)')

>>

2. 当k分别取1, 2, 3, 4的时候,在同一个窗口的四个子图画出函数*sin(*)

y t k t

每个子图加上适当的标题;

>> t = 0:pi/100:2*pi;

% 画第一幅子图

plot(t,y)

title('y = t.*sin(4*t)')

>>

3. 画出三个两两相切的圆(假设这个三个圆的方程分别是:2222

x y

+-=);

(4)9

x y x y

1,(3)4

+=-+=,22

>> ezplot('x^2+y^2=1',[-7 7 -7 7]);

hold on

ezplot('(x-3)^2+y^2=4',[-7 7 -7 7]);

hold on

ezplot('x^2+(y-4)^2=9',[-7 7 -7 7]);

4. 在极坐标系下画出三叶形曲线3*sin(3*)

=,并尝试不同

的线形和颜色;

theta = 0:pi/100:2*pi;

>> r=3*sin(3*theta);

>> polar(theta,r,'r')

theta = 0:pi/100:2*pi; r=3*sin(3*theta); polar(theta,r,'g+')

5. 画出星形线33

==的图形;

x t y t

2*cos,2*sin

>> t = 0:pi/100:10*pi;

x=2*cos(t).^3;

>> y=2*sin(t).^3;

>> plot(x,y)

6. 画出双纽线22222

+=-的图形;

()3*()

x y x y

ezplot('(x^2+y^2)^2=3*(x^2-y^2)',[-3 3 -3 3])

7. 画出三维曲线图:*cos,*sin,

===

x t t y t t z t

>> t=0:pi/100:2*pi;