【MatLab fdatool】Butterworth 低通 & 带阻 滤波器 （MatLab &C++）

最新推荐文章于 2024-05-31 17:28:57 发布
最新推荐文章于 2024-05-31 17:28:57 发布
阅读量1.8k 收藏 9
点赞数
分类专栏： 滤波 文章标签： 数字信号处理
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_44410512/article/details/110148709
版权
MatLab:

1，终端中输入 fdatool
在这里插入图片描述
2，选择参数，如下
在这里插入图片描述

3，【Edit】 → \rightarrow →【convert structure…】
在这里插入图片描述

4，【Edit】 → \rightarrow →【convert to Single section】
5，【File】 → \rightarrow →【Generate MATLAB code 】 → \rightarrow →【Filter Design Function】 保存 我这里选择的是 Lowpass.m
6，【File】 → \rightarrow →【Export…】
在这里插入图片描述

Butterworth 时域滤波公式
y [ i ] = ∑ j = 0 N + 1 ( b [ j ] × x [ i − j ] ) − ∑ j = 1 N + 1 ( a [ j ] × y [ i − j ] ) y[i]=\sum_{j=0}^{N+1}(b[j]\times x[i-j])-\sum_{j=1}^{N+1}(a[j]\times y[i-j]) y[i]=j=0∑N+1​(b[j]×x[i−j])−j=1∑N+1​(a[j]×y[i−j])
其中 b为Num、a为Den、N为阶数

将Num Den 导出
编写c++语言代码：
产生频率为 50、300的正弦波
改变printf 生成 x和y

code :
#include <bits/stdc++.h>
using namespace std;
#define pi 3.1415926
#define S 1000
typedef struct {
	double InPut;
	double OutPut;
	double x[11];
	double y[11];
	double b[11]={0.00049945,0.0049945,0.022475,0.059934,0.10489,0.12586,0.10489,0.059934,0.022475,0.0049945,0.00049945};
	double a[11]={1,-1.9924,3.0195,-2.8185,2.0387,-1.0545,0.41445,-0.11572,0.022499,-0.0026689,0.00014876};
	int Begin=0;
	int End=0;
	int M=0;
	int N=11;
}IIR;

void IIRCalc(IIR &Lowpass){
	Lowpass.x[Lowpass.Begin]=Lowpass.InPut;
	Lowpass.y[Lowpass.Begin]=Lowpass.OutPut=0;
	Lowpass.End=(Lowpass.Begin)%Lowpass.N;
	for(int i=0;i<Lowpass.N;i++){
		if(Lowpass.End<0){
			Lowpass.End=(Lowpass.End+Lowpass.N)%Lowpass.N;
		}
		Lowpass.OutPut+=Lowpass.b[i]*Lowpass.x[Lowpass.End];
		Lowpass.OutPut-=Lowpass.a[i]*Lowpass.y[Lowpass.End];
		--Lowpass.End;
	}
	Lowpass.y[Lowpass.Begin]=Lowpass.OutPut;
	++Lowpass.Begin;
	if(Lowpass.Begin>=Lowpass.N){
		Lowpass.Begin=0;
	}
}
IIR Lowpass;
int main(){
	freopen("out.txt", "w", stdout);
	for(int i=0;i<S;i++){
		if(sin(2*pi*100*i/1000)>0) Lowpass.InPut=5;
		else Lowpass.InPut=0;
		IIRCalc(Lowpass);
		printf("%lf ,", Lowpass.OutPut);
	}
}

将 x y的数据复制到 MATLAB中，查看结果

