巴特沃斯低通滤波器(C语言)(共5页).doc

《巴特沃斯低通滤波器(C语言)(共5页).doc》由会员分享，可在线阅读，更多相关《巴特沃斯低通滤波器(C语言)(共5页).doc（5页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上能实行，c语言不能显示图像，但是你要懂得如何用c调用matlab，过程比较复杂，这个程序任何参数都可以用，不限制#include#include#include#include#includedoublecalcN(doublefp,doubleap,doublefs,doubleas)doublei=pow(10,(as/10)-1;doubletmp=ap/10;doublej=pow(10,tmp)-1;doublek=log10(i/j);doublet=log10(fs/fp);doubleN=0.5*k/t+1;returnN;intmain(void)i

2、ntfp,ap,fs,as;printf(请输入fp:);scanf(%d,&fp);printf(请输入ap:);scanf(%d,&ap);printf(请输入fs:);scanf(%d,&fs);printf(请输入as:);scanf(%d,&as);intf=calcN(fp,ap,fs,as);printf(N=%dn,f);charlin40;char*str=N=0;lin0=0;charin4;itoa(f,in,10);strcat(lin,str);strcat(lin,in);strcat(lin,;);/Engine*ep;if(!(ep=engOpen(NULL)p

3、rintf(openfailed!pleasetryagainn);engSetVisible(ep,0);/mxArray*xx=mxCreateDoubleMatrix(1,N,mxREAL);/mxArray*yy=mxCreateDoubleMatrix(1,N,mxREAL);/memcpy(mxGetPr(xx),N*sizeof(double);/memcpy(mxGetPr(yy),N*sizeof(double);/engPutVariable(ep,xx,xx);/engPutVariable(ep,yy,yy);engEvalString(ep,n=0:0.01:2;);

4、engEvalString(ep,lin);engEvalString(ep,z,p,k=buttap(N););engEvalString(ep,b,a=zp2tf(z,p,k););engEvalString(ep,H,w=freqs(b,a,n););engEvalString(ep,magH=(abs(H).2;);engEvalString(ep,plot(w,magH););engEvalString(ep,axis(0201);engEvalString(ep,xlabel(w/wc););engEvalString(ep,ylabel(|H(jw)|2););engEvalString(ep,title(Butterworthanalogfilterprototype););engEvalString(ep,grid;);getchar();/mxDestroyArray(xx);/mxDestroyArray(yy);/engClose(ep);return0;专心-专注-专业