ERA matlab算法

function [ABCD]=era（hnNTsdef）;

% Eigensystem Realization Algorithm （ERA）
%
% Author: Samuel da Silva - UNICAMP
% e-mail: samsilva@fem.unicamp.br
% Date: 2006/10/20
% 
% [ABCD]=era（hnNTsdef）;
% 
% Inputs:
%    h: discrete-time impulse response
%    n: order of the system
%    N: number of samples to assembly the Hankel matrix
%    Ts: sample time
%    def: if = 1: the output will be the discrete-time state-space model
%         if = 2: the output will be the continuous-time state-space model
%          
% Otputs:
%    [ABCD]: state-space model
%  
% Note: For now it works to SISO systems and it is necessary the control toolbox
%
% References: Juang J. N. and Phan M. Q. “Identification and Control of
% Mechanical Systems“ Cambridge University Pre