sembucc.m

这是matlab在地球物理数据处理方面的源码

function semb = sembUCC(traces,t,x,v1,v2,n,window)

% function out = sembUCC(traces,t,x,v1,v2,n,window)
% This function calculates the unnormalized crosscorrelation type of 
% velocity spectra plot.
% traces = gather, t = timebase in seconds, x = offsets 
% v1 = initial velocity in m/s, v2 = final constant nmo velocity
% n = number of trial velocities to calculate.
% window = number of time samples of window over which coherency is measured.

[m,N] = size(traces);
delv = (v2-v1)/(n-1);
subplot(1,2,1)
imagesc(x,t,traces)
title('gather')

semb = zeros(m,n);
vel = v1:delv:v2;


for i = 1:n
   junk = zeros(m+window,N);
   edge = round(window/2);
	junk(edge+1:m+edge,:) = nmo(traces,t,x,vel(i)*ones(m,1));
	
   for k = 1:m
      semb(k,i) = 0.5*sum(sum(junk(k:k+window-1,:),2).^2-sum(junk(k:k+window-1,:).^2,2));     
    end % k  
end % i

subplot(1,2,2)
imagesc(vel,t,semb), xlabel('Stacking Velocity m/s')
ylabel('Time seconds')


out = 1;

return;