function?[ba]=Hankel（ngkT0）?
%用Hankel矩陣法求系統傳函；
%輸入參數：n：系統階數；gk：系統脈沖響應序列；T0：采樣間隔
%返回參數：b:?系統脈沖傳函分子系數；a:系統脈沖傳函分母系數

H=[];
for?i=1:n???????????????????%構造Hankel矩陣H
????for?j=1:n
????????temp（j）=gk（i+j）;????????
????end
????H=[H;temp];
end

g1=gk（n+2:2*n+1）*（-1）;
a=flipud（H\g1）;????????????%獲得脈沖傳函分母系數a

g2=gk（2:n+1）;
for?i=1:n?????????????????%根據a構造矩陣A
????for?j=1:n
????????if?i????????????A（ij）=0;
????????else?if?i==j
????????????????A（ij）=1;
????????????else?A（ij）=a（i-j）;
????????????end
????????end
????end
end
b=A*g2;??????????????????%獲得脈沖傳函分子系數b
%調整系數
b=b.‘;??????????????????
a=[1?a.‘];

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2019-05-09?00:38??系統辨識第一次作業\1st\
?????文件?????????740??2009-09-17?20:40??系統辨識第一次作業\1st\Hankel.m
?????文件?????????684??2019-05-09?00:30??系統辨識第一次作業\1st\HankelTestMain.m
?????文件???????40376??2019-05-09?00:38??系統辨識第一次作業\1st\impulse_si.mdl
?????文件???????47635??2019-05-09?00:14??系統辨識第一次作業\1st\ori_c&d_impulse_response.jpg
?????文件???????12147??2019-05-08?21:09??系統辨識第一次作業\1st\param.mat
?????文件?????????328??2019-05-08?23:21??系統辨識第一次作業\1st\test.m
?????文件???????18378??2019-05-08?19:43??系統辨識第一次作業\1st\untitled.slx
?????文件???????59592??2019-05-09?00:14??系統辨識第一次作業\1st\辨識后系統脈沖響應對比.jpg
?????文件???????56660??2019-05-09?00:20??系統辨識第一次作業\1st\辨識后系統脈沖響應對比0.01.jpg
?????文件???????57052??2019-05-09?00:23??系統辨識第一次作業\1st\辨識后系統脈沖響應對比0.1.jpg
?????文件???????60099??2019-05-09?00:25??系統辨識第一次作業\1st\辨識后系統脈沖響應對比0.23.jpg
?????文件???????59096??2019-05-09?00:26??系統辨識第一次作業\1st\辨識后系統脈沖響應對比t1.3.jpg
?????文件??????677672??2019-05-09?00:37??系統辨識第一次作業\脈沖響應法系統辨識.pptx