一阶RC高通滤波器详解(仿真+matlab+C语言实现)
如果本文帮到了你,帮忙点个赞;
如果本文帮到了你,帮忙点个赞;
如果本文帮到了你,帮忙点个赞;
HPF 一阶RC高通滤波器详解(仿真+matlab+C语言实现)
LPF 一阶RC低通滤波器详解(仿真+matlab+C语言实现)
预备知识
高通滤波器(HPF-high pass filter
)可以滤除频率低于截止频率的信号,类似的还有低通滤波器,带通滤波器,带阻滤波器。一阶RC高通滤波器的电路如下图所示;
关于电容
首先对电容的几个公式做一下补充;
电容大小CCC满足;
C=qV⋯①
C=\cfrac{q}{V} \cdots①
C=Vq⋯①
其中 qqq 是电容所带的电荷量,VVV 是电容两端的电势差;
另外,电流相当于单位时间流过导体的电荷量;因此电流 iii 满足;
i=dqdt⋯② i = \cfrac{dq}{dt}\cdots②i=dtdq⋯②
根据①,②可以得到电容大小 CCC 和电容的电流 iii 以及两端电压 VVV 的关系;
i(t)=Cdv(t)dti(t) = C\cfrac{dv(t)}{dt}i(t)=Cdtdv(t)
HPF的推导
由以上电路可知,假设电流为 i(t)i(t)i(t),则可知
{Vout=Ri(t)⋯③i(t)=CdQc(t)dt⋯④
\begin{cases}
V_{out} = Ri(t) \cdots③ \\
i(t) = C\cfrac{dQ_c(t)}{dt} \cdots④\\
\end{cases}
⎩⎨⎧Vout=Ri(t)⋯③i(t)=CdtdQc(t)⋯④
电容两端的电压为 Vc(t)V_c(t)Vc(t) 根据基尔霍夫定律,满足;
−Vin+VC+Vout=0-V_{in} + V_C + V_{out} = 0−Vin+VC+Vout=0
所以结合①,③,④可以得到;
Qc(t)=C(Vin(t)−Vout(t))⋯⑤Q_c(t) = C( V_{in}(t) - V_{out}(t)) \cdots ⑤Qc(t)=C(Vin(t)−Vout(t))⋯⑤
根据 ③,④,⑤ 可以得到以下关系;
Vout=C(dVindt−dVoutdt)⏞I(t)R=RC(dVindt−dVoutdt)⋯⑥
V_{out} = \overbrace{ C( \cfrac{dV_{in}}{dt} - \cfrac{dV_{out}}{dt}) }^{I(t)} R = RC(\cfrac{dV_{in}}{dt} - \cfrac{dV_{out}}{dt}) \cdots ⑥
Vout=C(dtdVin−dtdVout)I(t)R=RC(dtdVin−dtdVout)⋯⑥
将方程进行离散化,如果输入VinV_{in}Vin和输出输入VoutV_{out}Vout按照 △T\bigtriangleup_{T}△T的时间采样,那么可以将输入和输出序列化,则
VinV_{in}Vin序列化为:
(x1,x2,x3⋯ ,xn−1,xn)(x_{1},x_{2},x_{3}\cdots,x_{n-1},x_{n})(x1,x2,x3⋯,xn−1,xn)
VoutV_{out}Vout序列化为:
(y1,y2,y3⋯ ,yn−1,yn)(y_{1},y_{2},y_{3}\cdots,y_{n-1},y_{n})(y1,y2,y3⋯,yn−1,yn)
根据⑥式可以进行离散化,因此最终滤波输出的序列 yiy_{i}yi 如下所示;
yi=RC(xi−xi−1△T−yi−yi−1△T)⋯⑦
y_{i} = RC(\cfrac{x_{i}-x_{i-1}}{\bigtriangleup_{T}} -\cfrac{y_{i}-y_{i-1}}{\bigtriangleup_{T}} )\cdots⑦
yi=RC(△Txi−xi−1−△Tyi−yi−1)⋯⑦
将⑦再进一步简化得到;
yi=αyi−1+α(xi−xi−1)
y_{i} = \alpha y_{i-1} + \alpha(x_i - x_{i-1})
yi=αyi−1+α(xi−xi−1)
其中 α=RCRC+△T\alpha = \cfrac{RC}{RC+\bigtriangleup_{T}}α=RC+△TRC
所以换成得到;
RC=△T(α1−α)⋯⑧RC = \bigtriangleup_{T}(\cfrac{\alpha}{1-\alpha}) \cdots⑧RC=△T(1−αα)⋯⑧
另外截止频率和低通滤波器的相同;
fc=12πRCf_c = \cfrac{1}{2\pi RC}fc=2πRC1
将⑧式代入可以得到截止频率和 α\alphaα 的关系;
fc=1−α2πα△Tf_c = \cfrac{1-\alpha}{2\pi \alpha \bigtriangleup_{T}}fc=2πα△T1−α
这个公式便于简化后面程序以及截止频率的计算。
simulink 仿真
这里根据公式⑥构建simulink
的子模块subsystem
;
Vout=RC(dVindt−dVoutdt)
V_{out} = RC(\cfrac{dV_{in}}{dt} - \cfrac{dV_{out}}{dt})
Vout=RC(dtdVin−dtdVout)
具体如下所示;
整体的仿真如下图所示;
其中Sine Wave
频率设置为2*pi*40
,频率为40
赫兹;
其中Sine Wave1
频率设置为2*pi*4
,频率为4
赫兹;
所以这里需要使得2*pi*4
的信号衰减,所以根据,截止频率 fcf_cfc 的计算公式,可以改变增益的值,具体如下所示;
这里RC
增益为0.005
,因此
fc=12πRC=12π∗0.005≈31.8f_c = \cfrac{1}{2\pi RC} = \cfrac{1}{2\pi *0.005} \approx 31.8fc=2πRC1=2π∗0.0051≈31.8
simulink 运行结果
matlab 实现
matlab
根据以下这个公式进行数字滤波器的设计;
yi=αyi−1+α(xi−xi−1)
y_{i} = \alpha y_{i-1} + \alpha(x_i - x_{i-1})
yi=αyi−1+α(xi−xi−1)
另外 α\alphaα 的值如何确定需要参考⑧式;
Serial = 0:0.1:100;
Fs = 1;
Phase = 0;
Amp = 1;
N0 = 2*pi*Fs*Serial - Phase;
X0 = Amp*sin(N0);
subplot(4,1,1);
plot(X0);
Fs = 0.02;
N1 = 2*pi*Fs*Serial - Phase;
X1 = Amp*sin(N1);
subplot(4,1,2);
plot(X1);
X2=X0+X1;
subplot(4,1,3);
plot(X2);
len = length(X2);
X3=X2;
p=0.75;
for i=2:len
X3(i) = p*X3(i-1)+p*(X2(i)- X2(i-1))
end
subplot(4,1,4);
plot(X3);
简单地分析一下,代码中的X1
,X2
,X3
;
- X1频率为
1
- X2频率为
0.02
{α=p=0.75△T=0.1fc=1−α2πα△T \begin{cases}
\alpha = p=0.75 \\
\bigtriangleup_{T} = 0.1 \\
f_c = \cfrac{1-\alpha}{2\pi \alpha \bigtriangleup_{T}} \end{cases}⎩⎪⎪⎪⎨⎪⎪⎪⎧α=p=0.75△T=0.1fc=2πα△T1−α
因此可以得到截止频率如下;
fc=0.252π∗0.75∗0.1≈0.53f_c=\cfrac{0.25}{2\pi *0.75* 0.1} \approx 0.53fc=2π∗0.75∗0.10.25≈0.53
matlab 运行结果
C语言实现
typedef struct
{
int16_t Input[2];
int16_t Output[2];
int32_t FilterTf;
int32_t FilterTs;
int32_t Ky;
} high_filter;
void high_filter_init(high_filter *v);
int16_t high_filter_calc(high_filter *v);
其中;
FilterTs
为采样时间 △T\bigtriangleup_{T}△T;FilterTf
为RC
时间常数;Input[0]
表示 xix_ixi;Input[1]
表示 xi−1x_{i-1}xi−1;Output[0]
表示 yiy_iyi;Output[1]
表示 yi−1y_{i-1}yi−1;Ky
表示RCRC+△T\cfrac{RC}{RC+\bigtriangleup_{T}}RC+△TRC;
参考公式如下所示;
void high_filter_init(high_filter *v){
v->Ky = v->FilterTf*1024/(v->FilterTs + v->FilterTf);
}
int16_t high_filter_calc(high_filter *v){
int32_t tmp = 0;
tmp = ((int32_t)v->Ky*v->Output[1] + v->Ky*(v->Input[0] - v->Input[1]))/1024;
if(tmp>32767){
tmp = 32767;
}
if( tmp < -32768){
tmp = -32768;
}
v->Output[0] = (int16_t)tmp;
v->Output[1] = v->Output[0];
v->Input[1] = v->Input[0];
return v->Output[0];
}
一阶RC高通滤波器详解(仿真+matlab+C语言实现)的更多相关文章
- 一阶RC低通滤波器详解(仿真+matlab+C语言实现)
文章目录 1 预备知识 2 simulink 仿真 3 simulink 运行结果 4 matlab实现 5 matlab运行结果 6 C语言实现 7 C语言运行结果 如果本文帮到了你,帮忙点个赞: ...
- 《算法详解:C++11语言描述》已出版
经过漫长的编写.修订和印刷过程,书籍<算法详解:C++11语言描述>终于出版了!目前本书已在各大电商平台上架,搜索书名即可找到对应商品.本书的特色在于: 介绍最新的C++11.C++14和 ...
- 机器学习笔记----Fuzzy c-means(FCM)模糊聚类详解及matlab实现
前言:这几天一直都在研究模糊聚类.感觉网上的文档都没有一个详细而具体的讲解,正好今天有时间,就来聊一聊模糊聚类. 一:模糊数学 我们大家都知道计算机其实只认识两个数字0,1.我们平时写程序其实也是这样 ...
- 遗传学详解及Matlab算法实现
遗传学算法概述 从之前转载的博客<非常好的理解遗传算法的例子>中可以知道,遗传学算法主要有6个步骤: 1. 个体编码 2. 初始群体 3. 适应度计算 4. 选择运算 5. 交叉运算 6. ...
- 详解keil采用C语言模块化编程时全局变量、结构体的定义、声明以及头文件包含的处理方法
一.关于全局变量的定义.声明.引用: (只要是在.h文件中定义的变量,然后在main.c中包含该.h文件,那么定义的变量就可以在main函数中作为全局变量使用) 方法1: 在某个c文件里定义全局变量后 ...
- 1、WWDC2014 详解OSX/iOS8/Swift语言
OS X 10新特性 1.系统图标扁平化. 2.系统菜单栏可定制. 3.新的通知中心. 4.新的Spotlight. 5.新的iCloud Drive,新增文档同步功能(学Windows). 6.新的 ...
- ESP8266使用详解--基于Lua脚本语言
这些天,,,,今天终于看到了希望,,,天道酬勤 先说实现的功能...让ESP8266连接无线网,然后让它建立服务器,,我的客户端连接上以后,发给客户端发数据模块打印到串口,,往ESP8266串口里发数 ...
- 【实验课选题详解】用C语言实现万年历
题目要求 编程实现万年历,要求: 可根据用户输入或系统日期进行初始化,如果用户无输入则显示系统日期所在月份的月历,并突出显示当前日期: 可根据用户输入的日期查询,并显示查询结果所在月份的月历,突出显示 ...
- 「实验课选题详解」用C语言实现万年历
题目要求 编程实现万年历,要求: 可根据用户输入或系统日期进行初始化,如果用户无输入则显示系统日期所在月份的月历,并突出显示当前日期: 可根据用户输入的日期查询,并显示查询结果所在月份的月历,突出显示 ...
随机推荐
- Coin Change UVA
Suppose there are 5 types of coins: 50-cent, 25-cent, 10-cent, 5-cent, and 1-cent. We want to makech ...
- JavaScript函数作用域和声明提前(3.10.1 page.57)
<h4>3.函数作用域和声明提前</h4> <p> <!--<script type="text/javascript">-- ...
- Python爬取养眼图片
1.准备 各位绅士们,你可能会觉得疫情在家无聊,那么现在我们的Python语言可以满足你们的需求.项目需要的工具(1)Python3(2)requests库requests库可以通过代码pip ins ...
- 虚拟机体验NAS私人云全揭秘:序言——虚拟机体验NAS私人云缘由
"世界在新冠肺炎疫情后将永远改变",对于2020春天在全球蔓延的新冠肺炎疫情,美国前国务卿基辛格做了这样的评价.确实,也改变了我们.春节期间,本着少添乱的原则,响应国家号召,自我隔 ...
- Springboot:员工管理之修改员工(十(8))
构建员工修改请求 com\springboot\controller\EmployeeController.java /*调转到员工修改页 携带员工信息 restful风格*/ @GetMapping ...
- python爬虫面试题集锦及答案
1.爬取数据后使用哪个数据库存储数据的,为什么? - 2.你用过的爬虫框架或者模块有哪些?优缺点? - 3.写爬虫是用多进程好?还是多线程好? - 4.常见的反爬虫和应对方法? - 5.需要登录的网页 ...
- 获取 ProgramData 之类的环境变量(文件夹路径)的值
GetEnvironmentVariable Recognized Environment Variables https://docs.microsoft.com/en-us/windows/dep ...
- Silverlight Tools Beta2更新了中文语言支持
1,似乎是微软偷偷摸摸更新的......刚才无意间发现,已经下载安装并测试,已在中文版的VS2008安装成功.注意下载页面的安装说明: http://www.microsoft.com/downloa ...
- nginx响应超时upstream timed out 问题处理
环境介绍 服务器:centos7.2 应用:tomcat集群 服务:nginx 代理 问题描述: 这段时间,听项目组项目经理和业务需求人员跟我反馈,线上业务人员在操作业务交易时,有时会出现nginx错 ...
- MySQL Change Data Directory
为什么80%的码农都做不了架构师?>>> Stop MySQL using the following command: sudo /etc/init.d/mysql stop ...