5.5 图形窗功能简介
图 5.5-1 图形窗工具条专用按键
【例5.5-1】
图5.5-2 个性化的图形
(1)
t=(0:pi/100:pi)'; y1=sin(t)*[1,-1]; y2=sin(t).*sin(9*t); t3=pi*(0:9)/9; y3=zeros(size(t3)); plot(t,y1,t,y2,t3,y3)
(2)
21
图 5.5-3 编辑工作模式下的图形窗界面
(3)
图 5.5-4 横坐标上限设置示例
(4)
图 5.5-5 包络线的设置示例
(5)
图 5.5-6 调制曲线的设置示例
22
(6)
图 5.5-7 调制曲线的设置示例
习题5
?x?acost1. 已知椭圆的长、短轴a?4,b?2,用“小红点线”画如图p5-1所示的椭圆?。
y?bsint?(提示:参量t;点的大小;axis equal) 3210y-1-2-3-4-3-2-10x1234图 p5-1
2. 根据表达式??1?cos?绘制如图p5-2的心脏线。(提示:polar;注意title中特殊字
符;线宽;axis square。可以用plot试试。) ?=1-cos?12090 2 1.5150 1 0.518003060210330240270300图 p5-2 23
3. A,B,C三个城市上半年每个月的国民生产总值表p5.1。试画出如图p5-3所示的三城市上
半年每月生产总值的累计直方图。(提示:bar(x,Y,'style'); colormap(cool); legend。)
表p5.1 各城市生产总值数据(单位:亿元) 城市 A B C 1月 170 120 70 6005004003002001000 ABC2月 120 100 50 3月 180 110 80 4月 200 180 100 5月 190 170 95 6月 220 180 120 123456图 p5-3
4. 二阶线性系统的归一化(即令?n?1)冲激响应可表示为:
1?? t?esin(? t)???y(t)??te?t?1e?(???)t?e?(???)t?2??0???1????1,其中??|1??2|,?为阻尼??1系数。(1)希望在同一张图上,绘制t?[0,18]区间内??0.2:0.2:1.4不同取值时的各条曲线(参见图p5-4)。在此图上,??1的各条曲线为细蓝线;??1为粗黑线;??1为细红线;并且对最上方及最下方的两条曲线给出??0.2和??1.4的醒目标志。(2)读者运行题下程序exmp504.m,可以发现该程序画出的曲线中没有“粗黑线”。你能讲出原因吗?如何对exmp504.m作最少的修改(比如只改一条指令),就可画出所需图形。(提示:该题深层次地暴露数值计算可能存在的隐患。)
24
0.8??=0.20.6? = 0.20.40.2? = 1.40-0.2-0.4024681012141618图 p5-4
% exmp504.m 供第4道习题使用的程序 clc,clf,clear; t=(0:0.05:18)'; N=length(t);
zeta=0.2:0.2:1.4; L=length(zeta); y=zeros(N,L); hold on for k=1:L
zk=zeta(k);
beta=sqrt(abs(1-zk^2));
if zk<1 %满足此条件,绘蓝色线 y=1/beta*exp(-zk*t).*sin(beta*t); plot(t,y,'b') if zk<0.4
text(2.2,0.63,'\\zeta = 0.2') end
elseif zk==1 %满足此条件,绘黑色线 y=t.*exp(-t);
plot(t,y,'k','LineWidth',2)
else %其余,绘红色线
y=(exp(-(zk-beta)*t)-exp(-(zk+beta)*t))/(2*beta); plot(t,y,'r') if zk>1.2
text(0.3,0.14,'\\zeta = 1.4') end end end
text(10,0.7,'\\Delta\\zeta=0.2') axis([0,18,-0.4,0.8]) hold off box on grid on
5. 用绿实线绘制x?sin(t),y?cos(t),z?t的三维曲线,曲线如图p5-5所示。(提
示:参变量;plot3;线色线粗。)
25