dpsssave.m

多频道功率谱计算程序

function dpsssave(N,W) 
% A simple utility function to calculate dpss and save them to disk 
% with filename and format as used by the various mtm routines. 
%  
% Alternatively, a script something like the following could be 
% used to totally automate the process - run this one overnight 
% if you are going to compute many large dpss! 
% 
% cd e:\matlab\mtm    % Substitute your desired directory here. 
% [E,V]=dpsscalc(256,2);   save s256_2.mat 
% [E,V]=dpsscalc(256,5/2); save s256_25.mat 
% [E,V]=dpsscalc(256,3);   save s256_3.mat 
% [E,V]=dpsscalc(256,7/2); save s256_35.mat 
% [E,V]=dpsscalc(512,2);   save s512_2.mat 
% [E,V]=dpsscalc(512,5/2); save s512_25.mat 
% [E,V]=dpsscalc(512,3);   save s512_3.mat 
% [E,V]=dpsscalc(512,7/2); save s512_35.mat 
% 
% etc., etc. 
% 
% Written by Eric Breitenberger, version date 10/1/95. 
% Please send comments and suggestions to [e-mail] 
% 
 
 
dir=input('Directory where DPSS are to be stored: ', 's'); 
eval(['cd ' dir]); 
 
if rem(W,1)==0 
  filename=['s' int2str(N) '_' int2str(W) '.mat'];  
else 
  filename=['s' int2str(N) '_' int2str(10*W) '.mat'];  
end 
 
disp('Calculating the DPSS...') 
 
eval(['save ' filename])