【正文】 function a=xtxs(x,p)n=length(x)。m=0。X=zeros(np1,p)。for i=1:np1。 for j=1:p。 X(i,j)=x(j+m)。 end m=m+1。endi=np。for j=1:p1。 X(i,j)=x(m+j)。endx(i,p)=x(m+p)。x1=zeros(np)。for i=1:np。 x1(i)=x(i+p)。endi=np。x1(i)=x(i+p)。a=(X\x139。)。%xk=x*a%k=xk39。（五）.數(shù)據(jù)載入子程序function T=gett(ha,tv)%get the string or value of haif tv==1 T=str2double(get(ha,39。string39。))。elseif tv==2 T=get(ha,39。value39。)。elseif tv==3 T=get(ha,39。string39。)。End最后通過鍵入道瓊斯工業(yè)指數(shù)的信源子程序，進(jìn)行預(yù)測。如圖42圖42點(diǎn)擊“指導(dǎo)老師”或者“設(shè)計(jì)人員”菜單項(xiàng)都可彈出相應(yīng)的姓名。如圖44圖44圖42預(yù)測的是用10階的預(yù)測器預(yù)測未來10周的股票最高價(jià)趨勢。在預(yù)測類型的下拉菜單中可以選擇“最高價(jià)預(yù)測”、“最低價(jià)預(yù)測”、“開盤預(yù)測”或“收盤預(yù)測”四個(gè)選項(xiàng)如圖45，圖45同時(shí)在預(yù)測階數(shù)和預(yù)測周數(shù)文本框中填寫數(shù)據(jù)得到不同的預(yù)測值。本設(shè)計(jì)模型的信源長度為104，隨著階數(shù)的增大，預(yù)測值精度不斷提高，但當(dāng)階數(shù)增大到26時(shí)，其誤差到達(dá)最小，如果繼續(xù)增大階數(shù)，已經(jīng)不再影響其精度。驗(yàn)證方法為：在道瓊斯指數(shù)的四個(gè)信源中取一種信源，其長度為104，截取最后的20個(gè)信源，把前84個(gè)信源載入預(yù)測器中。在周數(shù)相同的情況下，選擇不同的階數(shù)來預(yù)測后20個(gè)信源，把預(yù)測出來的結(jié)果與截取的后20個(gè)信源相比較，再求其誤差就可觀察該模型的精確度。結(jié)論本設(shè)計(jì)內(nèi)容在理解信號(hào)與系統(tǒng)基本原理的前提下，利用MATLAB設(shè)計(jì)了一個(gè)線性預(yù)測系統(tǒng)，該系統(tǒng)利用一個(gè)離散時(shí)間有限脈沖響應(yīng)（FIR）濾波器來解決屬于預(yù)測建模等問題。這是一個(gè)基于MATLAB計(jì)算機(jī)仿真的股票線性預(yù)測模型，它用股票的開盤、收盤、最高、最低四種價(jià)位為源信號(hào)進(jìn)行預(yù)測，期望用選擇濾波器的階數(shù)來調(diào)整它的精確度，做到預(yù)測誤差最小。參考文獻(xiàn)① 作者：張威，《MATLAB基礎(chǔ)與編程入門》，西安電子大學(xué)出版社出版社，2007年第1版，第4頁。② 作者：黃忠霖，《控制系統(tǒng)MATLAB計(jì)算及仿真》，國防工業(yè)出版社,2007年底版，第134頁。③ 作者：吳兆熊，《數(shù)字信號(hào)處理》，西安交通大學(xué)出版社，2006年第1版，第82頁。④ 作者：王嘉梅，《基于MATLAB的數(shù)字信號(hào)處理與實(shí)踐開發(fā)》，西安電子大學(xué)出版社，2007年弟1版，第157頁。附錄High,Low,Open,Close21原始數(shù)據(jù)function zz=High()High=[ ... ... ... ... ... ... ... ... ... ... ... ... ... ... ]。zz=High39。function zz=Low()Low=[ ... ... ... ... ... ... ... ... ... ... ... ... ... ... ]。zz=Low。function zz=Close21()Close21=[ ... ... ... ... ... ... ... ... ... ... ... ... ... ... ]。zz=Close2139。function zz=Open()Open=[ ... ... ... ... ... ... ... ... ... ... ... ... ... ... ]。zz=Open39。