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