examp58.m

来自「ADSP TOOLBOX: Version 2.0 and gui m-file」· M 代码 · 共 14 行

clc,echo on 
%EXAMPLE 58 
xp=[1 2 -1 0 2 3]; hp=2:-1:-3;  % One period of xp and hp
N=length(xp);                   % Length of xp and hp 
YP=fft(hp).*fft(xp);            % Product of FFT of x and h
yp=ifft(YP);yp=real(yp)         % Real result of periodic convolution 
Nr=N+N-1;                       % Length of regular convolution 
z=zeros(1,Nr-N);                % padding zeros 
xz=[xp z];hz=[hp z];            % Zero-padded x and h 
Y=fft(xz).*fft(hz);
y=real(ifft(Y))                 % Display regular conv
echo off %end of example