矩形窗、汉宁窗、升余弦窗等Nyquist窗函数matlab仿真代码

@ 矩形窗、汉宁窗、升余弦窗等Nyquist窗函数matlab仿真代码

抑制窄带干扰的一种很自然的思路，就是利用频谱的旁瓣较低的窗函数代替 DFT 中的矩形窗。这种技术只需要对接收机稍加改进，而且不需要对原有的通信协议做任何的变动。因此我们称它为“接收机窗函数技术”。

matlab仿真多种窗函数的时域图和频谱图，容易发现：时域上窗函数越长、边缘越平缓，则频域旁瓣越低。



MATLAB代码如下：

%%%%%%%%%%%%%%%%%%%%%绘制不同Nquist窗对应的频谱%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%对比不同的长度的窗函数，可知窗函数时域越长，边沿越平缓，频谱旁瓣越低
%N长矩形窗，2N长矩形窗，汉宁窗，升余弦窗时域频域对比
%  Autor:lyq 2023-11-14
clc;clear all;close all;
%长度为N的矩形窗
N = 8;
w1 = ones(1,N);%长度为N的矩形窗
figure 
subplot(2,1,1);stem(w1);title('N长矩形窗时域图')
subplot(2,1,2);stem(abs(fft(w1,64)));title('N长矩形窗频域图')%注：DFT的点数会影响频谱周期

%长度为2N的矩形窗
w2 = ones(1,2*N);%长度为2N的矩形窗
figure 
subplot(2,1,1);stem(w2);title('2N长矩形窗时域图')
subplot(2,1,2);stem(abs(fft(w2,64)));title('2N长矩形窗频域图')

%汉宁窗 
w3 = hann(2*N);%返回一个长度为2N点的汉宁窗
% wvtool(w3);%wvtool工具绘制指定窗口的时域图和频域图，频域是傅里叶变换的幅度的平方
figure 
subplot(2,1,1);stem(w3);title('2N长汉宁窗时域图')
subplot(2,1,2);stem(abs(fft(w3,64)));title('2N长汉宁窗频域图')

%升余弦窗函数
Tu = 8;
r = 0.5;%滚降系数
Tv = r*Tu;%扩展宽度
w4 = [(1-sin(pi*((-Tv:-1)-Tv/2)/(r*Tu)))/2,ones(1,Tu-Tv+1),1-(1-sin(pi*((-Tv+1:-1)-Tv/2)/(r*Tu)))/2];
figure 
subplot(2,1,1);stem(w4);title(['滚降系数',num2str(r),'升余弦窗时域图'])
subplot(2,1,2);stem(abs(fft(w4,64)));title(['滚降系数',num2str(r),'升余弦窗频域图'])

简易仿真，欢迎交流指正！