求解CPPIStr函数代码（固定比例组合保险策略模拟）

问题：固定比例组合保险策略模拟（问题在代码括号内）
CPPIStr函数源代码：

1. function [F,E,A,G,SumTradeFee,portFreez]=CPPIStr(PortValue,Riskmulti,GuarantRatio,TradeDayTimeLong,...
   
2.     TradeDayOfYear,adjustCycle,RisklessReturn,TradeFee,SData)
   
3. %input
   
4. %PortValue,Riskmulti,GuarantRatio,TradeDayTimeLong,TradeDayOfYear,adjustCycle,Risklessreturn,TradeFee...
   
5. %SData is simulation index data
   
6. %output
   
7. %F,E,A,G,SumTradeFee
   
8. %portFreez default is 0, if portFreez=1,portfolio freez there would
   
9. %have no risk--investment
   
10. %%
    
11. %初始交易费用（交易佣金）为0；
    
12. SumTradeFee=0;
    
13. %F,E,A,G的初始化，长度为N+1            ？？（ 长度为什么要N+1）？？
    
14. F=zeros(1,TradeDayTimeLong+1);
    
15. E=zeros(1,TradeDayTimeLong+1);
    
16. A=zeros(1,TradeDayTimeLong+1);
    
17. G=zeros(1,TradeDayTimeLong+1);
    
18. %给定F,E,A,G的初始值
    
19. %初始组合资产
    
20. A(1)=PortValue;
    
21. %初始安全底线
    
22. F(1)=GuarantRatio*PortValue*exp(-RisklessReturn*TradeDayTimeLong/TradeDayOfYear);
    
23. %初始风险资产
    
24. E(1)=max(0,Riskmulti*(A(1)-F(1)));
    
25. %无风险资产
    
26. G(1)=A(1)-E(1);
    
27. %%
    
28. %是否进行风险资产平仓
    
29. %portFreez=0 正常,portFreez=1;平仓
    
30. portFreez=0;
    
31. %if portFreez=1,portfolio freez there would have no risk--investment
    
32. %%
    
33. %开始逐日模拟，循环计算
    
34. %根据T-1日情况与T日市场行情，计算T日产品净值
    
35. for i=2:TradeDayTimeLong+1
    
36.     E(i)=E(i-1)*(1+(SData(i)-SData(i-1))/(1+SData(i-1)));             ？ ？为什么这么算？？
    
37.     G(i)=G(i-1)*(1+RisklessReturn/TradeDayOfYear);
    
38.     A(i)=E(i)+G(i);
    
39.     F(i)=GuarantRatio*PortValue*exp(-RisklessReturn*(TradeDayTimeLong-i+1)/TradeDayOfYear);     
    
40.     %判断是否进行调仓，调仓周期为adjustCycle
    
41.     %mod函数求余数的意思，若adjustCycle=20, i为20的整数倍时
    
42.     %mod(i,adjustCycle)=0
    
43.     if mod(i,adjustCycle)==0
    
44.         temp=E(i);                  ？这个temp放这里干嘛用？                        E(i)=max(0,Riskmulti*(A(i)-F(i)));
    
45.         SumTradeFee=SumTradeFee+TradeFee*abs(E(i)-temp);               ？  ？abs里不是为0了吗？？
    
46.        G(i)=A(i)-E(i)-TradeFee*abs(E(i)-temp);
    
47.     end
    
48.     %判断是否平仓，若风险资产为0，组合平仓冻结
    
49.     if E(i)==0
    
50.         A(i)=G(i);
    
51.         portFreez=1;
    
52.     end
    
53. end

复制代码

函数测试代码：（？？最后产品的收益率在哪里看到的？？）

1. %set value
   
2. PortValue=100;%产品组合初始价值
   
3. Riskmulti=2;%产品风险乘数为2
   
4. GuarantRatio=1;%产品保本率为100%
   
5. TradeDayTimeLong=250;%产品期限为250个交易日
   
6. TradeDayOfYear=250;%模拟假设一年交易日为250个
   
7. adjustCycle=10;%调整周期为每10个交易日调整一次
   
8. RisklessReturn=0.05;%无风险产品收益率为5%
   
9. TradeFee=0.005;%风险资产的交易费用为0.5%
   
10. 
    
11. %to generate Brow random number
    
12. %预期收益率年化转日化
    
13. Mean=1.2^(1/TradeDayOfYear)-1;
    
14. %预期波动率年化转日化
    
15. Std=0.3/sqrt(TradeDayOfYear);
    
16. %初始价格
    
17. Price0=100;
    
18. SData=RandnPrice(Price0,Mean,Std,TradeDayOfYear);
    
19. %将初始价格并入随机价格序列
    
20. SData=[Price0;SData];
    
21. %[X0,X1,~,Xn]
    
22. 
    
23. %调用CPPIStr函数
    
24. [F,E,A,G,SumTradeFee,portFreez]=CPPIStr(PortValue,Riskmulti,GuarantRatio,TradeDayTimeLong,...
    
25.     TradeDayOfYear,adjustCycle,RisklessReturn,TradeFee,SData);
    
26. %to plot
    
27. figure;
    
28. %子图1，模拟的风险资产的价格序列
    
29. subplot(2,1,1)
    
30. plot(SData)
    
31. legend('Hs300-Simulation')
    
32. xlabel('t');
    
33. ylabel('price');
    
34. %set(gca,'XTick',10:10:100) 设定刻度
    
35. %set(gca,'YTick',10:10:100)
    
36. %CPPI策略的运行情况
    
37. subplot(2,1,2)
    
38. plot(A,'-.')
    
39. hold on
    
40. plot(E,'-o')
    
41. plot(F,'-k')
    
42. plot(G,'-x')
    
43. %标记线型
    
44. legend('PortValue','RiskAssect','GuarantLine','RisklessAssect')
    
45. xlabel('t');
    
46. ylabel('price')
    
47. axis([0 300 0 150])
    
48. %总的交易费用
    
49. SumTradeFee

复制代码