哈尔滨工程大学数字信号处理五Word文档格式.docx

哈尔滨工程大学数字信号处理五Word文档格式.docx

《哈尔滨工程大学数字信号处理五Word文档格式.docx》由会员分享，可在线阅读，更多相关《哈尔滨工程大学数字信号处理五Word文档格式.docx（24页珍藏版）》请在冰点文库上搜索。

（3）函数hamming（）用于产生汉明窗，调用格式为：

w=hamming（N）

（4）函数bartlett（）用于产生巴特利窗，调用格式为：

w=bartlett（N）

（5）函数blackman（）用于产生布莱克曼窗，调用格式为：

w=blackman（N）

（6）函数traing（）用于产生traing窗，调用格式为：

w=traing（N）

（7）函数kaiser（）用于产生kaiser窗，调用格式为：

w=kaiser（N）

（8）函数chebwin（）用于产生切比雪夫窗，调用格式为：

w=chebwin（N）

二、实验内容

1.用MATLAB编程绘制各种窗函数的形状

程序集图形文件如下：

（1）矩形窗

>

N=10;

w=boxcar（N）;

nn=[0:

N-1];

plot（nn,w）;

（2）汉宁窗

w=hanning（N）;

（3）汉明窗

w=hamming（N）;

（4）巴特利窗

w=bartlett（N）;

（5）布莱克曼窗

w=blackman（N）;

（6）traing窗

w=triang（N）;

（7）kaiser窗

w=kaiser（N）;

（8）切比雪夫窗

w=chebwin（N）;

2.用MATLAB编程绘制各种窗函数的幅频相应，程序集图形文件如下：

[X,W]=dtft（w,50）;

subplot（211）

plot（W/2/pi,abs（X））

xlabel（'

w'

）,ylabel（'

abs'

）

subplot（212）

plot（W/2/pi,180/pi*angle（X））

phase'

（5）生布莱克曼窗

xlabel（'

（8）切比雪夫窗

3.绘制矩形窗信号的幅频响应，窗长度分别为，N=10，N=20，N=50，N=100

N=10时

subplot（111）

N=20时

N=20;

[X,W]=dtft（w,100）;

N=50时

N=50;

N=100时

N=100;

4.已知周期信号x（t）=0.75+3.4cos2πft+2.7cos4πft+1.5sin3.5πft+2.57πft,其中f=25/16Hz,若截断时间长度分别为信号周期的0.9和1.1倍，试绘制和比较采用下面窗函数提取的x（t）的频谱。

fs=10;

%取抽样频率为10Hz

T0=2.56;

%函数x（t）的周期为T0=2.56

f=25/16;

Ts=1/fs;

N1=0.9*T0/Ts;

%截断时间为信号周期的0.9倍

N1=fix（N1）;

%对N1朝零方向取整

N1-1];

w=boxcar（N1）;

x=0.75+3.4*cos（2*pi*f*nn/fs）+2.7*cos（4*pi*f*nn/fs）+1.5*sin（3.5*pi*f*nn/fs）+2.5*sin（7*pi*f*nn/fs）;

x=w.*x'

;

%用窗函数提取x（t）

[X,W]=dtft（x,500）;

subplot（211）,plot（W/2/pi,abs（X））

title（'

t=0.9*T0'

）,xlabel（'

N2=1.1*T0/Ts;

%截断时间为信号周期的1.1倍

N2=fix（N2）;

N2-1];

w=boxcar（N2）;

x=0.75+3.4*cos（2*pi*f*nn/fs）+2.7*cos（4*pi*f*nn/fs）

+1.5*sin（3.5*pi*f*nn/fs）+2.5*sin（7*pi*f*nn/fs）;

subplot（212）,plot（W/2/pi,abs（X））

t=1.1*T0'

w=hanning（N1）;

%用窗函数提取x（t）

w=hanning（N2）;

w=hamming（N1）;

w=hamming（N2）;

%函数x（t）的周期为T0=2.56

%截断时间为信号周期的0.9倍

%对N1朝零方向取整

w=bartlett（N1）;

%截断时间为信号周期的1.1倍

w=bartlett（N2）;

w=blackman（N1）;

w=blackman（N2）;

w=triang（N1）;

w=triang（N2）;

w=kaiser（N1）;

w=kaiser（N2）;

w=chebwin（N1）;

w=chebwin（N2）;