山东科技大学数字信号处理实验五Word下载.docx

山东科技大学数字信号处理实验五Word下载.docx

《山东科技大学数字信号处理实验五Word下载.docx》由会员分享，可在线阅读，更多相关《山东科技大学数字信号处理实验五Word下载.docx（14页珍藏版）》请在冰豆网上搜索。

.'

）;

xlabel（'

n'

ylabel（'

x1（n）'

n=0:

10;

x1=（n+1）.*（n>

=0&

n<

=3）;

x2=（8-n）.*（n>

=4&

=7）;

x=x1+x2;

x2（n）'

x1=（4-n）.*（n>

x2=（n-3）.*（n>

x3（n）'

n=-5:

0.2:

5;

x=cos（pi*n/4）;

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）;

t'

x5（n）'

（4）8点

7;

X=fft（x,8）;

subplot（2,1,1）;

stem（n,X,）;

X1（n）'

subplot（2,1,2）;

stem（n,abs（X））;

）

|X1（e^j^\omega）|'

title（'

幅频特性'

16点

15;

X=fft（x,16）;

stem（n,X）;

32点：

31;

X=fft（x,32）;

（5）

X2

8点

X2（n）'

|X2（e^j^\omega）|'

X3

X3（n）'

|X3（e^j^\omega）|'

（6）

4点

3;

X=fft（x,4）;

X4（n）'

|X4（e^j^\omega）|'

（7）

fs=64;

N5=16;

n5=0:

N5-1;

k5=0:

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）;

N5=32;

N5=64;

可先截取M点进行DFT,再将截取长度扩大1倍截取,比较结果,如果二者的差别满足分析误差要求,则可以近似表示该信号的频谱,如果不满足误差要求就继续将截取长度加倍,重复比较,直到结果满足要求