用Levy或hankel矩阵等方法系统辨识程序

系统辨识的方法：包括面积法1、面积法2、Levy法、Hankel矩阵、脉冲响应

clear;
%给定测试系统（分子、分母系数）
num=[17.5 7.5 1];
den=[4 15 17.5 7.5 1];
%--------------产生辨识所用阶跃响应数据h----------------
T0=0.01;                    %采样时间间隔1s           
t=0:T0:100;     
L=length（t）;                %数据长度
h=（step（numdent））.‘;      %根据传递函数求系统的阶跃响应
figure（1）;
stem（h‘.‘）;

%-----------------------开始辨识-------------------------
K=h（L）;                  %系统增益
ht=h/K;                  %无因次序列
n=length（den）-1;
m=length（num）-1;
%求各阶面积根据4.2.21式
for i=1:n+m
    temp1=0;
    for j=1:L
        temp1=temp1+（1-ht（j））*（-1*t（j））^（i-1）/factorial（i-1）*T0;
    end
    if i>1
        temp2=0;
        for ii=1:i-1
            temp=0;
            for j=1:L
                 temp=temp+（1-ht（j））*（-1*t（j））^（ii-1）/factorial（ii-1）*T0;
            end
            temp2=temp2+A（i-ii

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1523  2008-12-09 22:19  第二次作业\AreaI.m

     文件       1345  2009-09-19 11:36  第二次作业\AreaII.m

     文件        740  2009-09-17 20:40  第二次作业\Hankel.m

     文件        768  2008-12-04 16:32  第二次作业\HankelTestMain.m

     文件       2957  2009-09-19 20:39  第二次作业\Levy.m

     文件        711  2008-12-04 15:55  第二次作业\PulseTF.m

     文件        699  2008-12-04 16:32  第二次作业\PulseTFTestMain.m

     文件       1331  2008-12-08 11:15  第二次作业\SysFT.m

     目录          0  2009-10-17 09:23  第二次作业

----------- ---------  ---------- -----  ----

                10074                    9