MATLAB code:
x=[0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,5.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,];
y=[0.000000 ,0.002497 ,0.032445 ,0.196948 ,0.740984 ,1.931973 ,3.694213 ,5.337129 ,5.885702 ,4.878327 ,2.802004 ,0.703358 ,-0.537560 ,-0.582657 ,0.443640 ,2.164129 ,4.042906 ,5.420045 ,5.702908 ,4.707834 ,2.829188 ,0.841061 ,-0.478677 ,-0.649599 ,0.363548 ,2.169119 ,4.104456 ,5.449851 ,5.673689 ,4.670028 ,2.830362 ,0.869649 ,-0.464139 ,-0.662928 ,0.345608 ,2.169334 ,4.117860 ,5.456904 ,5.667553 ,4.661501 ,2.830322 ,0.875946 ,-0.460722 ,-0.665756 ,0.341554 ,2.169250 ,4.120819 ,5.458558 ,5.666250 ,4.659574 ,2.830251 ,0.877336 ,-0.459921 ,-0.666356 ,0.340639 ,2.169201 ,4.121472 ,5.458945 ,5.665974 ,4.659139 ,2.830222 ,0.877642 ,-0.459735 ,-0.666483 ,0.340432 ,2.169184 ,4.121616 ,5.459035 ,5.665916 ,4.659041 ,2.830212 ,0.877709 ,-0.459691 ,-0.666509 ,0.340386 ,2.169179 ,4.121647 ,5.459056 ,5.665904 ,4.659019 ,2.830209 ,0.877724 ,-0.459681 ,-0.666515 ,0.340375 ,2.169177 ,4.121654 ,5.459061 ,5.665901 ,4.659014 ,2.830208 ,0.877727 ,-0.459679 ,-0.666516 ,0.340373 ,2.169177 ,4.121656 ,5.459062 ,5.665901 ,4.659013 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,2.830208 ,0.877728 ,-0.459678 ,-0.666516 ,0.340372 ,2.169177 ,4.121656 ,5.459063 ,5.665901 ,4.659012 ,];
subplot(6,1,1);
plot(x);

subplot(6,1,2);
y1=abs(fft(x)).*2/1000;
plot(y1);

y3 = filter(loss,x);
subplot(6,1,3);
plot(y3);

subplot(6,1,4);
y4=abs(fft(y3)).*2/1000;
plot(y4);

subplot(6,1,5);
plot(y);

subplot(6,1,6);
y5=abs(fft(y)).*2/1000;
plot(y5);

在这里插入图片描述

同样 带阻滤波器如下：

在这里插入图片描述
在这里插入图片描述

code

#include <bits/stdc++.h>
using namespace std;
#define pi 3.1415926
#define S 50000
double x[200000];
double y[200000];
double b[20]={ 0.36979, -1.1967, 3.3981,-5.7896,8.6699,-9.2278,8.6699,-5.7896,3.3981,-1.1967,0.36979};
double a[20]={1,-2.6071, 5.8006,-7.9918,9.6486,-8.442,6.5155,-3.631,1.7741,-0.52858,0.13674};

int main(){
  freopen("out.txt", "w", stdout);
	for(int i=15;i<S;i++){
		x[i]=sin(2*pi*5000*i/S)+sin(2*pi*6000*i/S)+sin(2*pi*7000*i/S)+sin(2*pi*8000*i/S)+sin(2*pi*9000*i/S)+sin(2*pi*10000*i/S)+sin(2*pi*15000*i/S);
		for(int j=0;j<11;j++)
			y[i]+=b[j]*x[i-j];
		for(int j=1;j<11;j++)
			y[i]-=a[j]*y[i-j];
		printf("%lf ,", y[i]);
	}
}

MATLAB 显示 同上

