prony 工具箱in matlab

Prony Toolbox is a software tool in MATLAB which performs Prony analysis. Prony Toolbox (PTbox) is designed based on several considerations including data preprocessing, model order selection, model order selection criteria, signal subspace selection, signal and noise separation, root inspection and assessing residuals. The PTbox provides flexibility to compare and display analysis results simultaneously for several parameter variations.

function [iappaia_listtau_listomega_listSUB_INDenergyp]=applyprony（test_timetest_dataNSUB_Ncriteria_val）
% Implements MATLAB Signal Processing Tooltbox prony function 
% performprony GUI

test_t_increment=test_time（2）-test_time（1）;
fs=1/test_t_increment; % Sampling frequency=5 e8
[inumzidenz]=prony（test_dataNN）;
% Compute and plot the impulse response of the Prony approximating system
[iapptapp]=impz（inumzidenzlength（test_time）fs）;
[rpk]=residuez（inumzidenz）;
a_list=r（:）;
spoles=log（p（:））*fs;
tau_list=1./real（spoles）;
omega_list=imag（spoles）;
n_size=size（test_time1）; 
n=0:1:n_size;
for ct=1:N
    energy（ct）=abs（a_list（ct）^2）*（sum（abs（p（ct）.^n）.^2））;
end

if （criteria_val==1）
    [Mag ISort]=sort（abs（a_list））;
else
    [EnISort]=sort（energy）;
end
ISort=ISort（end:-1:1）;
FULL_IND=[ISort（:）]‘ ;
ai=zeros（size（test_time））;
SUB_IND=FULL_IND（1:SUB_N）;
for i=SUB_IND
   ai=ai+a_list（i）*exp（spoles（i）*tapp）;
end
ai=real（ai）;
% Remove bug PA should be same if modes are equal to model order
if（SUB_N==N）
    ai（1）=iapp（1）;
end
% Update the SUB_N for prony fit when no. of modes is not specified
if（SUB_N==0）
    SUB_IND=FULL_IND;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2596  2003-07-06 13:53  Prony Toolbox Download\about.html
     文件        1201  2003-06-14 18:27  Prony Toolbox Download\applyprony.m
     文件       52488  2003-05-22 15:29  Prony Toolbox Download\comparesessions.fig
     文件       24992  2003-07-06 14:13  Prony Toolbox Download\comparesessions.m
     文件        2988  2003-07-06 15:17  Prony Toolbox Download\demos.html
     目录           0  2003-08-09 17:22  Prony Toolbox Download\Demos\
     文件        3456  2003-07-14 01:10  Prony Toolbox Download\Demos\bbnoise_demo.mat
     文件         318  2003-07-14 01:10  Prony Toolbox Download\Demos\create_demo.m
     文件        3456  2003-07-14 01:03  Prony Toolbox Download\Demos\pure_demo.mat
     文件        3456  2003-07-14 01:07  Prony Toolbox Download\Demos\puretone_demo.mat
     文件       11947  2003-07-06 15:18  Prony Toolbox Download\design.html
     文件       22638  2002-10-31 21:25  Prony Toolbox Download\ecelogo.jpg
     文件         472  2003-07-06 22:17  Prony Toolbox Download\export_file.mat
     文件      105312  2003-07-06 18:40  Prony Toolbox Download\exportresults.fig
     文件       71116  2003-07-06 23:24  Prony Toolbox Download\exportresults.m
     文件         637  2003-07-07 14:46  Prony Toolbox Download\fft_analysis.m
     文件       12450  2003-07-06 15:19  Prony Toolbox Download\gui.html
     文件      235384  2003-04-27 01:17  Prony Toolbox Download\imageicons.mat
     文件        2389  2003-02-21 13:08  Prony Toolbox Download\mmv2struct.m
     目录           0  2003-08-09 17:22  Prony Toolbox Download\PCRdata\
     目录           0  2003-08-09 17:22  Prony Toolbox Download\PCRdata\60Ar\
     文件       40272  2003-05-15 16:42  Prony Toolbox Download\PCRdata\60Ar\60Ar_i100.mat
     文件       40272  2003-05-19 16:15  Prony Toolbox Download\PCRdata\60Ar\60Ar_i1000.mat
     文件       40272  2003-05-15 16:48  Prony Toolbox Download\PCRdata\60Ar\60Ar_i1500.mat
     文件       40272  2003-05-15 16:43  Prony Toolbox Download\PCRdata\60Ar\60Ar_i400.mat
     文件       40272  2003-05-15 16:44  Prony Toolbox Download\PCRdata\60Ar\60Ar_i700.mat
     文件       40272  2003-05-15 16:51  Prony Toolbox Download\PCRdata\60Ar\60Ar_v100.mat
     文件       40272  2003-05-15 16:50  Prony Toolbox Download\PCRdata\60Ar\60Ar_v1000.mat
     文件       40272  2003-05-15 16:49  Prony Toolbox Download\PCRdata\60Ar\60Ar_v1500.mat
     文件       40272  2003-05-15 16:51  Prony Toolbox Download\PCRdata\60Ar\60Ar_v400.mat
     文件       40272  2003-05-15 16:50  Prony Toolbox Download\PCRdata\60Ar\60Ar_v700.mat
............此处省略32个文件信息