-0.0007 0 0.0007 -0.0366 0 0.0366 -0.0007 0 0.0007 用 mesh 命令绘制上例中的网格曲面。 [X,Y] = meshgrid(-2:.2:2, -2:.2:2); Z = X .* exp(-X.^2 - Y.^2); mesh(Z)
上机练习（7）
x0.（ 1）按照的步长间隔 ? x ? 1 绘制函数 y ? xe ? 在0
≤x≤1时的曲线。
（ 2）用图形表示离散函数 y ? n ? 6 ? 1 ，其中 n 为[0,12]的自然数。 （ 3）分别采用
? x ?
10? 、 1 的步长，绘制连续
100?x2y21调制波形 y sin( t ) sin( 9 t ) 的图像。 ?（ 4）给出一系列的 a 值，采用函数 2 ? 2 ? 1 画
a25?a一组椭圆。 （
5）用曲面图命令
surf
表现函数 z ? x 2 ? y 2
的图像。
（ 6）绘制颜色为蓝色，数据点用五角星标识的下述函数在(0,5)上的虚线图。
y?xesinx用函数 pie 和 pie3 绘制饼图。 x=[200,360,120,400,320];
subplot(2,2,1), pie(x,[0 0 0 1 0])
subplot(2,2,2), pie3(x,[0 0 0 1 0])
subplot(2,2,3), pie(x(2:5))
subplot(2,2,4),
x=[0.1,0.12,0.21,0.34,0.11]; pie3(x ,{'A','B','C','D','E'})
绘制图形，并用函数 xlabel、title 和 legend 命令进行标注。t=0:0.1:4*pi; y=sin(t); y1=cos(t); plot(t,y,':',t,y1,'r*') xlabel('x 轴 (0--4\\pi)','fontsize',12,'fontweight','bold') ylabel('y 轴','fontsize',12,'fontweight','bold') title('
绘
制
正
弦
波
和
余
弦
Pos=1','fontsize',10,'fontweight','bold','fontangle','italic') text(pi,0,'\\leftarrowsin(\\pi)=0') text(pi,-1,'\\leftarrowcos(\\pi)=-1')
text(pi/2,0.9,['\%uparrowsin(\\pi/2)=',num2str(sin(pi/2))])
波
text(0,-0.6,['绘图日期：',date]) text(0,-0.8,['MATLAB 版本：',version]) legend('正弦波','余弦波') figure(2) plot(t,y,':',t,y1,'r*') title('
绘
制
正
弦
波
和
余
弦
波
Pos=0','fontsize',10,'fontweight','bold','fontangle','italic') legend('正弦波','余弦波',0) grid on figure(3) plot(t,y,':',t,y1,'r*') title('
绘
制
正
弦
波
和
余
弦
波
Pos=-1','fontsize',10,'fontweight','bold','fontangle','italic')
创建一个图形对象，练习和体会 set 和 get 函数的用法。 x=0:0.1:2*pi;
y1=sin(x); %创建一个正弦函数
H1_sin=plot(x,y1,'*') %画图并返回图像的句柄值 set(H1_sin,'color',[1,0.5,0],'Linewidth',1) %调用函数 set 设置正弦曲线的颜色和线宽 y2=cos(x); hold on
H1_cos=plot(x,y2)
set(H1_cos,'color',[0.25 0 1])
title('正弦和余弦函数曲线','fontsize',16,'color','blue')
H_sin_color=get(H1_sin,'color') %调用 get 函数返回具有句柄 H1_sin 的对象的颜色 程序执行结果如下 H1_sin = 3.0010 H1_cos = 101.0009 H_sin_color =
1.0000 0.5000 0
参数ωn ＝10， x0 ＝1， v0 ＝0，计算的终止时间 t=2。试求ξ0.1 到 1运动方程的解，并画出波形。 2. MATLAB 编程 编写 M 文件 ex1.m
%首先清空 MATLAB 的工作空间 clear; %给定初值 wn=10; tf=2; x0=1; v0=0;
从 %计算不同的ξ 值所对应的振型 for j=1:10; eta(j)=0.1*j; wd(j)=wn*sqrt(1-eta(j)^2); %求振幅 A
a=sqrt((wn*x0*eta(j)+v0)^2+(x0*wd(j))^2)/wd(j); %求相位角
phi=atan2(wd(j)*x0,v0+eta(j)*wn*x0); %设定自变量数组 t t=0:tf/1000:tf; %求过渡过程
x(j,:)=a*exp(-eta(j)*wn*t).*sin(wd(j)*t+phi); end
%在同一个图形窗口中绘制不同的ξ 值所对应的振型 plot(t,x(1,:),t,x(2,:),t,x(3,:),t,x(4,:),... t,x(5,:),t,x(6,:),t,x(7,:),t,x(8,:),... t,x(9,:),t,x(10,:)) grid on
%新建一个图形窗口，绘制三维网格图 figure mesh(x)