%pq=3peqe=5M=30;
N?=?500;%樣本點數
B?=?1;%分母項，確定為AR模型
A?=[1-1.30.86?-0.676];%分子項，可根據它求出實際的功率譜
yangben=filter（B?A0.1*randn（1N））;%濾波產生樣本

r=zeros（1N）;%為總體最小二乘演化公式z，詳見式3.4.45
pp=0;%有效秩初始化
for?p=1:35
????for?i=1:N-p
????????r（p）=r（p）+yangben（i）*yangben（i+p）;%混入噪聲時的自相關函數
????end
????r（p）=1/（N-p）*r（p）;
end
for?l=1:30
????re（l:）=[r（l+5）?r（l+4）?r（l+3）?r（l+2）?r（l+1）?r（l）];%生成矩陣RE即為3.4.43
end
[USV]?=svd（re）;%奇異值分解
for?i=1:6
????s=S/S（11）;%歸一化奇異值
end
for?i=1:6
????if?s（ii）>0.05
????????pp=pp+1;%確定有效秩
????end
end
pp
ss1=zeros（pp+1）;
for?tt=1:6-pp
????for?j=1:pp
????????ss1=ss1+S（jj）^2*V（tt:tt+ppj）*V（tt:tt+ppj）‘;%式3.4.56，求出矩陣SP
????end
end
a=inv（ss1）;
b=a（1:pp+11）/a（11）%式3.4.59，求出未知參數的最小二乘估計
%?bfi=b/b（1）
[h1?w1]?=?freqz（B?A?100）;
[h2?w2]?=?freqz（1?b‘?100）;
plot（w1?abs（h1）?‘r‘??w2?abs（h2））;%畫圖，比較最小二乘估計與實際模型的差別。
legend（‘估計值‘‘實際值‘2）
%?zplane（BA）;

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

?????文件???????1052??2018-01-04?21:46??信號處理代碼\arma.m

?????文件????????845??2018-01-04?21:46??信號處理代碼\bandp.m

?????文件???????2041??2018-01-04?21:46??信號處理代碼\fmlin.m

?????文件????????283??2018-01-04?21:46??信號處理代碼\kalman_filter.m

?????文件???????2409??2018-01-04?21:47??信號處理代碼\mfiles\altes.m

?????文件???????2884??2018-01-04?21:47??信號處理代碼\mfiles\ambifunb.m

?????文件???????5455??2018-01-04?21:47??信號處理代碼\mfiles\ambifuwb.m

?????文件???????1620??2018-01-04?21:47??信號處理代碼\mfiles\amexpo1s.m

?????文件???????1628??2018-01-04?21:47??信號處理代碼\mfiles\amexpo2s.m

?????文件???????1654??2018-01-04?21:47??信號處理代碼\mfiles\amgauss.m

?????文件???????1601??2018-01-04?21:47??信號處理代碼\mfiles\amrect.m

?????文件???????1634??2018-01-04?21:47??信號處理代碼\mfiles\amtriang.m

?????文件???????1971??2018-01-04?21:47??信號處理代碼\mfiles\anaask.m

?????文件???????2060??2018-01-04?21:47??信號處理代碼\mfiles\anabpsk.m

?????文件???????2100??2018-01-04?21:48??信號處理代碼\mfiles\anafsk.m

?????文件???????1397??2018-01-04?21:48??信號處理代碼\mfiles\anapulse.m

?????文件???????2217??2018-01-04?21:48??信號處理代碼\mfiles\anaqpsk.m

?????文件???????1878??2018-01-04?21:48??信號處理代碼\mfiles\anasing.m

?????文件???????1390??2018-01-04?21:48??信號處理代碼\mfiles\anastep.m

?????文件???????5976??2018-01-04?21:48??信號處理代碼\mfiles\atoms.m

?????文件???????8126??2018-01-04?21:48??信號處理代碼\mfiles\Contents.m

?????文件???????3805??2018-01-04?21:48??信號處理代碼\mfiles\contwtgn.m

?????文件???????2688??2018-01-04?21:48??信號處理代碼\mfiles\contwtgnmir.m

?????文件???????2242??2018-01-04?21:48??信號處理代碼\mfiles\correlmx.m

?????文件???????5573??2018-01-04?21:55??信號處理代碼\mfiles\CVS\Entries

?????文件?????????12??2018-01-04?21:55??信號處理代碼\mfiles\CVS\Repository

?????文件?????????57??2018-01-04?21:55??信號處理代碼\mfiles\CVS\Root

?????文件???????1478??2018-01-04?21:48??信號處理代碼\mfiles\d2statio.m

?????文件???????1392??2018-01-04?21:48??信號處理代碼\mfiles\disprog.m

?????文件???????1134??2018-01-04?21:48??信號處理代碼\mfiles\divider.m

............此處省略122個文件信息