_-Y-_-Y-_
关注
  • 0
    点赞
  • 踩
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
带阻滤波MATLAB代码
11-20
带阻滤波器MATLAB的S函数代码，测试过
带阻滤波器matlab代码-IIR_Butterworth_Filter_C_Sharp:适用于C＃的IIR巴特沃思滤波器
05-27
带阻滤波器matlab代码C＃代码用于计算Butterworth滤波器的系数并过滤数据 该代码计算带通，带阻，低通和高通巴特沃斯滤波器的系数。 它还过滤数据，但不应用零相位延迟。 每个过滤器函数将返回2行x N个系数的2D向量，其中行1 =分子，行2 =分子。 方法“ Check_stability_iir”可用于检查过滤器的稳定性。 Please, keep in mind that if the filter is unstable, numerical instability leading to numerical overflow might happen when the order selected is extremely high. 如果发生这种情况，程序可能会在分母处分配默认值10 ^ 10。 带通：函数为“ double [] [] Lp2bp（double W_f1，double W_f2，int order_filt）”。 前两个自变量是归一化的两个截止频率（f1 / SF，f2 / SF），这意味着截止频率必须在（0,1）的区间内，最后一个自变量是阶数。 带
matlab 设计巴特沃茨滤波器，然后在c++中使用
random_liu_2013
07-23 607
一、用matlab设计巴特沃茨滤波器 %1 低通 fd = 10;%hz 截止频率 fs = 167;%hz 采样频率 wn = fd/(0.5*fs); %butter 第一个参数为阶数，下面使用的是一阶 [a b]= butter(1,wn,'low') %2 高通 fd = 10;%hz截止频率 fs = 167;%hz采样频率 wn = fd/(0.5*fs); [c d]= butter(1,wn,'high') %小测试 t = 0:1/fs:1; x = cos(2*...
【matlab】fdatool设计IIR和FIR滤波器
最新发布
weixin_43853598的博客
05-31 443
【matlab】fdatool设计IIR和FIR滤波器
巴特沃斯带阻数字滤波器设计matlab程序及仿真图_滤波器_IIR_Butterworth_LPF_数字滤波器_iir_slip
09-30
1.设计IIR巴特沃斯低通滤波器2.将滤波器用于加噪声信号的处理3.改变参数指标比较不同
用Matlab设计Butterworth滤波器后在c/c++中使用
易
10-08 985
设置好参数后，点最下方的“Design Filter”，生成参数。 然后在菜单File下Export导出参数： 导出对话框 点Export后即可在工作区看到导出的数据Hd. 最后使用命令行将参数转换为C语言可用的参数表： [b,a]=sos2tf(Hd.sosMatrix,Hd.ScaleValues) 可得到b和a参数表： b = 1.0e-03 * 0.0231 0...
基于谐波检测中的数字低通滤波器的MATLAB设计
01-19
由于电力系统中非线性电子元件的大量使用，使得谐波污染问题日益严重，当谐波含量超过一定限度时就可能对电网和用户造成极大的危害，且增加线路损耗，降低线路传输...本文利用MATLAB中的FDAtool实现了对数字低通滤波
用Matlab的FDAtool生成IIR滤波器参数以及参数生成C 语言文件
04-28
"Matlab的FDAtool生成IIR滤波器参数以及参数生成C语言文件" 本文主要讲述了使用Matlab的FDAtool生成IIR滤波器参数以及参数生成C语言文件的方法。在数字滤波器设计中，IIR滤波器是一种常用的滤波器类型，它的设计...
用Matlab的FDAtool生成IIR滤波器参数以及参数生成C 语言文件.pdf
07-04
在这个过程中，模拟低通原型滤波器的设计指标被用来创建数字滤波器的传输函数。 在示例中，设计了一个用于去除50Hz电频干扰并保留语音信号（80-3200Hz）的带通滤波器，采样率为8kHz。通过调整FDAtool的参数，确定了...
基于MATLAB的FIR带阻数字滤波器的设计.pdf
10-31
在本文中，我们使用MATLAB的FDATool设计了一个FIR带阻数字滤波器，用于消除50Hz的噪声。实验结果表明，该滤波器可以很好地消除50Hz的噪声。 数字滤波器的设计原理是对提出的设计要求，通过计算使物理可实现的数字...
将Matlab的fdatool生成的Coefficients导出
08-03
Matlab的fdatool是一款强大的滤波器设计和分析工具，它允许用户根据特定需求创建自定义滤波器。本文将详细介绍如何利用fdatool生成的系数（Coefficients）并将其导出到Altera的FIR（Finite Impulse Response）...
用于设计高阶巴特沃斯 IIR 和均衡滤波器的 C++ 类_C++_代码_相关文件_下载
07-12
Butterworth 滤波器设计是 C++ 类的集合和配套的单元测试套件，用于使用双线性变换设计高阶 Butterworth IIR 和 EQ 滤波器。 生成的滤波器系数被拆分为级联双二阶部分，以便在任何双二阶或二阶部分 (SOS) 实现中轻松使用。 特征 低通、高通、带通和带阻 IIR 和 EQ 滤波器设计 低和高搁架滤波器设计 高阶参数提升/削减 EQ 滤波器设计 *Biquad 和 Biquad Chain 实现（用于过滤具有级联 biquad 部分的音频缓冲区） 紧凑、可读且注释良好的代码库 更多详情、使用方法，请下载后阅读README.md文件
matlab代码仿真设计的导通低通带通带阻四种fir数字滤波器
09-21
matlab代码仿真设计的导通低通带通带阻四种fir数字滤波器，亲测可用，并且做了非常详细的代码说明，
带阻滤波器课程设计
10-28
详细描述了带阻滤波原理，给出了不同种类的带阻滤波程序，同时对程序运行结果进行了仔细的分析。
c实现巴特沃斯低通滤波器
10-25
用c实现巴特沃斯低通滤波器，代码通俗易懂，易于上手。
FIRFilter(C# FIR滤波器——含低通、高通、带通与带阻).rar
12-19
最近需要用到Fir滤波器，在网上也看了不少资料，发现一个稍微能用的（https://blog.csdn.net/BIGFatming/article/details/92386914），主要代码也是直接copy的，但是在使用过程中发现，三角窗的实现好像不对，而且只实现了低通的，根据该文章内容，我自己重写了一个，包括三角窗函数的实现，添加了高通、带通、带阻的实现。环境Visual Studio2015，c#。使用方法：https://blog.csdn.net/hopeless123/article/details/103616011
C++ Butterworth N阶滤波器设计
huangjunsheng123的博客
03-07 2414
多通道信号实时数字滤波
FIR滤波器，低通、高通、带通、带阻VC实现
热门推荐
hitwhzhongqiu的专栏
12-15 2万+
1.前言：数字信号处理相关知识准备 通常来说，一种理想滤波器的频率响应是很容易理解的，如图所示。   图1 滤波器频响 以低通为例，滤波器频率响应函数为 。 所谓滤波器处理的过程，简单来说，可以用公式 来表示，由卷积的性质可以知道，该公式的另一种形式为 其中x(n)为要处理的数据序列，h(n)为逼近滤波器的时域响应 其中，hd(n)为对应不同类型滤波器的单位冲击响应，
matlab高通滤波函数,Butterworth高通滤波器matlab程序(原创)
weixin_32204175的博客
03-16 2628
用matlab编写的Butterworth高通滤波程序，可以运行，结果如截图源代码：clear all;close all;P=input('input image data:');subplot(221);imshow(P);title('原始图像');d=size(P);if(d(3)>1)P=rgb2gray(P);endsubplot(222);imshow(P);title('灰度...
在matlab中不适用fdatool工具设计FIR带阻滤波器
05-31
好的，以下是在MATLAB中使用fir1函数自己手动设计FIR带阻滤波器的步骤： 1. 确定需要设计的滤波器的特性，包括截止频率、通带和阻带范围等参数。 2. 使用fir1函数计算出滤波器的系数。fir1函数的第一个参数是滤波器的阶数，第二个参数是截止频率，第三个参数是滤波器类型（此处为带阻滤波器），第四个参数是滤波器系数。 3. 将需要进行滤波的信号输入fir1函数中，并将输出结果保存到另一个变量中。 以下是一个示例代码： ```matlab % 设计一个FIR带阻滤波器 fs = 1000; % 采样率 f1 = 50; % 带阻滤波器的下限截止频率 f2 = 100; % 带阻滤波器的上限截止频率 Wp = [f1 f2]/(fs/2); % 截止频率归一化 Ws = [f1-10 f2+10]/(fs/2); % 通带和阻带范围 Rp = 1; % 通带最大允许波动 Rs = 40; % 阻带最大允许波动 [n, Wn] = buttord(Wp, Ws, Rp, Rs); % 计算最优的滤波器阶数和截止频率 b = fir1(n, Wn, 'stop'); % 计算滤波器系数 % 调用fir1函数进行滤波 x = randn(1, 1000); % 随机生成一个信号 y = filter(b, 1, x); % 进行滤波 ``` 希望这可以帮助你手动设计一个FIR带阻滤波器。如果你有任何问题，请随时问我。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值