信号与系统作业.docx
《信号与系统作业.docx》由会员分享，可在线阅读，更多相关《信号与系统作业.docx（10页珍藏版）》请在冰豆网上搜索。
信号与系统作业
信号与系统实验
姓名:
***
学号：
***************
专业：
*************
实验一：
heaviside
functiony=heaviside（t）
y=zeros（size（t））;
y（t>0）=1;
y（t==0）=nan;
f=Heaviside（t）
plot（t,f）
axis（[-1,3,-0.2,1.2]）
实验二：
t=-5:
0.05:
5
f=sign（t）
subplot（2,1,1）
title（'符号函数sgn（t）'）
plot（t,f）,axis（[-5,5,-1.1,1.1]）
ff=1/2+1/2*f
subplot（2,1,2）
title（'单位阶跃信号u（t）'）
plot（t,ff）,axis（[-5,5,-1.1,1.1]）
实验三：
t=-3*pi:
pi/100:
3*pi;
ft=sinc（t/pi）;
plot（t,ft）
title（'³éÑùº¯Êýsa（t）'）
xlabel（'t'）
ylabel（'sa（t）'）
实验四：
k=0:
40;
subplot（2,1,1）
stem（k,cos（k*pi/8）,'filled'）
title（'cos（k*pi/8）'）
subplot（2,1,2）
stem（k,cos（2*k）,'filled'）
title（'cos（2*k）'）
实验五：
symst
f1=sym（'（t-4）*（heaviside（t）-heaviside（t-4））'）
subplot（2,3,1）,ezplot（f1）
f2=sym（'sin（2*pi*t）'）
subplot（2,3,4）,ezplot（f2,[-4,4]）
y1=subs（f1,t,-t）
f3=f1+y1
subplot（2,3,2）,ezplot（f3）
f4=-f3
subplot（2,3,3）,ezplot（f4）
f5=f2*f3
subplot（2,3,5）,ezplot（f5）
f6=f1*f2
subplot（2,3,6）,ezplot（f6）
实验六：
clc
clf
N=256;%时间分隔的点数
R=10/N;t=-2:
R:
2;
f=heaviside（t+1）-heaviside（t-1）;
wf=2*pi*50;%频率宽度
Nf=64;%频谱的点数
W=linspace（0,wf,Nf）;%W为频率正半轴的采样点
dw=wf/（Nf-1）;
F=f*exp（-j*t'*W）*R;%求F（jw）
F=real（F）;
W=[-fliplr（W）,W（2:
Nf）];%补上负频率
F=[fliplr（F）,F（2:
Nf）];%号补上负频率区的F（jw）的值
subplot（211）;plot（t,f,'linewidth',1.4）;
xlabel（'t'）;ylabel（'f（t）'）;
title（'f（t）=u（t+1）-u（t-1）'）;
subplot（212）;plot（W,F,'linewidth',1.4）;
title（'f（t）的傅里叶变换F（w）'）;
实验七：
symstwf
f=exp（-2*t）*sym（'Heaviside（t）'）;
F=fourier（f）
omiga=-10:
10;
subplot（3,1,1）;ezplot（f,[0:
2,0:
1.2]）;
title（'单边指数信号'）
subplot（3,1,2）;ezplot（abs（F）,[-10:
10]）;
title（'幅频特性图'）
subplot（3,1,3）;subplot（3,1,3）;plot（omiga,angle（subs（F,'w',omiga）））
title（'相频特性图'）
实验八：
symstw
F=1/（1+w^2）;
f=ifourier（F,w,t）
ezplot（f）
实验九：
a=[12-321];
b=[10-4];
p=roots（a）
q=roots（b）
holdon
plot（real（p）,imag（p）,'x'）;
plot（real（q）,imag（q）,'o'）;
title（'H（s）的零极点图'）;
gridon;
ylabel（'虚部'）;
xlabel（'实部'）
实验十：
a=[12-321];
b=[10-4];
sjdt（a,b）
实验十一：
a=[12321]
b=[10-4]
p=roots（a）;
pxm=max（real（p））;
ifpxm>=0
'系统不稳定'
else
freqs（b,a）
end
升级会员 