clear;clc;
%levy法，其實(shí)就是計(jì)算VSTU?和解一個(gè)行列式
%二階，改不改差不多
%需輸入的參數(shù)：w?Re?Im?L
%--------------生成測(cè)試數(shù)據(jù)w、Gw-------------
%給定測(cè)試系統(tǒng)（分子、分母系數(shù)）:???傳函??（1.18s+20）/（10.9s^2+371s+1）
num=[1.18?20];
den=[10.9?371?1];
T0=1;??????????????????????????????????????????????????????????????????????%采樣時(shí)間間隔1s
L=100;?????????????????????????????????????????????????????????????????????%數(shù)據(jù)長(zhǎng)度
k=1:T0:L;??????
w=2*pi/L*k;????????????????????????????????????????????????????????????????%頻率，0到2pi內(nèi)均勻取100份，最大是不是可以取更大？
%--------------獲得系統(tǒng)頻率特性--------------
for?i=1:L
????Gw（i）=（1.2j*w（i）+20）/（10*（j*w（i））^2+370j*w（i）+1）;??????????????????????%傳函參數(shù)為什么與頻響一樣？
end
figure（1）
plot（abs（Gw））;
figure（2）
stem（abs（Gw）‘.‘）;
%------------------用Levy法辨識(shí)---------------------
%系統(tǒng)分子、分母階次m、n
m=length（num）-1;n=length（den）-1;???????