clc;clear all; close all;
sample=512; %采样点
trace=64; %道数
f=35;
for m=1 : trace
% n=-255:1:256; %模型1
% n=-400+m:1:111+m; %模型2
%模型3
mm=round((m-32)^2/40);
n=-400+mm:1:111+mm;
C=randn(1,sample);
fn=(1-2*(0.001*n*pi*f).^2).*exp(-(0.001*n*pi*f).^2);
Q(m,:)=fn;
% N(m,:)=C;
B=C*0.03+fn;
A(m,:)=B;
% w=awgn(B(:,m),1,'measured');%添加SNR=1的白噪声
% A(:,m)=w;
end
snr=d_snr(Q,A); %理论
for m=1 : trace
figure(1);
plot(A(m,:)+m,1:512); hold on;
end
hold off; axis tight;
title('地震子波');
xlabel('道数');
ylabel('时间/ms');
for i = 1:sample
k(i)= (sum(A(:,i)))^2;
end
Es = sum(k(1:sample))/trace;
E = sum(sum(A.^2));
En = E - Es;
SNR = 10*log10(Es/En); %能量叠加法
评论0