文档介绍:数据可视化(即画图)部分作业
x二acost
、短轴a=4,b=2,用“小红点线”画椭圆<7.。(参见图p5-l)
[y二bsint
程序:
t=2*pi*(0:199)/199;
a=4;b=2;
x=a*c数据可视化(即画图)部分作业
x二acost
、短轴a=4,b=2,用“小红点线”画椭圆<7.。(参见图p5-l)
[y二bsint
程序:
t=2*pi*(0:199)/199;
a=4;b=2;
x=a*cos(t);y=b*sin(t);
plot(x,y,'r・','MarkerSize',15)
axisequal
-3-2-101234
图p5-1
=1-cos9绘制如图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
,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=[170120180200190220;
1201001**********;
70508010095120]
bar(x,Y','grouped');
bar(x,Y','stacked');colormap(cool);
%legend('A','B','C','Location','NorthWest')legend('A','B','C',2)
图p5-3
(即令®二1)冲激响应可表示为:
n
y(t)=<
1
2p
e-匚tsin(Bt)p
te-1
-(匚一卩)t——e-(匚+卩)t]
o<C<1
C=1,其中Phj1—匚21,:为阻尼匚>1
系数。试在同一张图上,画出C=::,时间区间为
te[0,18]。此外,要求(1)匚<1的各条曲线用蓝色,匚>1的用红色,匚二1的用黑
色,且线较粗;(2)C==。(本题
要求图形参见图p5-4)
%,clf,clear;t=(0::18)';
N=length(t);zeta=::;
L=length(zeta);y=zeros(N,L);
holdon
fork=1:L
zk=zeta(k);
beta=sqrt(abs(1-zkA2));
ifzk-1<-2*eps
y=1/beta*exp(-zk*t).*sin(