科技创新中期检查报告中期成果

《科技创新中期检查报告中期成果》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、多次波压制的L1范数自适应相减方法研究中期成果：(成员：陈阳康，崔宇佳，金兆宇，孙佳传，张亦卓)初步程序框架：主函数：Clearfidl=fopen('data・bin','r*);fid2=fopen('multiple'z'r*);fid3=fopen(1primary1z1w1);shot=5;trace=201;sample=750;win=200;nf=15;sample_shift=15;cl(l:win/2)=0;c2(1:win/2)=0;fori=l:win/2;cl(i)=(win/2-i)/(win/2);c2(

2、i)=i/(win/2);endfork=l:shot;D0=zeros(sample,trace);M0=zeros(sample,trace);D0=fread(fidl,[sample,trace],1float');M0=fread(fid2,[sample,trace]z'float');%subplot(121)zimagesc(DO)%subplot(122)zimagesc(MO)D=D01;M=M01;M2=zeros(trace,sample);forj=l:trace;forkk=l:sample;if(kk+s

3、ample_shift<=sample)M2(jzkk)=M(jAkk+sample_shift);endendendM=M2;M=M/max(max(M));D=D/max(max(D));[DI]=match_win(win,cl,c2,Dztrace,sample,nf);elfsubplot(121),imagesc(DO);subplot(122),imagesc(DI!;pause(0.1);fwrite(fid3,DI*,1float1);fprintf(1%dshotsubtractiondone!1,k);e

4、nd%imagesc([D1zM1zDIT]);sgray(15);gridon子函数match_win：function[DI]=match_win(win,cl,c2,D,Mztrace,sample,nf)DI(1:trace,1:sample)=0;kl=l;k2=kl+win-l;while(k2<=sample)%fprintf(11Matchbetween%d--〉%dkl,k2);dd(:,1:win)=match(D(:,kl:k2)ZM(:,kl:k2),trace,winznf);if(kl==l)fork

5、k=kl:k2;if(kk<=win/2)DI(:zkk)=dd(:rkk);elseDl(:,kk)=cl(kk-win/2)*dd(:,kk);endendelseforkk=kl:k2;if(kksample)

6、%fprintf(12Matchbetween%d-->%d1,kl,sample);winl=sample-kl+l;if(nf>winl)nf=floor(winl/4);endddl(:,1:winl)=match(D(:,kl:sample)zM(:,kl:sample)ftrace,winl,nf);forkk=kl:sample;if(kk

7、dreturn;子函数match：function[Dl]=match(D,trace,sample,nf);forii=l:trace;fori=l:nf;forj=l:sample;if(j-i+l>0&&j-i+l<=sample);a(j,i)=M(ii,j-i+1);endendendAA=a1*a;b=a**D(iiz:)*;AA=AA+m=x(max(AA))*0-05*eye(nf);if(sum(AA)~=0)f=AAb;elsef=0;endy=a*f；if(sum(AA)~=0)Dl(ii,:)=D(ii,:)

8、-y(:)1;elseD1=D;endendreturn;效果图：左图为压制前的地震波时距曲线图；右图为压制后的地震波时距曲线图。仔细观察可以看出多次波明就受到了较人程度压制。