clear;
%給定測試系統（分子、分母系數）
num=[17.5?7.5?1];
den=[4?15?17.5?7.5?1];
%--------------產生辨識所用階躍響應數據h----------------
T0=0.01;????????????????????%采樣時間間隔1s???????????
t=0:T0:100;?????
L=length（t）;????????????????%數據長度
h=（step（numdent））.‘;??????%根據傳遞函數求系統的階躍響應
figure（1）;
stem（h‘.‘）;

%-----------------------開始辨識-------------------------
K=h（L）;??????????????????%系統增益
ht=h/K;??????????????????%無因次序列
n=length（den）-1;
m=length（num）-1;
%求各階面積根據4.2.21式
for?i=1:n+m
????temp1=0;
????for?j=1:L
????????temp1=temp1+（1-ht（j））*（-1*t（j））^（i-1）/factorial（i-1）*T0;
????end
????if?i>1
????????temp2=0;
????????for?ii=1:i-1
????????????temp=0;
????????????for?j=1:L
?????????????????temp=temp+（1-ht（j））*（-1*t（j））^（ii-1）/factorial（ii-1）*T0;
????????????end
????????????temp2=temp2+A（i-ii

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1523??2008-12-09?22:19??第二次作業\AreaI.m

?????文件???????1345??2009-09-19?11:36??第二次作業\AreaII.m

?????文件????????740??2009-09-17?20:40??第二次作業\Hankel.m

?????文件????????768??2008-12-04?16:32??第二次作業\HankelTestMain.m

?????文件???????2957??2009-09-19?20:39??第二次作業\Levy.m

?????文件????????711??2008-12-04?15:55??第二次作業\PulseTF.m

?????文件????????699??2008-12-04?16:32??第二次作業\PulseTFTestMain.m

?????文件???????1331??2008-12-08?11:15??第二次作業\SysFT.m

?????目錄??????????0??2009-10-17?09:23??第二次作業

-----------?---------??----------?-----??----

????????????????10074????????????????????9