学习目标:
掌握快速傅里叶变换,数字滤波器等理论知识。
学习内容:
傅里叶变换
数字滤波器
学习时间:
2020年8月28日
学习产出:
傅里叶变换阅读知乎文章:
傅里叶分析之掐死教程
转载CSDN的一篇文章:
数字滤波器的理论知识
例子:
对一个数字滤波器进行傅里叶变换,在matlab R2019b中实现
function example
%UNTITLED1 Summary of this function goes here
% Detailed explanation goes here
% 生成原始信号
fs = 3000;
x = 0:1/fs:20;%采样频率必须至少是信号中最大频率分量频率的两倍
data = 0.5*sin(2*pi*300*x);
function my_plot(data,fs)
% 对数据data进行傅里叶变换
N = length(data);
xdft = fft(data);
xdft = xdft(1:floor(N/2)+1);
freq = 0:fs/N:fs/2;
plot(freq, (2/N)*abs(xdft))
end
% 对原始信号进行傅里叶变换并观察频域情况
subplot(211)%绘制子图,意思是大图含有2行1列共2个子图,正在绘制是第一个
my_plot(data ,fs)
% 对原始信号进行滤波
N = 4; % Order
Fc1 = 5; % First Cutoff Frequency
Fc2 = 60; % Second Cutoff Frequency
% Construct an FDESIGN object and call its BUTTER method.
h = fdesign.bandpass('N,Fc1,Fc2', N, Fc1, Fc2, fs);
Hd = butter(h);
filter_data = filter(Hd, data);
% 对观察后的信号进行傅里叶变化并观察频域情况
subplot(212)
my_plot(filter_data ,fs)
end
运行之后效果图: