风险管理hansenmodel作业.docVIP

风险管理作业■■■■Hansens Basic Model 1、给定校准值： 0 = 0.99 3 = 0.025 0 = 0.36 利用稳态方程组，解出稳态值，并且校准参数A。对比书上结果（P108） 欧拉方程： c Y 1=0 匐才(心+1-6)] AC =(1 —0)(1 —的亡 Cznl 约束条件： y 严 陰=(1-5)K + Z-G t=a 卞 随机过程：+严0人+￡阿 Y 稳态方程组:f l=0?(a肓+ 1-5) 稳态方程组: K Y = Ka^~a 7 8K=Y-C 校准：力=0.025 0 = 0.99 a = 0.36 People spent a third of their time working, so H=l/3. MATLAB求解方程组： alpha=0.36;beta=0.99;delta=0.025;H=l/3; syms AC K Y; S1 =beta*(alpha*(Y/K)+1 -delta)-1; S2=(l-H)*(l -alpha)* Y?A*C*H; S3=(KAalpha)*(HA(l -alpha))-Y; S4=Y-C-delta*K; [C,A,K,Y]=vpasolve(S 1 ,S2,S3,S4,C, A,K, Y) C =0.91810915771217469189791046671099 A=1.7213622291021671826625386996904 K = 12.663084512717455889799353341169 =1.2346862705301110891428943002402 校准后得到A=1.72,再把A带回方程组校准H。 alpha=0.36;beta=0.99;delta=0.025;A=1.72; syms HCK Y; S1 =beta* (alpha* (Y/K) +1 -delta)-1; S2=(l -H)*(l -alpha)*Y^A*C*H; S3=(KAalpha)*(HA( 1 -alpha))-Y; S4=Y-C-delta*K; [C,H,K,Y]=vpasolve(Sl,S2,S3,S4,C,H,K,Y) C =0.91859378752813681894435247137647 H =0.33350928547435097653438264791977 K = 12.669768803213156434665914117592 = 1.2353380076084657298110003243163P108 R=vpa(alpha* Y/K) R =0.03510101010101010101010101010101 H二0.3335,符合实际情况，运算结果和书上一致。 2、给定Q = 0.95,利用公式计算均衡运动方程组，对比书上结果(P109)。 均衡运动方程组： 孔=% ? K+% ? A ^t=Vrk ?斤/+匕2 ?2 简化后的均衡运动方程组: xt = Pxt_{+Qzt X =Z_i+S?z, 其中, xt = Kf+] Xt-\ ~ ?)1 =妆R =Vck%S =丿?7 Xt-\ ~ ? )1 = 妆 R = Vck % S = 丿 ? 7 方程组系数为: p = % 代入求解过程: Q = VkA A = [0 -K o or, B = [0 (1 — 5)K a o -ir, F = G = H=O 丿= (0, —1,0,07) K =(o, 1, 0, 0) L = M=0 N =p 求解系数： 0 = (F-JC]A)P2-(JC-}B-G + KC]A)P-KC}B + H 求根公式tP，0Pl r = —C?(AP + B) [N(F- JC-A) + (JR + FP + G- KC~}A)]- Q _ q =(Q D — L)N + KL D-M S = -C[\AQ + D) MATLAB 求解： A=[0 -K 0 or； B=[0;((l -delta)*K);alpha;-1 ]; C=[l -1 (-l/(l-H))0;Y-C0 0;-l O(l-alpha) 0;l 00-1]; D=[0 0 1 0]1; F=[0]; G=F; H=F; J=[0-l 0 (beta*R)]; K=[0 10 0]; L=F;M=F; N 二[0.95]; syms PQ R S; S5=(F-J*(CA(-1 ))*A)*(PA2)-(J*(CA(-1 ))*B-G+K*(CA(-1 ))*A)*P-K*(CA(-1 ))*B+H; P=vpasolve(S5,P) P = 0.95367388822440608100088193742432 1.0591681523