过程辨识中频率响应的levy法

clear;clc;
%levy法，其实就是计算VSTU 和解一个行列式
%二阶，改不改差不多
%需输入的参数：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
L=100;                                                                     %数据长度
k=1:T0:L;      
w=2*pi/L*k;                                                                %频率，0到2pi内均匀取100份，最大是不是可以取更大？
%--------------获得系统频率特性--------------
for i=1:L
    Gw（i）=（1.2j*w（i）+20）/（10*（j*w（i））^2+370j*w（i）+1）;                      %传函参数为什么与频响一样？
end
figure（1）
plot（abs（Gw））;
figure（2）
stem（abs（Gw）‘.‘）;
%------------------用Levy法辨识---------------------
%系统分子、分母阶次m、n
m=length（num）-1;n=length（den）-1;