MATLAB教程2012a第5章习题解答-yxs分析
数据和函数可视化
演习问题5和答案
知道椭圆的长、短轴，用“小红虚线”画椭圆。 (参照图p5-1 )
)解答
clf
a=4； b=2；
t=0:pi/80:2*pi；
x=a*cos(t；
y=b*sin(t )；
打印(x，y，' r.'，' MarkerSize '，15 ) )。
axis equal
xlabel(x ) ) )。
ylabel(y ) )。
shg
图p5-1
根据公式绘制图p5-2的心脏线。 (提示：使用极坐标绘图命令polar )
)解答
clf
theta=0:pi/50:2*pi；
rho=1- cos (theta )；
h=polar(Theta，rho，'-r ' )； %极坐标绘图命令。 h是画的线条图案。
使用set(h，' LineWidth '，4 ) set设置h图形对象的“线宽”
axis square%保证坐标的舍入性
title'rho=1-costheta'%使用特殊字符映射
图p5-2
a、b、c三个城市上半年月度国民生产总值见表p5.1。 试着制作如图p5-3所示的三个城市上半年每月生产总值的累计直方图。
表p5.1各城市生产总值数据(单位：亿元) )。
城市1月2月3月4月5月6月a 170120180200190220 b 120100110180170180 c 70508010095120
目的，目的
使用MATLAB帮助系统学习使用直方图命令polar。
bar命令的常见格式之一： bar(x，y，' style ' )。 x是自变量向量； y是与x行数相同的矩阵，y的各行为“一组”数据； style取stacked时，与同一组数据中的各要素对应的长方体相互重叠。
在本例中，y列表示城市。 制图时，各列颜色取自《色图》； 图形窗口的默认颜色图是jet。
)解答
x=(1:6 )；
y=[170、120、180、200、190、220； 120、100、110、180、170、180； 70、50、80、100、95、120 (；
%以下两个命令是为了向读者展示x、y中对应的数据关系而编写的。 不影响画画。
disp ([空白(5)、(x )、空白(4)、(YA )、空白(4)、YB )、空白(4)、(YC ) ]
磁盘([ x，Y] ) )。
%
bar(x，y，'堆叠' )； %层叠直条
颜色映射(颜色)； 采用cool色图
Egend(a )、b )、c )、2 )； %此命令的第四个输入量用于控制图例的位置。
axis ([ 0，7，0，600 ]
shg
x YA YB YC
1 170 120 70
2 120 100 50
3 180 110 80
4 200 180 100
5 190 170 95
6 220 180 120
图p5-3
二阶线性系统的归一化(即令)冲激响应可以表示为：
其中是衰减系数。 )1)在同一图上，想画出区间内取不同值时的各曲线(参照图p5-4 )。 在此图中，的每条曲线都是一条细线。 粗黑线； 细红线； 然后，在最上面和最下面的两条曲线上赋予和的明显标记。 )读者运行问题程序exmp504.m后，会发现该程序绘制的曲线没有“粗黑线”。 能告诉我理由吗？ 通过对exmp504.m进行最小限度的修改，例如可以通过仅更改一个命令来绘制所需的图元。 提示：这个问题深刻暴露了数值计算中可能存在的危险性。 ）
)解答
clc、clf、clear；
t=(0:0.05:18 )； n=长度(t；
泽塔=0.2:0.2:1.4； 可能画不出%黑线。 3
%zeta=linspace(0.2、1.4、7 )；
L=Length)；
y=Zeros(n，l )；
霍尔德on
for k=1:L
ZK=zeta (k；
贝塔=sqrt(ABS(1-ZK^2)；
if zk1
y=1/beta*exp(-ZK*t ).*sin ) beta*t )；
打印(t，y，' b ' ) )。
if zk0.4
text ) 2.2，0.63，'zeta=0.2 ' )
结束
elseif zk==1
y=t.*exp(-t )；
打印(t，y，' k '，'线宽度'，2 ) )。
else