zeros_wav.m

著名的seismiclab的代码 是地震学研究人员必备的工具

function [z] = zeros_wav(w);
%ZEROS_WAV: computes the zeros of a wavelet
%
%  [z] = zeros(w);
%
%  IN   w: a wavelet;
%
%  OUT  z: complex zeros
%
%  Example:  
%            w = [2 1 2 2.3 3 -10,1]; 
%            z = zeros_wav(w);
%            polar_plot(z); % polar_plot.m is also in SeismicLab.
%
%  SeismicLab
%  Version 1
%
%  written by M.D.Sacchi, last modified December 10,  1998.
%  [e-mail]
%
%  Copyright (C) 1998 Seismic Processing and Imaging Group
%                     Department of Physics
%                     The University of Alberta 


[N,M] = size(w);
if N == 1; wr = fliplr(w); z = roots(wr); end;
if M == 1; wr = flipud(w); z = roots(wr); end;
 
return