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

1、精品文档能实行， c 语言不能显示图像，但是你要懂得如何用c 调用 matlab ，过程比较复杂，这个程序任何参数都可以用，不限制#include<engine.h>#include<math.h>#include<stdio.h>#include<stdlib.h>#include<string.h>double calcN(double fp,double ap, double fs, double as)double i = pow(10,(as/10) - 1;double tmp = ap/10;double j = pow

2、(10,tmp) - 1;double k = log10(i/j);double t = log10(fs/fp);double N = 0.5*k/t + 1;return N;int main(void)int fp,ap,fs,as;printf("请输入 fp:");scanf("%d", &fp);随意编辑精品文档printf("请输入 ap:");scanf("%d", &ap);printf("请输入 fs:");scanf("%d", &am

3、p;fs);printf("请输入 as:");scanf("%d", &as);int f = calcN(fp, ap, fs, as);printf("N=%dn",f);char lin40;char *str = "N=0"lin0='0'char in4;itoa(f,in,10);strcat(lin,str);strcat(lin,in);strcat(lin, "");/Engine *ep;if(!(ep = engOpen(NULL)随意编辑精品文档

4、printf("open failed! please try againn");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(

5、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,

6、 "plot(w,magH);");engEvalString(ep, "axis(0 2 0 1)");engEvalString(ep, "xlabel('w/wc');");engEvalString(ep, "ylabel('|H(jw)|2');");engEvalString(ep, "title('Butterworth analog filter prototype');");随意编辑精品文档engEvalString(ep, "grid;&