%----------------------------------------------------------------------------
% exa031101_conv.m, for example 3.11.1
% 
% 不使用“fftfilt"文件，自己編程實(shí)現(xiàn)一個(gè)長(zhǎng)序列和一個(gè)短序列的卷積。
%----------------------------------------------------------------------------
clear;

% 生成濾波器系數(shù)h和混有高斯白噪的正弦信號(hào)x
h=fir1(10,0.3,hanning(11));
N=500;p=0.05;f=1/16;
u=randn(1,N)*sqrt(p);%
s=sin(2*pi*f*[0:N-1]);
x=u(1:N)+s;

% 將x分為長(zhǎng)度為L(zhǎng)的小段
L=20;M=length(h);
y=zeros(1,N+M-1);
tempy=zeros(1,M+L-1);
tempX=zeros(1,L);
for k=0:N/L-1
   tempx(1:L)=x(k*L+1:(k+1)*L);
   tempy=conv(tempx,h);
   y=y+[zeros(1,k*L),tempy,zeros(1,N-(k+1)*L)];
end

subplot(211);plot(x)
subplot(212);plot(y(1:N))