addzeros.m

来自「用matlab编写的数字信号处理程序」· M 代码 · 共 36 行

clc,clear all,close all;
n=[0:1:999];
xn=10*cos(0.2*pi*n)+1.2*cos(0.25*pi*n)+4.6*cos(0.6328*pi*n);    %先定义一个长序列，处理时再截取

%先采集128点，且加汉明窗，做DFT%%%%%%%%%%%%%%%%%%%%%%%%%
w1=hamming(128);
xn1=xn(1:128);
xn11=xn1.*w1';                 %xn11就是加了汉明窗的128点数据

%补齐0点得到256点数据，在做DFT%%%%%%%%%%%%%%%
xn2=zeros(1,128);
xn22=[xn11,xn2];
X2=fft(xn22,256);
subplot(2,1,1);
plot([0:255],X2);
xlabel('n 256-point');
ylabel('dft result --w');

%补齐0点得到512点数据%%%%%%%%%%%%%%%%%%%%%
xn3=zeros(1,384);
xn33=[xn11,xn3];
X3=fft(xn33,512);
subplot(2,1,2);
plot([0:511],X3);
xlabel('n 512-point');
ylabel('dft result --w');