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;???????