数据可视化部分答案.doc

文档介绍：数据可视化(即画图)部分作业
已知椭圆的长、短轴,用“小红点线”画椭圆。(参见图p5-1)
程序:
t=2*pi*(0:199)/199;
a=4;b=2;
x=a*cos(t);y=b*sin(t);
plot(x,y,'r.','MarkerSize',15)
axis equal
图 p5-1
根据表达式绘制如图p5-2的心脏线。(提示:采用极坐标绘线指令polar)
clf
t = 0:.01:2*pi;
P=1-cos(t);
pline=polar(t,P,'r'),
set(pline,'LineWidth',5)
title('P=1-cos\theta')
图 p5-2
A,B,。试画出如图p5-3所示的三城市上半年每月生产总值的累计直方图。(提示:使用指令bar)
各城市生产总值数据(单位:亿元)
城市
1月
2月
3月
4月
5月
6月
A
170
120
180
200
190
220
B
120
100
110
180
170
180
C
70
50
80
100
95
120
clear
x=1:6;
Y=[170 120 180 200 190 220;
120 100 110 180 170 180;
70 50 80 100 95 120]
bar(x,Y','grouped');
bar(x,Y','stacked');
colormap(cool);
% legend('A','B','C','Location','NorthWest')
legend('A','B','C',2)
图 p5-3
二阶线性系统的归一化(即令)冲激响应可表示为:,其中,为阻尼系数。试在同一张图上,画出不同取值时的各条曲线,时间区间为。此外,要求(1)的各条曲线用蓝色,的用红色,的用黑色,且线较粗;(2)
和的两条曲线要有醒目的阻尼系数标志。(本题要求图形参见图p5-4)
%
clc,clf,clear;
t=(0::18)';
N=length(t);
zeta=::;
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<-2*eps
y=1/beta*exp(-zk*t).*sin(beta*t);
plot(t,y,'b')
if zk<
text(,,'\zeta = ')
end
elseif (zk-1)<2*eps%
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>
text(,,'\zeta = ')
end
end
e