原创 简单说说卷积的应用

getFilter

function Hd = getFilter
%GETFILTER Returns a discrete-time filter object.

% MATLAB Code
% Generated by MATLAB(R) 9.4 and DSP System Toolbox 9.6.
% Generated on: 11-Aug-2018 21:15:57

Fpass = 2500;  % Passband Frequency
Fstop = 3000;  % Stopband Frequency
Apass = 1;     % Passband Ripple (dB)
Astop = 60;    % Stopband Attenuation (dB)
Fs    = 8192;  % Sampling Frequency

h = fdesign.lowpass('fp,fst,ap,ast', Fpass, Fstop, Apass, Astop, Fs);
15
16Hd = design(h, 'equiripple', ...
17    'MinOrder', 'any', ...
18    'StopbandShape', 'flat');

DrawFFT

function [  ] = DrawFFT( x, Fs )
% DrawFFT 对输入信号进行快速傅里叶变换
% 输入参数：x ：输入信号；   Fs：采样频率
% by Yang Xiaokang , 2016-02

L = length(x);
NFFT = 2^nextpow2(L);                  %确定FFT变换的长度
y = fft(x, NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);    %频率向量
figure;
plot(f, 2*abs(y(1:NFFT/2+1)));        %绘制频域图像
title('幅度谱');
xlabel('Frequency (Hz)');
ylabel('|y(f)|');
end

% 作者：TaylorSe
% 链接：https://www.zhihu.com/question/46038488/answer/144244073
% 来源：知乎
% 著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。

个人博客