西南石油信号分析课程设计测井专业.docxVIP

f=35; N=-30:30; tn=N.*0.004; bn=(1-2.*pi.*pi.*f.*f.*tn.*tn).*exp(-pi.*pi.*f.*f.*tn.*tn); g=zeros(1,300); g(35)=0.4; g(80)=0.4; g(101)=-0.15; g(123)=-0.5; g(155)=0.35; g(203)=-0.1; g(227)=-0.2; g(257)=-0.2; s=conv(bn,g); s(1:30)=[]; s(301:330)=[]; n=rand(1,300)-0.5; x=s+n; p1=sum(s.*s); p2=sum(n.*n); p=p1/p2; fprintf(%d,p); subplot(511); plot(N.*0.004,bn);title(地震子波);xlabel(t);ylabel(b(t)); subplot(512); A=0:299; plot(A.*0.004,g);title(反射系数); xlabel(t);ylabel(g(t)); subplot(514); plot(A.*0.004,n);title(随机干扰); xlabel(t);ylabel(n(t)); subplot(513); plot(A.*0.004,s); title(合成地震记录); xlabel(t);ylabel(s(t)); subplot(515); plot(A.*0.004,x);title(地震记录); xlabel(t);ylabel(x(t)); N=260;b=100;A1=1;A2=0.7;A3=0.5;A4=0.8;f1=35;f2=45;f3=70;f4=100; i=0:N-1; ti=i.*0.004; xi=A1.*exp(-b.*ti.*ti).*sin(2.*pi.*f1.*ti)+A2.*exp(-b.*ti.*ti).*sin(2.*pi.*f2.*ti)+A3.*exp(-b.*ti.*ti).*sin(2.*pi.*f3.*ti)+A4.*exp(-b.*ti.*ti).*sin(2.*pi.*f4.*ti); k=0:N-1; wi=exp(-j*2*pi/N); ik=i*k; wiik=wi.^ik; xkk=abs(xi*wiik); figure(1); subplot(311); plot(i.*0.004,xi);title(原始信号得出的结果); xlabel(t);ylabel(x(t)); subplot(312); plot(k.*1./260./0.004,xkk); xlabel(f);ylabel(|X(f)|); subplot(313); stem(k.*1./260./0.004,xkk); xlabel(f);ylabel(|X(f)|); xx=zeros(1,100); Xi=[xi,xx]; XKK=abs(fft(Xi)); figure(2); subplot(311); A=0:359; plot(A.*0.004,Xi);title(后面补零100个得出的结果); xlabel(t);ylabel(x(t)); subplot(312); plot(A.*1./360./0.004,XKK); xlabel(f);ylabel(|X(f)|); subplot(313); stem(A.*1./360./0.004,XKK); xlabel(f);ylabel(|X(f)|); xxx=zeros(1,300); y=[xi,xxx]; Y=abs(fft(y)); figure(3); subplot(311); B=0:559; plot(B.*0.004,y); title(后面补零300个得出的结果); xlabel(t);ylabel(x(t)); subplot(312); plot(B.*1./560./0.004,Y); xlabel(f);ylabel(|X(f)|); subplot(313); stem(B.*1./560./0.004,Y); xlabel(f);ylabel(|X(f)|); Xii=[xx,xi]; XKKK=abs(fft(Xii)); figure(4); subplot(311); plot(A.*0.004,Xii);title(前面补零100个得出的结果); xlabel(t);ylabel(x(t)); subplot(312); plot(A.*1./360./0.004,XKKK); xlabel(f);ylabel(|X(f)|); subplot(313); stem(A.*1./360./0.004,