polar_plot.m

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

function polar_plot(z);
%POLAR_PLOT: Plot zeros of a wavelet in polar coordinates
%
%  polar_plot(z)
%
%  IN   z: zeros of the wavelet, they can be computed
%          with zeros_wav.m
%
%  NOTE: some zeros migh end up outside the plot (check axis)
%
%  Example:    w = conv([2 -1 0.6+i*1 ],[1 -2 0+1.1*i] );
%              z = zeros_wav(w);
%              figure(1); polar_plot(z);
%
%  SeismicLab
%  Version 1
%
%  written by M.D.Sacchi, last modified March 4, 1998.
%  [e-mail]
%
%  Copyright (C) 1999 Seismic Processing and Imaging Group
%                     Department of Physics
%                     The University of Alberta
%
 

 
x=cos(0:0.1:2*pi); y = sin(0:0.1:2*pi);
 
plot(x,y);hold on;plot(real(z),imag(z),'s'); 
axis([-4.5 4.5 -4.5 4.5]); axis off
hold off;
return