swavelet3d.m

一种新的时频分析方法的matlab源程序。

echo on;
% --- The program SWAVELET3D computes and visualizes a 3-D Wavelet
% of the specific data set.
%
% Data loaded is:	'p71301.dat', 'p713010.dat','p713020.dat',
%	'p713030.dat', 'p713040.dat', 'p713050.dat','p713060.dat',
%	'p713070.dat', 'p713080.dat', 'p713090.dat','p713100.dat'.
% See marked lines in the code to change cases.
%
% Calling sequence-
% swavelet3d
%
% N. E. Huang (NASA GSFC)		14 Feb. 2001
% --- swavelet3d.m --- Version  14.Feb.2001 ---------------------

clear;						% Start Fresh

load p71301.dat; 				% Load Data
load p713010.dat; 	
load p713020.dat; 					
load p713030.dat; 
load p713040.dat; 	
load p713050.dat;
load p713060.dat; 
load p713070.dat; 	
load p713080.dat; 
load p713090.dat; 	
load p7130100.dat;

w1=wt2(p71301, 5.5, 4, 512, 100);	  % Compute the Wavelet for each line
w10=wt2(p713010, 5.5, 4, 512, 100);
w20=wt2(p713020, 5.5,  4, 512, 100);
w30=wt2(p713030, 5.5,  4, 512, 100);
w40=wt2(p713040, 5.5,  4, 512, 100);
w50=wt2(p713050, 5.5,  4, 512, 100);
w60=wt2(p713060, 5.5,  4, 512, 100);
w70=wt2(p713070, 5.5,  4, 512, 100);
w80=wt2(p713080, 5.5, 4, 512, 100);
w90=wt2(p713090, 5.5,  4, 512, 100);
w100=wt2(p7130100, 5.5,  4, 512, 100);

% Construct the 3-d Wavelt; the axies and the Wavelet volume

wx=linspace(0, 26.54, 170)';

k=linspace( 512/(26.54*4), 1/26.54,100)';



W=zeros(100,170, 11);			% Make & Assemble Data Volume

W(:, :, 1)=w1;
W(:, :, 2)=w10;
W(:, :, 3)=w20;
W(:, :, 4)=w30;
W(:, :, 5)=w40;
W(:, :, 6)=w50;
W(:, :, 7)=w60;
W(:, :, 8)=w70;
W(:, :, 9)=w80;
W(:, :, 10)=w90;
W(:, :, 11)=w100;



save Wvlt W wx k;

tt=linspace(15, 25, 11)';

% --- Change this Value -----> x.xxx <---------- Change this Value -- %

p = patch(isosurface(wx,k,tt,W, 0.01));

% ------------------------------------------------------------------- %
isonormals(wx,k, tt, W, p)
set(p, 'FaceColor', 'cyan', 'EdgeColor', 'none');
daspect([.5 .2 0.25])
campos([-45 10 90]);
lighting phong;
grid off;
camlight left

% --- Change the Title ---> Case   Frames  Line No. Contour Level --- %
title( 'X Distance: 0-26 cm; Y Distance: 15-25 cm; Wavenumber: 0-4.5 1/cm')
toptitle('3D Wavelet Spectrum: NEH-07 F130 L193 to 492 C.01')
% ------------------------------------------------------------------- %

% --- swavelet3d.m Ends Normally --- %