minus.m

地震、测井方面matlab代码,解释的比较详细

function ds=minus(i1,i2)
% Function adds a constant or matrix to traces of seismic dataset 
% Written by: E. R.: September 11, 2005
% Last updated:

if isstruct(i1)  &  strcmp(i1.type,'seismic')  &  isnumeric(i2)
   ds=i1;
   sz=size(i2);
   [nsamp,ntr]=size(i1.traces);
   if prod(sz) == 1  |  all(sz == [nsamp,ntr])
      ds.traces=ds.traces-i2;

   else
      if all(sz == [1,ntr])
         for ii=1:nsamp
            ds.traces(ii,:)=ds.traces(ii,:)-i2;
         end
      
      elseif all(sz == [nsamp,1])
         for ii=1:ntr
            ds.traces(:,ii)=ds.traces(:,ii)-i2;
         end
      
      else
         error('Operator "-" is not defined for this size of scaler.')
      
      end
   end   
         

elseif  isstruct(i2)  &  strcmp(i2.type,'seismic')  &  isnumeric(i1)
   ds=i2;
   sz=size(i1);
   [nsamp,ntr]=size(i2.traces);
   if prod(sz) == 1   |  all(sz == [nsamp,ntr])
      ds.traces=i1-ds.traces;

   else
      
      if all(sz == [1,ntr])
         for ii=1:nsamp
            ds.traces(ii,:)=i1-ds.traces(ii,:);
         end
      
      elseif all(sz == [nsamp,1])
         for ii=1:ntr
            ds.traces(:,ii)=i1-ds.traces(:,ii);
         end

      else
         error('Operator "-" is not defined for this size of scaler.')
      
      end
   end   
         

else
   error('Operator "-" is not defined for these arguments.')
end