levy法编制辨识程序

levy法编制辨识程序，基本精确辨识了系统各个参数，因为误差准则的原因，系数略有差别。

clear;
%w Re Im L
%--------------生成测试数据w、Gw-------------
%给定测试系统（分子、分母系数）:（1.18s+20）/（10.9s^2+371s+1）
num=[1.18 20];
den=[10.9 371 1];
T0=1;                       %采样时间间隔1s          
k=1:T0:100;      
L=length（k）;                %数据长度
w=2*pi/L*k;
%获得系统频率特性
for i=1:L
    Gw（i）=（1.2j*w（i）+20）/（10*（j*w（i））^2+370j*w（i）+1）;
end
figure（2）
stem（abs（Gw）‘.‘）;
%------------------用Levy法辨识---------------------
%系统分子、分母阶次m、n
m=length（num）-1;
n=length（den）-1;
N=2*n+1;              %j=1:N
Re=real（Gw）;
Im=imag（Gw）;
%计算VSTU
for j=1:N
    V（j）=0;
    for i=1:L
        V（j）=V（j）+w（i）^（j-1）;
    end
end
for j=1:N
    S（j）=0;
    for i=1:L
        S（j）=S（j）+w（i）^（j-1）*Re（i）;
    end
end
for j=1:N
    T（j）=0;
    for i=1:L
        T（j）=T（j）+w（i）^（j-1）*Im（i）;
    end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2945  2011-09-21 08:33  Levy.m