搞定傅里叶和数字滤波器理论知识

学习目标:

掌握快速傅里叶变换,数字滤波器等理论知识。

学习内容:

傅里叶变换
数字滤波器

学习时间:

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

运行之后效果图:
上面的代码运行后实现的效果图

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值