开发工具:
文件大小: 3kb
下载次数: 0
上传时间: 2008-11-16
详细说明： 用C语言实现巴特沃兹滤波器的原理 #include #include "math.h" void gainc(b,a,n,ns,x,y,len,sign) int n,ns,len,sign; double b[],a[],x[],y[]; { int i,j,k,n1; double ar,ai,br,bi,zr,zi,im,re,den,numr,numi,freq,temp; double hr,hi,tr,ti; n1=n+1; for (k=0;k0 ;i-- ) { re=br; im=bi; br=(re+b[j*n1+i])*zr-im*zi; bi=( re+b[j*n1+i])*zi+im*zr; } ar=0.0; ai=0.0; for (i=n;i>0 ;i-- ) { re=ar; im=ai; ar=(re+a[j*n1+i])*zr-im*zi; ai=(re+a[j*n1+i])*zi+im*zr; } br=br+b[j*n1+0]; ar=ar+1.0; numr=ar*br+ai*bi; numi=ar*bi-ai*br; den=ar*ar+ai*ai; hr=numr/den; hi=numi/den; tr=x[k]*hr-y[k]*hi; ti=x[k]*hi+y[k]*hr; x[k]=tr; y[k]=ti; } switch(sign) { case 1: { temp=sqrt(x[k]*x[k]+y[k]*y[k]); if (temp!=0.0) { y[k]=atan2(y[k],x[k]); } else { y[k]=0.0; } x[k]=temp; break; } case 2: { temp=x[k]*x[k]+y[k]*y[k]; if (temp!=0.0) { y[k]=atan2(y[k],x[k]); } else { temp=1.0e-40; y[k]=0.0; } x[k]=10.0*log10(temp); } } } } ...展开收缩
(系统自动生成,下载前可以参看下载内容)
下载文件列表
相关说明
- 本站资源为会员上传分享交流与学习,如有侵犯您的权益,请联系我们删除.
- 本站是交换下载平台，提供交流渠道，下载内容来自于网络，除下载问题外，其它问题请自行百度。
- 本站已设置防盗链，请勿用迅雷、QQ旋风等多线程下载软件下载资源，下载后用WinRAR最新版进行解压.
- 如果您发现内容无法下载，请稍后再次尝试；或者到消费记录里找到下载记录反馈给我们.
- 下载后发现下载的内容跟说明不相乎，请到消费记录里找到下载记录反馈给我们，经确认后退回积分.
- 如下载前有疑问，可以通过点击"提供者"的名字，查看对方的联系方式，联系对方咨询.
