地探学院论格式1.docVIP

record=zeros(1,180); w1=zeros(1,180); f=25; fs=1000; number=100; t=-number/2+1:number/2; t1=t+50; t2=t1+80; a(t1)=8*cos(2*(pi* f*t/ fs)).*exp(-25*t/fs); b(t1)= 1/6* a(t1); for i=1:100 w1(1,i)=b(1,i); end plot(t1,b(t1)); hold on; w2=zeros(1,180); aa(t2)=- 8*cos(2*(pi* f*t/ fs)).*exp(-25*t/fs); bb(t2)= 1/6* aa(t2); for i=1:100 w2(1,i+80)=bb(1,i+80); end plot(t2,bb(t2)); hold off; for i=1:180 record(1,i)=w1(1,i)+w2(1,i); end figure; plot(record); y=fft(record,180); mag=abs(y*2/180); f=(0:179)*fs/180; xp=angle((y)*180/pi); subplot(232); plot(f,xp); axis([0,500,-5,5]); xl