?? exa031101_conv.m
字號:
%----------------------------------------------------------------------------
% exa031101_conv.m, for example 3.11.1
%
% 不使用“fftfilt"文件,自己編程實現一個長序列和一個短序列的卷積。
%----------------------------------------------------------------------------
clear;
% 生成濾波器系數h和混有高斯白噪的正弦信號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分為長度為L的小段
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))
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -