| 查看: 957 | 回复: 1 | ||
| 【奖励】 本帖被评价1次，作者baobiao007增加金币 0.8 个 | ||
[资源]
【分享】频率域镶边法带通滤波c程序【原创】
|
||
|
#include #include #include #include"FFT.h" //效果似乎不如时域褶积法 /*构建镶边带通滤波器 f0--通频带中心频率 hf--通频带半宽度 num-滤波器长度 df--频率采样间隔 */ void Filter(float h[],int num,float f0,float hf,float df) { int i; double pi=3.1415926; float bf=hf;//镶边宽度 float f1,f2,f3,f4,f; f2=f0-hf; f3=f0+hf; f1=f2-bf; f4=f3+bf; for(i=0; i<=num/2; i++) { f=i*df; if(f h[ i ]=pow(sin(pi*(f-f1)/(2.0*(f2-f1))),2); else if(f>f2 && f else if(f>=f3 && f<=f4) h[ i ]=pow(sin(pi*(f-f4)/(2.0*(f4-f3))),2); else h[ i ]=0.0; } for(i=num/2+1; i } void main() { const int M=128;//滤波器长度 float dt=0.002;//时间采样间隔 float df=1.0/(M*dt);//频率采样间隔 float Inisigr[M]={0};//要滤波的信号 float Inisigi[M]={0}; float Afsig2r[M]={0};//滤波后的信号 float Afsig2i[M]={0}; float Afsig[M]={0};//理想情况下滤波后的信号 float hr[M],hi[M];//滤波器频谱 int i,k; k=(int)(log(M)/log(2.0)+0.5); 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 f1,f2,f3; f1=20.0; f2=50.0; f3=80.0; for(i=0; i Inisigr[ i ]=Afsig[ i ]+sin(2*pi*f1*i*dt)+sin(2*pi*i*f3*dt); fprintf(fp1,"%f\n",Inisigr[ i ] ); } //形成滤波器 Filter(hr,M,50.0,10.0,df); //对要滤波的信号进行FFT fft(Inisigr,Inisigi,k,1); //频谱相乘,滤去低于40hz高于60hz的信号 for(i=0; i Afsig2r[ i ]=hr[ i ]*Inisigr[ i ]; Afsig2i[ i ]=Inisigi[ i ]*hr[ i ]; } //FFT反变换 fft(Afsig2r,Afsig2i,k,-1); //输出结果 for(i=0; i fprintf(fp2,"%f\n",Afsig[ i ]); fprintf(fp3,"%f\n",Afsig2r[ i ]); } fclose(fp1); fclose(fp2); fclose(fp3); } 滤波前：  滤波后：  [ Last edited by baobiao007 on 2011-2-14 at 22:55 ] |
回复此楼» 猜你喜欢
地球化学论文润色/翻译怎么收费?
已经有272人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问考博报名：普通招考和统招是一个概念吗？
已经有22人回复
- 求助：各位大侠，小弟想问下鱼苗解剖的通过观察内脏知道病情的问题
已经有14人回复
- 怎样通过1HNMR或13CNMR计算物质的含量？
已经有8人回复
- 有通知后，要作的第一件事是什么？
已经有38人回复
- 【化学】【2011-07-18】中国化学会第四届全国多酸化学学术研讨会会议通知（第一轮）
已经有60人回复
- 有自主权的学校自主评审教授、副教授的程序介绍
已经有8人回复
- 【课件】矿井通风与安全（中国矿业大学）
已经有4人回复
- 【讨论】化学专业出身，如何跨行进入高分子领域
已经有16人回复
- 【分享】2010年同济医学院药学院学术夏令营通知（完结）
已经有7人回复
- 收到中国矿业大学面试通知，学校咋样？
已经有44人回复
2楼2017-05-20 13:45:17
