数字处理实习报告.doc

附1：C语言和F语言程序 褶积滤波 F语言 PROGRAM EXAM1 PARAMETER(PI=3.1415926535) PARAMETER(DT=0.002) DIMENSION H(-50:50),X(100),Z(200) REAL X,H,Z INTEGER I OPEN(1,FILE=INPUT1.DAT,FORM=FORMATTED,STATUS=OLD) READ(1,*)(X(I),I=1,100) DO 10 I=-50,50 IF(I.NE.0)THEN H(I)=SIN(2*PI*70*I*DT)/(PI*I*DT) ELSE H(I)=140 ENDIF 10 CONTINUE CALL CON(X,H,Z,100,101,200) OPEN(2,FILE=CONLOW.DAT,FORM=FORMATTED) WRITE(2,*) Z DO 20 I=-50,50 IF(I.NE.0)THEN H(I)=(SIN(2*PI*80*I*DT)- SIN(2*PI*10*I*DT))/(PI*I*0.002) ELSE H(I)=140 ENDIF 20 CONTINUE 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点） #includestdio.h #include math.h #include stdlib.h #define pi 3.1415926 void 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,x[i]); fclose(fp1); //输出地震数据 //for(i=0;im;i++) //printf(%d %f\n,i,x[i]); //低通滤波器h(t) for(i=-50;i0;i++) h[i+50]=(float)(sin(2*pi*70*i*dt)/(pi*i*dt)); for(i=1;i=50;i++) h[i+50]=(float)(sin(2*pi*70*i*dt)/(pi*i*dt)); h[50]=140; //褶积 conv(x,100,h,101,y,200); //输出低通滤波后的值y(t) for(i=50;i150;i++) fprintf(fp2,%f\n,y[i]); //带通滤波器h(t) for(i=-50;i0;i++) h[i+50]=(float)((sin(2*pi*80*i*dt)-sin(2*pi*10*i*dt))/(pi*i*dt)); for(i=1;i=50;i++) h[i+50]=(float)((sin(2*pi*80*i*dt)-sin(2*pi*10*i*dt))/(pi*i*dt)); h[50]=140; //褶积 conv(x,m,h,n,y,l); //输出带通滤波后的值y(t) for(i=50;i150;i++) { fprintf(fp2,%f\n,y[i]); } fclose(fp2); } 快变滤波 F语言 PROGRAM EXAM2 PARAMETER (PI=3.1415926535) PARAMETER (DT=0.002) COMPLEX C(128),D(128),Y(128)