数字处理实习报告.doc

附1：C语言和F语言程序1. 褶积滤波F语言 PROGRAM EXAM1PARAMETER(PI=3.1415926535) PARAMETER(DT=0.002)DIMENSION H(-50:50),X(100),Z(200) REAL X,H,Z INTEGER IOPEN(1,FILE=INPUT1.DAT,FORM=FORMATTED,STATUS=OLD)READ(1,*)(X(I),I=1,100)DO 10 I=-50,50IF(I.NE.0)THEN H(I)=SIN(2*PI*70*I*DT)/(PI*I*DT)ELSE H(I)=140 ENDIF10 CONTINUE CALL CON(X,H,Z,100,101,200) OPEN(2,FILE=CONLOW.DAT,FORM=FORMATTED)WRITE(2,*) ZDO 20 I=-50,50IF(I.NE.0)THENH(I)=(SIN(2*PI*80*I*DT)- SIN(2*PI*10*I*DT)/(PI*I*0.002)ELSE H(I)=140 ENDIF 20CONTINUE CALL CON(X,H,Z,100,101,200) OPEN(3,FILE=CONBAND.DAT,FORM=FORMATTED)WRITE(3,*) z CLOSE(1)CLOSE(2) CLOSE(3)END C语言/结果是两道，1道（低通结果200点），2道（带通结果200点）#include#include #include #define pi 3.1415926void main() void conv(float x,int m,float h,int n,float y,int l); FILE *fp1,*fp2; float dt; float *x,*h,*y; int i,m,n,l; dt=0.002; m=100; n=101; l=m+n-1; h=(float *)calloc(101,sizeof(float); x=(float *)calloc(100,sizeof(float); y=(float *)calloc(200,sizeof(float); fp1=fopen(INPUT1.DAT,r); fp2=fopen(y1.dat,w); /读入地震数据x(t) for(i=0;im;i+) fscanf(fp1,%f,&xi); fclose(fp1); /输出地震数据 /for(i=0;im;i+) /printf(%d %fn,i,xi); /低通滤波器h(t) for(i=-50;i0;i+) hi+50=(float)(sin(2*pi*70*i*dt)/(pi*i*dt); for(i=1;i=50;i+) hi+50=(float)(sin(2*pi*70*i*dt)/(pi*i*dt); h50=140; /褶积 conv(x,100,h,101,y,200); /输出低通滤波后的值y(t) for(i=50;i150;i+) fprintf(fp2,%fn,yi); /带通滤波器h(t) for(i=-50;i0;i+) hi+50=(float)(sin(2*pi*80*i*dt)-sin(2*pi*10*i*dt)/(pi*i*dt); for(i=1;i=50;i+) hi+50=(float)(sin(2*pi*80*i*dt)-sin(2*pi*10*i*dt)/(pi*i*dt); h50=140; /褶积 conv(x,m,h,n,y,l); /输出带通滤波后的值y(t) for(i=50;i150;i+) fprintf(fp2,%fn,yi); fclose(fp2);2. 快变滤波F语言 PROGRAM EXAM2PARAMETER (PI=3.1415926535)PARAMETER (DT=0.002)COMPLEX C(128),D(128),Y(128)DIMENSION X(128),H(128) OPEN(1,FILE=INPUT1.DAT,FORM=FORMATTED,STATUS=OLD) READ(1,*)(X(I),I=1,100)DO 10 I=101,128 X(I)=0.0 10 CONTINUE DO 20 I=1,128 C(I)=CMPLX(X(I),0.0) 20 CONTINUECALL FFT(128,C,1.0)DF=1.0/(DT*128)DO 30 I=1,64 F=DF*I IF(F.LE.70)THEN H(I)=1.0 ELSE H(I)=0.0 ENDIF30CONTINUE DO 40 I=65,128 F=I*DF H(I)=H(129-I) 40 CONTINUEDO 50 I=1,128 Y(I)=CMPLX(H(I),0.0) D(I)=C(I)*Y(I) 50 CONTINUE CALL FFT(128,D,-1.0)OPEN(2,FILE=LOWPASS.DAT,FORM=FORMATTED) DO 110 I=1,128X(I)=REAL(D(I) WRITE(2,*)X(I) 110 CONTINUECLOSE(2)DO 60 I=1,64 F=I*DF IF(F.GE.10).AND.(F.LT.80) THEN H(I)=1.0 ELSE H(I)=0.0 ENDIF 60 CONTINUE DO 70 I=65,128 H(I)=H(129-I) 70 CONTINUE DO 80 I=1,128 Y(I)=CMPLX(H(I),0.0) D(I)=C(I)*Y(I) 80 CONTINUE CALL FFT(128,D,-1.0)OPEN(3,FILE=BANDPASS.DAT,FORM=FORMATTED)DO 90 I=1,128X(I)=REAL(D(I)WRITE(3,*)X(I) 90 CONTINUE CLOSE(3)CLOSE(1)END C语言/得到两个文件，低通滤波结果y21（一道128点），带通滤波y22（一道128点）#include#include #include #define pi 3.1415926void main() void fft(double sr,double sx,int m0,int inv); FILE *fp1,*fp2,*fp3; float dt,df,*h,*s,f; double *sr,*sx,*yr,*yx; int i,m,kc; dt=0.002; m=128; df=1/(m*dt); kc=70/df; s=(float *)calloc(128,sizeof(float); sr=(double *)calloc(128,sizeof(double); sx=(double *)calloc(128,sizeof(double); h=(float *)calloc(128,sizeof(float); yr=(double *)calloc(128,sizeof(double); yx=(double *)calloc(128,sizeof(double); fp1=fopen(INPUT1.DAT,r); fp2=fopen(y21.dat,w); fp3=fopen(y22.dat,w); /地震数据x(t) for(i=0;i100;i+) fscanf(fp1,%f,&si); fclose(fp1); for(i=0;i100;i+) sri=(double)si;printf(%fn,sri); for(i=100;i128;i+)sri=0.0; for(i=0;i128;i+)sxi=0.0; fft(sr,sx,7,1); /低通滤波器h(t) for(i=0;i128;i+)hi=0; for(i=0;i=kc;i+) hi=1.0; for(i=(127-kc);i128;i+) hi=1.0; /频域滤波 for(i=0;i128;i+) yri=sri*hi; yxi=sxi*hi; /反变换到时域 fft(yr,yx,7,-1); /输出 for(i=0;i128;i+) fprintf(fp2,%fn,yri); fclose(fp2); /带通滤波器h(t) for(i=0;i10/df)&(i=80/df) hi=1.0; else hi=0; for(i=128/2+1;i128;i+) hi=h127-i; /频域滤波 for(i=0;i128;i+) yri=sri*hi; yxi=sxi*hi; /反变换到时域 fft(yr,yx,7,-1); /输出 for(i=0;i128;i+) fprintf(fp3,%fn,yri); fclose(fp3);3、褶积滤波和递归滤波的比较F语言（1）褶积滤波 PROGRAM EXAM31DIMENSION X(50),H(20),Z(69),C(39),A(20),Y(88)OPEN(1,FILE=INPUT3.DAT,FORM=FORMATTED,STATUS=OLD)READ(1,*)(X(I),I=1,50)DATA H /1.0,5.254,11.6,13.7,8.47,0.775,-3.325,-2.764, $-0.364,1.099,0.97,0.15,-0.37,-0.344,-0.06,-0.126,0.122, $0.025,-0.042,-0.043/cccc 非零相位褶积滤波CALL CON(X,H,Z,50,20,L)OPEN(2,FILE=CON_FILTER.DAT,FORM=FORMATTED)DO 10 I=1,69WRITE(2,*)Z(I)10 CONTINUECLOSE(2)CCCC 零相位褶积滤波 DATA A /-0.043,-0.042,0.025,0.122,-0.126,-0.06,-0.344,-0.37, $0.15,0.97,1.099,-0.364,-2.764,-3.325,0.775,8.47,13.7,11.6, &5.254,1.0/ CALL CON(H,A,C,20,20,L)CALL CON(X,C,Y,50,39,NN)OPEN(3,FILE=ZERO_CON.DAT,FORM=FORMATTED)WRITE(3,*)(Y(I),I=1,88)CLOSE(3)CLOSE(1)END（2）递归滤波 PROGRAM EXAM32DIMENSION X(50),Y(50),A(5),B(4)OPEN(1,FILE=INPUT3.DAT,FORM=FORMATTED,STATUS=OLD)READ(1,*)(X(I),I=1,50)CCC 非零相位递归滤波 DO 10 I=1,50 Y(I)=0.010 CONTINUE DATA A /1,4,6,4,1/DATA B /-1.254,0.987,-0.341,0.0524/DO 20 I=1,50 DO 30 J=1,5 IF(I+1-J).GT.0)THEN Y(I)=Y(I)+A(J)*X(I-J+1) ELSE ENDIF 30 CONTINUE DO 40 K=1,4 IF(I-K).GT.0)THEN Y(I)=Y(I)-B(K)*Y(I-K) ELSE ENDIF40 CONTINUE 20 CONTINUE OPEN (2,FILE=RECURE_FILTER.DAT,FORM=FORMATTED)WRITE(2,*)(Y(I),I=1,50)CLOSE(2)CCCC 零相位递归滤波 DO 50 I=1,49 X(I)=0.050 CONTINUE DO 60 I=50,1,-1 DO 70 J=1,5 IF(I+J-1).LE.50)THEN X(I)=X(I)+A(J)*Y(I+J-1) ELSE ENDIF70 CONTINUE DO 80 K=1,4 IF(I+K).LE.50)THEN X(I)=X(I)-B(K)*X(I+K) ELSE ENDIF80 CONTINUE60 CONTINUE OPEN (3,FILE=ZERO_RECURE.DAT,FORM=FORMATTED)WRITE(3,*)(X(I),I=1,50)CLOSE(3)CLOSE(1)ENDC语言（1） 褶积滤波/自相关函数求出的结果和用褶积求出的结果不一样/输出两个文件，y311zheji（一道69点），y312zhejizero（一道88点）#include#include #include #define pi 3.1415926void main() void conv(float x,int m,float h,int n,float y,int l); void corr(float x,float y,int n); FILE *fp1,*fp2,*fp3; float dt; float *x,*y,*h1; float h20=1.0,5.254,11.6,13.7,8.47,0.775,-3.325,-2.764,-0.364,1.009,0.97,0.15,-0.37,-0.344,-0.06,-0.126,0.122,0.025,-0.042,-0.043; float h020=-0.043,-0.042,0.025,0.122,-0.126,-0.06,-0.344,-0.37,0.15,0.97,1.099,-0.364,-2.764,-3.325,0.775,8.47,13.7,11.6,5.254,1.0; int i,m,n,l; dt=0.002; m=50; n=20; l=m+n-1; x=(float *)calloc(50,sizeof(float); y=(float *)calloc(l,sizeof(float); h1=(float *)calloc(39,sizeof(float); fp1=fopen(INPUT3.DAT,r); fp2=fopen(y311zheji.dat,w); fp3=fopen(y312zhejizero.dat,w); /地震数据 for(i=0;im;i+) fscanf(fp1,%f,&xi); fclose(fp1); /输出地震数据 /for(i=0;im;i+) /printf(%d %fn,i,xi); /非零相位褶积滤波 conv(x,m,h,n,y,l); for(i=0;il;i+) fprintf(fp2,%fn,yi); fclose(fp2); /零相位褶积滤波 /滤波器h1(t) /corr(h,h1,20); conv(h,20,h0,20,h1,39); /零相位褶积滤波 conv(x,m,h1,39,y,88); for(i=0;i88;i+) fprintf(fp3,%fn,yi/0.002); fclose(fp3);/线性褶积子程序void conv(float x,int m,float h,int n,float y,int l) int k,i; for(k=0;kl;k+) yk=0.0; for(i=0;i=0&k-i=n-1)yk=yk+xi*hk-i*0.002; （2） 递归滤波/运行速度很慢？经正向递归滤波（从地震记录得到头部对整张记录递推），再经反向递归滤波（从输出记录的尾部递推），即得到零相位滤波结果/输出两个文件，y321（一道50点），y322zero（一道50点）#include#include #include void main() FILE *fp1,*fp2,*fp3; float dt=0.002; float *x,*y,*h; float a5=1,4,6,4,1; float b5=0,-1.254,0.987,-0.341,0.0524; int i,j, m=50; x=(float *)calloc(m,sizeof(float); y=(float *)calloc(m,sizeof(float); fp1=fopen(INPUT3.DAT,r); fp2=fopen(y321.dat,w); fp3=fopen(y322zero.dat,w); /地震数据 for(i=0;im;i+) fscanf(fp1,%f,&xi); fclose(fp1); /输出地震数据 /for(i=0;im;i+) /printf(%d %fn,i,xi); y0=x0; y1=x1; y2=x2; y3=x3; /非零相位递归滤波 for(i=4;i50;i+) yi=0.0; for(j=0;j5;j+)yi+=aj*xi-j-bj*yi-j; for(i=0;i=0;i-) xi=0.0;for(j=0;j5;j+)xi+=aj*yi+j-bj*xi+j; for(i=0;im;i+) fprintf(fp3,%fn,xi); fclose(fp3);4. 设计高通滤波因子 F语言 PROGRAM EXAM4PARAMETER(PI=3.1415926535)PARAMETER(DT=0.004)COMPLEX C(101)DIMENSION H(101),D(101)DO 10 I=-50,50 J=I+51 IF(I.NE.0)THEN H(J)=0.0-SIN(2*PI*30*I*DT)/(PI*I*DT) ELSE H(J)=-59 ENDIF10 CONTINUEDO 20 I=1,101C(I)=CMPLX(H(I),0.0)20 CONTINUECALL DFT(101,C,1.0)DO 30 I=1,101D(I)=ABS(C(I)30 CONTINUE OPEN(1,FILE=HIGHPASS.DAT,FORM=FORMATTED)WRITE(1,*) (D(I),I=1,101) CLOSE(1)ENDC语言/输出1道，样点数101#include#include #include #define pi 3.1415926void main() void dft(float x,float XR,float XI,int N); FILE *fp; float dt=0.004; float *xr,*xi,*h; int i,n=101; xr=(float *)calloc(101,sizeof(float); xi=(float *)calloc(101,sizeof(float); h=(float *)calloc(101,sizeof(float); fp=fopen(y4.dat,w); /高通滤波器h(t) for(i=-50;i0;i+) hi+50=(float)(-sin(2*pi*30*i*dt)/(pi*i*dt); for(i=1;i=50;i+) hi+50=(float)(-sin(2*pi*30*i*dt)/(pi*i*dt); h50=1/dt-60; /变换到频率域 dft(h,xr,xi,n); /输出 for(i=0;i101;i+) fprintf(fp,%fn,sqrt(xri*xri+xii*xii); fclose(fp);5. 频谱分析F语言 PROGRAM EXAM5 PARAMETER (PI=3.1415926535)PARAMETER (DT=0.004)COMPLEX C(128)DIMENSION X(128),A(128)CCCCCCCCCCCCCCCCCCCC周期波形CCCCCCCCCCCCCCCCCCCCCCCCDO 110 I=1,128 X(I)=SIN(2.0*PI*I/64)110 CONTINUE DO 120 I=1,128 C(I)=CMPLX(X(I)120CONTINUE CALL FFT(128,C,1.0)DO 130 I=1,128A(I)=ABS(C(I)*DT130CONTINUEOPEN(1,FILE=AMP1.DAT,FORM=FORMATTED)WRITE(1,*)(A(I),I=1,128)CLOSE(1)CCCCCCCCCCCCCCCCCC脉冲波形CCCCCCCCCCCCCCCCCCCCCCCCCCCC DO 210 I=1,128 IF(I.EQ.1.0) THEN X(I)=1.0 ELSE X(I)=0.0 ENDIF210 CONTINUE DO 220 I=1,128 C(I)=CMPLX(X(I)220CONTINUE CALL FFT(128,C,1.0)DO 230 I=1,128A(I)=ABS(C(I)*DT230CONTINUEOPEN(2,FILE=AMP2.DAT,FORM=FORMATTED)WRITE(2,*)(A(I),I=1,128)CLOSE(2)CCCCCCCCCCCCCCCCC地震波CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC OPEN(3,FILE=WAVE.DAT,FORM=FORMATTED,STATUS=OLD)READ(3,*)(X(I),I=1,128)DO 320 I=1,128 C(I)=CMPLX(X(I)320CONTINUE CALL FFT(128,C,1.0)DO 330 I=1,128A(I)=ABS(C(I)*DT330CONTINUEOPEN(4,FILE=AMP3.DAT,FORM=FORMATTED) WRITE(4,*)(A(I),I=1,128)CLOSE(4)CLOSE(3)ENDC语言/振幅谱（sqrt(xri*xri+xii*xii)）/输出3道，样点数为：128,128,128，分别为周期波形，脉冲波形，地震波的振幅谱#include#include #include #define pi 3.1415926void main() void fft(float sr,float sx,int m0,int inv); FILE *fp,*fpp/*,*fp1,*fp2,*fp3*/; float dt=0.004,*xr,*xi; int i,m=128; xr=(float *)calloc(128,sizeof(float); xi=(float *)calloc(128,sizeof(float); fp=fopen(WAVE.DAT,r); fpp=fopen(y5.dat,w); /fp1=fopen(y51.dat,w); /fp2=fopen(y52.dat,w); /fp3=fopen(y53.dat,w); /周期波形 for(i=1;i=128;i+)xri-1=(float)sin(2*pi*i/64); for(i=1;i=128;i+)xii-1=0.0; fft(xr,xi,7,1); for(i=0;i128;i+) fprintf(fpp,%fn,sqrt(xri*xri+xii*xii); / fprintf(fp1,%fn,fabs(xri); / fclose(fp1); /脉冲波形 xr0=1; for(i=2;i=128;i+)xri-1=0.0; for(i=1;i=128;i+)xii-1=0.0; fft(xr,xi,7,1); for(i=0;i128;i+) fprintf(fpp,%fn,sqrt(xri*xri+xii*xii); / fprintf(fp2,%fn,xri); / fclose(fp2); /地震波 for(i=0;i128;i+) fscanf(fp,%f,&xri); fclose(fp); for(i=1;i=128;i+)xii-1=0.0; fft(xr,xi,7,1); for(i=0;i128;i+) fprintf(fpp,%fn,sqrt(xri*xri+xii*xii); / fprintf(fp3,%fn,fabs(xri); fclose(fpp);6. 分析补零对振幅谱的影响F语言 PROGRAM EXAM6PARAMETER (PI=3.1415926535)DIMENSION X(60),Y(64),Z(128)COMPLEX A(60),B(64),C(128)CCCCCCCCCCCCCCCCCCC周期波形CCCCCCCCCCCCCCCCCCCCCCDO 10 I=1,60 X(I)=SIN(2*PI*I/30) A(I)=CMPLX(X(I),0.0) B(I)=A(I) C(I)=A(I)10 CONTINUE DO 20 I=61,6420 B(I)=CMPLX(0.0,0.0)DO 30 I=61,12830 C(I)=CMPLX(0.0,0.0) CALL DFT(60,A,1.0)CALL FFT(64,B,1.0)CALL FFT(128,C,1.0)DO 40 I=1,60 X(I)=ABS(A(I)40CONTINUE DO 50 I=1,6450 Y(I)=ABS(B(I)DO 60 I=1,12860 Z(I)=ABS(C(I) OPEN(1,FILE=AMP_SPE1.DAT,FORM=FORMATTED)WRITE(1,*)(X(I),I=1,60)WRITE(1,*)(Y(I),I=1,64)WRITE(1,*)(Z(I),I=1,128)CLOSE(1)CCCCCCCCCCCCCCCCC非周期波形CCCCCCCCCCCCCCCCCCCCCCCCCCCCC OPEN(1,FILE=WAVE.DAT,FORM=FORMATTED,STATUS=OLD)READ(1,*)(X(I),I=1,100) DO 70 I=1,60 A(I)=CMPLX(X(I),0.0) B(I)=A(I) C(I)=A(I)70 CONTINUEDO 80 I=61,6480 B(I)=CMPLX(0.0,0.0)DO 90 I=61,12890 C(I)=CMPLX(0.0,0.0) CALL DFT(60,A,1.0)CALL FFT(64,B,1.0)CALL FFT(128,C,1.0)DO 100 I=1,60 X(I)=ABS(A(I)100CONTINUE DO 110 I=1,64110 Y(I)=ABS(B(I)DO 120 I=1,128120 Z(I)=ABS(C(I) OPEN(2,FILE=AMP_SPE2.DAT,FORM=FORMATTED)WRITE(2,*)(X(I),I=1,60)WRITE(2,*)(Y(I),I=1,64)WRITE(2,*)(Z(I),I=1,128)CLOSE(2)ENDC语言/输出两个文件，y61为周期波形的（60,64,128）振幅谱，y62为地震波的（60,64,128）振幅谱#include#include #include #define pi 3.1415926void main() void fft(double sr,double sx,int m0,int inv); void dft(float x,float xr,float xi,int N); FILE *fp,*fp1,*fp2,*fp3; float dt,df,*x,*xr,*xi; double *sr,*sx; int i; x=(float *)calloc(128,sizeof(float); xr=(float *)calloc(128,sizeof(float); xi=(float *)calloc(128,sizeof(float); sr=(double *)calloc(128,sizeof(double); sx=(double *)calloc(128,sizeof(double); fp=fopen(WAVE.DAT,r); fp1=fopen(y61.dat,w); fp2=fopen(y62.dat,w); /周期波形dft for(i=1;i=60;i+)xi-1=(float)sin(2*pi*i/30); dft(x,xr,xi,60); for(i=0;i60;i+) fprintf(fp1,%fn,fabs(xri); /周期波形补零至64 for(i=1;i=60;i+)sri-1=(double)sin(2*pi*i/30); for(i=61;i=64;i+)sri-1=0.0; for(i=1;i=64;i+)sxi-1=0.0; fft(sr,sx,6,1); for(i=0;i64;i+) fprintf(fp1,%fn,sqrt(sri*sri+sxi*sxi); /周期波形补零至128 for(i=1;i=60;i+)sri-1=(double)sin(2*pi*i/30); for(i=60;i=128;i+)sri=0.0; for(i=1;i=128;i+)sxi-1=0.0; fft(sr,sx,7,1); for(i=0;i128;i+) fprintf(fp1,%fn,sqrt(sri*sri+sxi*sxi); fclose(fp1); /地震波 for(i=0;i60;i+) fscanf(fp,%f,&xi); fclose(fp); dft(x,xr,xi,60); for(i=0;i60;i+) fprintf(fp2,%fn,fabs(xri); /地震波补零至64 for(i=0;i60;i+) fscanf(fp,%f,&xi); fclose(fp); for(i=0;i60;i+)sri=(double)xi; for(i=60;i=63;i+)sri=0.0; for(i=0;i=63;i+)sxi=0.0; fft(sr,sx,6,1); for(i=0;i64;i+) fprintf(fp2,%fn,sqrt(sri*sri+sxi*sxi); /地震波补零至128 for(i=0;i60;i+) fscanf(fp,%f,&xi); fclose(fp); for(i=0;i60;i+)sri=(double)xi; for(i=60;i=127;i+)sri=0.0; for(i=0;i128;i+)sxi=0.0; fft(sr,sx,7,1); for(i=0;i128;i+) fprintf(fp2,%fn,sqrt(sri*sri+sxi*sxi); fclose(fp2);7. 线性褶积与循环褶积F语言PROGRAM EXAM7PARAMETER(PI=3.1415926535)PARAMETER(DT=0.002)DIMENSION X(100),X1(101),X2(200),H(101),H2(200), $Y(200),Y1(101) ,Y2(200)OPEN(1,FILE=INPUT1.DAT,FORM=FORMATTED,STATUS=OLD) READ(1,*)(X(I),I=1,100)DO 40 I=1,100 X1(I)=X(I) X2(I)=X(I) 40CONTINUE X1(101)=0.0DO 30 I=101,200 X2(I)=0.0 30CONTINUE DO 10 I=-50,50 IF(I.NE.0)THEN J=I+51 H(J)=SIN(2*PI*70*I*DT)/(PI*I*DT) ELSE J=I+51 H(J)=140 ENDIF 10 CONTINUE DO 50 I=1,101 50 H2(I)=H(I) DO 60 I=102,200 H2(I)=0.060 CONTINUE CALL CON(X,H,Y,100,101,200)CALL CIRCON(X1,H,Y1,101)CALL CIRCON(X2,H2,Y2,200) OPEN(2,FILE=CONV.DAT,FORM=FORMATTED) WRITE(2,*)(Y(I),I=1,200)WRITE(2,*)(Y1(I),I=1,101)WRITE(2,*)(Y2(I),I=1,200) CLOSE(2)CLOSE(1)ENDC语言/输出3道，样点数为：200,101,200#include#include #include #define pi 3.1415926void main() void conv(float x,int m,float h,int n,float y,int l); void cir_conv(float x,float h,float y,int n); FILE *fp1,*fp2; float dt; float *x,*h,*y; int i,m,n,l; dt=0.002; m=