| 查看: 909 | 回复: 3 | ||
| 【奖励】 本帖被评价1次，作者baobiao007增加金币 0.8 个 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[资源]
【分享】时域褶积带通滤波c程序【原创，无重复】
|
||
|
#include #include #include #include"FFT.h" /*带通滤波因子 num--滤波因子长度（奇数） f0--通频带中心频率 hf--通频带半宽度 dt--时间采样间隔(与输入信号相同) */ void FilterH(float h[],int num,float f0,float hf,float dt) { int i,j; double pi=3.1415926; float t; float *w; //产生汉明时窗w[] w=(float *)calloc(num,sizeof(float)); for(i=0; i<=num/2; i++) w[ i ]=0.54+0.46*cos(2*pi*i/num); for(i=num/2+1; i //产生滤波因子 h[0]=4.0*hf; for(i=1; i h[ i ]=0.0; if(i<=num/2) { t=i*dt;j=i;} else { j=num-i; t= -j*dt;} h[ i ]=2.0*sin(2*pi*hf*t)*cos(2*pi*f0*t)*w[ j ]/(pi*t); } } /*输入信号与滤波因子褶积 h[]--滤波因子,x[]--输入信号 保证x[]更长,结果存放在y[],长度也为m */ void ConvFilter(float h[],int n,float x[],int m,float y[],float dt) { int i,j,k; int half; half=(n-1)/2; for(i=0; i y[ i ]=0.0; for(j= -half; j<=half; j++) { k=i-j; if(j<0 && k if(j>=0 && k>=0) y[ i ]+=h[j]*x[k]*dt; } } } //产生合成信号 //30 60 90hz正弦信号相加 void Generate(float x[],int num) { int i; double pi=3.1415926; float dt=0.002; float f1,f2,f3; f1=30.0; f2=60.0; f3=90.0; for(i=0; i } void main() { const int N=101;//滤波因子长度 const int M=128;//信号长度 float Inisig[M]={0};//要滤波的信号 float Afsig2[M]={0};//滤波后的信号 float Afsig[M]={0};//理想情况下滤波后的信号 float h[N];//滤波因子 int i; FILE *fp1,*fp2,*fp3; fp1=fopen("inisig.txt","w" );//合成信号 fp2=fopen("afsig.txt","w" );//理想50hz信号 fp3=fopen("afsig2.txt","w" );//滤波得到的50hz信号 //产生信号,20hz,50hz,80hz正弦信号叠加 double pi=3.1415926; float dt=0.002; float f1,f2,f3; f1=20.0; f2=50.0; f3=80.0; for(i=0; i Inisig[ i ]=Afsig+sin(2*pi*f1*i*dt)+sin(2*pi*i*f3*dt); } FilterH(h,N,50.0,10.0,dt); //滤去低于70hz和高于80hz信号 ConvFilter(h,N,Inisig,M,Afsig2,dt); for(i=0; i fprintf(fp1,"%f\n",Inisig[ i ]); fprintf(fp2,"%f\n",Afsig[ i ]); fprintf(fp3,"%f\n",Afsig2[ i ]); } fclose(fp1); fclose(fp2); fclose(fp3); } 滤波前：  滤波后：  [ Last edited by baobiao007 on 2011-2-14 at 20:42 ] |
回复此楼» 猜你喜欢
首次在隐伏矿体上方发现金、硫酸铅、硝酸铅、三氧化钨等含金属纳米微粒
已经有8人回复
-
求助：有没有大神有流体包裹体计算软件
已经有1人回复
-
地球物理学和空间物理学论文润色/翻译怎么收费?
已经有62人回复
-
282求调剂
已经有0人回复
-
香港城市大学 数据科学系+能源与环境学院 李堂华 课题组诚招全奖博士/博士后/研究助理
已经有51人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Project Euler 50 欧拉工程 50 题
已经有12人回复
- 外国导师国内兼职的可以通过csc么
已经有7人回复
- 关于2012年度国家自然科学基金政策重大调整的通知
已经有1272人回复
- 夫妻两个都是J1，一起面签会不容易通过吗？
已经有49人回复
- 《科学通报》投稿 word排版求助
已经有14人回复
- 联培顺利通过CSC录取。但是美国导师很可能要跳到别的大学了，咋办？
已经有51人回复
- 大家一般开始一个新课题，会通过哪些渠道去找综述呢？
已经有22人回复
- 【求助】Weickert的各向异性扩散方程的滤波方法的matlab程序代码
已经有8人回复
4楼2017-05-06 21:14:26
2楼2011-04-14 10:28:20
3楼2011-04-14 12:32:15
20