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

能实行，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)intfp,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)printf(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;);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;);getcha