Home > tidal_ellipse > example.m

example

PURPOSE ^

demonstrate how to use ap2ep and ep2ap

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

demonstrate how to use ap2ep and ep2ap

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 %demonstrate how to use ap2ep and ep2ap
0002 %
0003 Au=rand(4,3,2);         % so 4x3x2 multi-dimensional matrices are used for the
0004 Av=rand(4,3,2);         % demonstration.
0005 Phi_v=rand(4,3,2)*360;  % phase lags inputs are expected to be in degrees.
0006 Phi_u=rand(4,3,2)*360;
0007 
0008 figure(1)
0009 clf
0010 [SEMA ECC INC PHA w]=ap2ep(Au, Phi_u, Av, Phi_v, [2 3 1]);
0011 figure(2)
0012 clf
0013 [rAu rPhi_u rAv rPhi_v rw]=ep2ap(SEMA, ECC, INC, PHA, [2 3 1]);
0014 
0015 %check if ep2ap has recovered Au, Phi_u, Av, Phi_v
0016 max(max(max(abs(rAu-Au))))               %  = 9.9920e-16
0017 max(max(max(abs(rAv-Av))))               %  = 6.6613e-16
0018 max(max(max(abs(rPhi_u-Phi_u))))         %  = 4.4764e-13
0019 max(max(max(abs(rPhi_v-Phi_v))))         %  = 1.1369e-13
0020 max(max(max(max(abs(w-rw)))))            %  = 1.3710e-15
0021 % for the random realization I had, the differences are listed on the right
0022 % hand of the above column. What are yours?
0023 
0024 % The above example function calls have already plotted an ellipse for you.
0025 % To plot an ellipse separately, you may do
0026 %
0027 figure(3)
0028 clf
0029 plot(real(squeeze(w(2,3,1,:))), imag(squeeze(w(2,3,1,:))));
0030 
0031 %here squeeze is needed because w is a multiple dimensional array.
0032 
0033 
0034 %Zhigang Xu
0035 %Nov. 12, 2000

Generated on Wed 30-Nov-2005 15:40:57 by m2html © 2003