山东科技大学数字信号处理实验五.docx
《山东科技大学数字信号处理实验五.docx》由会员分享，可在线阅读，更多相关《山东科技大学数字信号处理实验五.docx（14页珍藏版）》请在冰豆网上搜索。
山东科技大学数字信号处理实验五
（3）
n=-1:
4;
x=[（n>=0）&（n<=3）];
stem（n,x,'.'）;
xlabel（'n'）;
ylabel（'x1（n）'）;
n=0:
10;
x1=（n+1）.*（n>=0&n<=3）;
x2=（8-n）.*（n>=4&n<=7）;
x=x1+x2;
stem（n,x,'.'）;
xlabel（'n'）;
ylabel（'x2（n）'）;
n=0:
10;
x1=（4-n）.*（n>=0&n<=3）;
x2=（n-3）.*（n>=4&n<=7）;
x=x1+x2;
stem（n,x,'.'）;
xlabel（'n'）;
ylabel（'x3（n）'）;
n=-5:
0.2:
5;
x=cos（pi*n/4）;
stem（n,x,'.'）;
xlabel（'n'）;
ylabel（'x4（n）'）;
t=-0.4:
0.0001:
0.4;
x1=cos（8*pi*t）;
x2=cos（16*pi*t）;
x3=cos（20*pi*t）;
x=x1+x2+x3;
plot（t,x）;
xlabel（'t'）;
ylabel（'x5（n）'）;
（4）8点
n=0:
7;
x=[（n>=0）&（n<=3）];
X=fft（x,8）;
subplot（2,1,1）;
stem（n,X,）;
xlabel（'n'）;
ylabel（'X1（n）'）;
subplot（2,1,2）;
stem（n,abs（X））;
xlabel（'n'）
ylabel（'|X1（e^j^\omega）|'）
title（'幅频特性'）;
16点
n=0:
15;
x=[（n>=0）&（n<=3）];
X=fft（x,16）;
subplot（2,1,1）;
stem（n,X）;
xlabel（'n'）;
ylabel（'X1（n）'）;
subplot（2,1,2）;
stem（n,abs（X））;
xlabel（'n'）
ylabel（'|X1（e^j^\omega）|'）
title（'幅频特性'）;
32点：
n=0:
31;
x=[（n>=0）&（n<=3）];
X=fft（x,32）;
subplot（2,1,1）;
stem（n,X）;
xlabel（'n'）;
ylabel（'X1（n）'）;
subplot（2,1,2）;
stem（n,abs（X））;
xlabel（'n'）
ylabel（'|X1（e^j^\omega）|'）
title（'幅频特性'）;
（5）
X2
8点
n=0:
7;
x1=（n+1）.*（n>=0&n<=3）;
x2=（8-n）.*（n>=4&n<=7）;
x=x1+x2;
X=fft（x,8）;
subplot（2,1,1）;
stem（n,X）;
xlabel（'n'）;
ylabel（'X2（n）'）;
subplot（2,1,2）;
stem（n,abs（X））;
xlabel（'n'）
ylabel（'|X2（e^j^\omega）|'）
title（'幅频特性'）;
X2
16点
n=0:
15;
x1=（n+1）.*（n>=0&n<=3）;
x2=（8-n）.*（n>=4&n<=7）;
x=x1+x2;
X=fft（x,16）;
subplot（2,1,1）;
stem（n,X）;
xlabel（'n'）;
ylabel（'X2（n）'）;
subplot（2,1,2）;
stem（n,abs（X））;
xlabel（'n'）
ylabel（'|X2（e^j^\omega）|'）
title（'幅频特性'）;
X3
8点
n=0:
7;
x1=（4-n）.*（n>=0&n<=3）;
x2=（n-3）.*（n>=4&n<=7）;
x=x1+x2;
X=fft（x,8）;
subplot（2,1,1）;
stem（n,X）;
xlabel（'n'）;
ylabel（'X3（n）'）;
subplot（2,1,2）;
stem（n,abs（X））;
xlabel（'n'）
ylabel（'|X3（e^j^\omega）|'）
title（'幅频特性'）;
X3
16点
n=0:
15;
x1=（4-n）.*（n>=0&n<=3）;
x2=（n-3）.*（n>=4&n<=7）;
x=x1+x2;
X=fft（x,16）;
subplot（2,1,1）;
stem（n,X）;
xlabel（'n'）;
ylabel（'X3（n）'）;
subplot（2,1,2）;
stem（n,abs（X））;
xlabel（'n'）
ylabel（'|X3（e^j^\omega）|'）
title（'幅频特性'）;
（6）
4点
n=0:
3;
x=cos（pi*n/4）;
X=fft（x,4）;
subplot（2,1,1）;
stem（n,X）;
xlabel（'n'）;
ylabel（'X4（n）'）;
subplot（2,1,2）;
stem（n,abs（X））;
xlabel（'n'）
ylabel（'|X4（e^j^\omega）|'）
title（'幅频特性'）;
8点
n=0:
7;
x=cos（pi*n/4）;
X=fft（x,8）;
subplot（2,1,1）;
stem（n,X）;
xlabel（'n'）;
ylabel（'X4（n）'）;
subplot（2,1,2）;
stem（n,abs（X））;
xlabel（'n'）
ylabel（'|X4（e^j^\omega）|'）
title（'幅频特性'）;
16点
n=0:
15;
x=cos（pi*n/4）;
X=fft（x,16）;
subplot（2,1,1）;
stem（n,X）;
xlabel（'n'）;
ylabel（'X4（n）'）;
subplot（2,1,2）;
stem（n,abs（X））;
xlabel（'n'）
ylabel（'|X4（e^j^\omega）|'）
title（'幅频特性'）;
（7）
fs=64;
N5=16;
n5=0:
N5-1;
k5=0:
N5-1;
T=1/fs;
t=n5*T;
x5=cos（8*pi*t）+cos（16*pi*t）+cos（20*pi*t）;
X5k=fft（x5,N5）
fudu5=abs（x5）
subplot（211）;
stem（n5,fudu5）;
subplot（212）;
stem（k5,X5k）;
fs=64;
N5=32;
n5=0:
N5-1;
k5=0:
N5-1;
T=1/fs;
t=n5*T;
x5=cos（8*pi*t）+cos（16*pi*t）+cos（20*pi*t）;
X5k=fft（x5,N5）
fudu5=abs（x5）
subplot（211）;
stem（n5,fudu5）;
subplot（212）;
stem（k5,X5k）;
fs=64;
N5=64;
n5=0:
N5-1;
k5=0:
N5-1;
T=1/fs;
t=n5*T;
x5=cos（8*pi*t）+cos（16*pi*t）+cos（20*pi*t）;
X5k=fft（x5,N5）
fudu5=abs（x5）
subplot（211）;
stem（n5,fudu5）;
subplot（212）;
stem（k5,X5k）;
可先截取M点进行DFT,再将截取长度扩大1倍截取,比较结果,如果二者的差别满足分析误差要求,则可以近似表示该信号的频谱,如果不满足误差要求就继续将截取长度加倍,重复比较,直到结果满足要求
升级会员 