【正文】 h(1) = alpha0+alpha1.*(rmean0pheiomiga... .*rvar0)^2+beta1.*rvar0+gamma0.*Ht(1)。 for k2 = 2:T2 h(k2) =alpha0+alpha1.*(r(k21)phei... omiga.*h(k21)).^2+beta1.*h(k21)+gamma0.*Ht(k2)。 end impact2 =zeros(T2,1)。 impact2(1)=r(1)pheiomiga.*(alpha0+alpha1... .*(rmean0pheiomiga.* rvar0))+beta1.*rvar0+gamma0.*Ht(1)。 for k3 = 2:T2 impact2(k3)=r(k3)pheiomiga.*(alpha0+... alpha1.*(r(k31)pheiomiga.*h(k31)))+... beta1.*h(k31)+gamma0.*Ht(k3)。 end logml=T2.*log(2)T2.*log(gamma(lamda+1))ones(1,T2)... *log(sqrt(h.*gamma(lamda)./gamma(3.*lamda)))... .*ones(1,T2)*abs(impact2.*sqrt(4.^lamda... .*gamma(3.*lamda)./h.*gamma(lamda))).^(1/lamda)。 ML2(K2)=logml。 if ML2(K2) mx mx = ML2(K2)。 vector2 = [alpha0,alpha1,beta1,gamma0,phei,omiga,lamda]。 end end vectorl(X2,:)=vector2。endvectorlb=mean(vectorl)39。Para2=vectorlb。%r0=[rmean0。r([1:(T21)],1)]。ht=zeros(T2,1)。ht(1)= vectorlb(1)+vectorlb(2).*(rmean0vectorlb(5)... vectorlb(6).*rvar0)^2+vectorlb(3).*rvar0+vectorlb(4).*Ht(1)。for Y2=2:T2 ht(Y2)= vectorlb(1)+vectorlb(2).*(r0(Y21)vectorlb(5)... vectorlb(6).*ht(Y21))^2+vectorlb(3).*ht(Y21)+... vectorlb(4).*Ht(Y2)。end%et=rvectorlb(5)vectorlb(6).*ht。oht=[Ht,ht]。oet=[Et,et]。aic1=log(et39。*et/T2)+2*6./T2。aic2=log(Et39。*Et/T1)+2*4./T1。sc1=log(et39。*et/T2)+6*log(T2)/T2。sc2=log(Et39。*Et/T1)+4*log(T1)/T1。aic=[aic1,aic2]。sc=[sc1,sc2]。附錄B：function output=MCsimulate(fd)%% MCsimulate % Simulate the future station of equation group1 and equation group1 with% Monte Carlo this praogram also draw the graphs of output% FD is the days that be predicted% OUTPUT is the output of simulation% format long sHt0=。 sEt0=。 sHt=zeros(fd,1)。sHt(1)=+*sEt0^2+*sHt0。 sEt=zeros(fd,1)。sEt(1)=normrnd(0,1)*sqrt(sHt(1))。 for K1=2:fd sHt(K1)=+*sEt(K11)^2+*sHt(K11)。 sEt(K1)=normrnd(0,1)*sqrt(sHt(K1))。 end Rt=+sEt。 slht0=。slet0=。slht=zeros(fd,1)。 slht(1)=+*slet0^2+*slht0+*sHt(1)。 miut=zeros(fd,1)。 miut(1)=+*slht(1)。 slet=zeros(fd,1)。slet(1)=normrnd(0,1)*sqrt(slht(1))。 rt=zeros(fd,1)。 rt(1)=miut(1)+slet(1)。 for K2=2:fd slht(K2)=+*slet(K21)^2+*... slht(K21)+*sHt(K2)。 miut(K2)=+*slht(K2)。 slet(K2)=normrnd(0,1)*sqrt(slht(K2))。 rt(K2)=miut(K2)+slet(K2)。 end output=[sHt,slht,Rt,rt]。% figure subplot(2,1,1),plot(sHt),axis([0 fd 0 7e5]) title(39。未來上證180條件方差預測39。) subplot(2,1,2),plot(slht),axis([0 fd 0 2e4]) title(39。未來中國銀行條件方差預測39。) figure subplot(2,1,1),plot(Rt),axis([0 fd ]) title(39。未來上證180收益預測39。) subplot(2,1,2),plot(rt),axis([0 fd ]) title(39。未來中國銀行收益預測39。)致 謝當海風與烈日帶著又一個夏季來臨的訊號時，我的大學生活也即將結(jié)束?；仡欉@段說長不長說短不短的時光，我思緒萬千。想著它即將離我遠去，我依稀有幾分不舍。慶幸，我?guī)缀蹩梢詰浵肫鹈刻熳鲞^些什么事，學的了哪些知識，我因自己的勤奮而倍感欣慰，我也因我學有所成而倍感充實。我最想感謝的是海南大學的圖書館，在此，我懷著最虔誠的真心，對陪我走過四年的海大圖書館致以最崇高的謝意，感謝它給我提供諸多優(yōu)秀的古今中外書籍，使我在大學里飽嘗知識的滋潤，并獲益匪淺。再要感謝我的論文指導老師徐艷教授，謝謝徐教授在這篇文章的撰寫過程中給予的幫助、支持與修改意見。歲月如風一般飄走，不因人的懷舊而停留。在經(jīng)歷了大學生活洗禮后，我將滿懷希望地投入到現(xiàn)實社會的考驗之中。