QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 3736|回复: 0
打印 上一主题 下一主题

[建模教程] 投资组合问题

[复制链接]
字体大小: 正常 放大
浅夏110 实名认证       

542

主题

15

听众

1万

积分

  • TA的每日心情
    开心
    2020-11-14 17:15
  • 签到天数: 74 天

    [LV.6]常住居民II

    邮箱绑定达人

    群组2019美赛冲刺课程

    群组站长地区赛培训

    群组2019考研数学 桃子老师

    群组2018教师培训(呼伦贝

    群组2019考研数学 站长系列

    跳转到指定楼层
    1#
    发表于 2020-6-16 10:14 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta |邮箱已经成功绑定
    基本的投资组合模型5 q0 j1 F( ^' \- T! n
    例5  美国某三种股票( A ,B , C)12年(1943-1954)的价格(已经包括了分红在内)每年的增长情况如表6所示(表中还给出了相应年份的500种股票的价格指数 的增长情况)。例如,表中第一个数据1.300的含义是股票 A在1943年的年末价值是 其年初价值的1.300倍,即收益为30%,其余数据的含义依此类推。假设你在1955年 时有一笔资金准备投资这三种股票,并期望年收益率至少达到15%,那么你应当如何 投资?当期望的年收益率变化时,投资组合和相应的风险如何变化?
    - }: R% J; u& W4 [- _9 i1 J
    7 t0 {8 R( {( }  y" ^9 s
    & J& y6 p7 F4 g2 t$ i5 Q" W8 K% l+ `- s$ @; x4 e

    9 z& f( R% `2 c: C1 Q0 y6 f9 I( J" D% \, h/ b
    (1)问题分析9 U$ e# E0 u/ T3 X# |- l- W. J% @

    9 c( I! J2 A( [: J本例的问题称为投资组合(portfolio)问题,早在1952年Markowitz就给出了这个模型的基本框架,而且这个模型后来又得到了不断的研究和改进。一般来说,人 们投资股票时的收益是不确定的,因此是一个随机变量,所以除了考虑收益的期望值外, 还应当考虑风险。风险用什么衡量?Markowitz建议,风险可以用收益的方差(或标 准差)来进行衡量:方差越大,则认为风险越大;方差越小,则认为风险越小。在一定 的假设下,用收益的方差(或标准差)来衡量风险确实是合适的。为此,我们先对表6 中给出的数据计算出三种股票收益的均值和方差(包括协方差)备用。
    / F0 D! w1 D3 e( G8 N; q* w+ k; w9 f% w  P
    一种股票收益的均值衡量的是这种股票的平均收益状况,而收益的方差衡量的是这 种股票收益的波动幅度,方差越大则波动越大(收益越不稳定)。两种股票收益的协方 差表示的则是它们之间的相关程度:
    6 }8 e& N% R; ?
    : y) `* C0 ^+ J  l$ g- @i)协方差为0时两者不相关。
    ( a- l/ ?5 L7 _) g: t' o7 L0 |- K6 b
    ii)协方差为正数表示两者正相关,协方差越大则正相关性越强(越有可能一赚 皆赚,一赔俱赔)。( R  Q, m9 a* a/ Y  o% C% p

    ) e7 `6 T( g3 L# K$ E5 R( Ciii)协方差为负数表示两者负相关,绝对值越大则负相关性越强(越有可能一个 赚,另一个赔)。 + f/ i" g* L: u: ^5 }1 f: f
    2 ]  R8 h5 E3 @* n) L& y- J1 o
    记股票 A , B,C  每年的收益率分别为 (注意表中的数据减去1以后才是年收益率),则    是一个随机变量。用E 和D分别表示随机变量的数学期望和方差算子,用cov表示两个随机变量的协方差(covariance),根据概率论的知 识和表6给出的数据,则可以计算出年收益率的数学期望为
    1 h  v+ M: c$ ]; u! \" J  m
    ; P/ K6 Y" d" C0 d# j; i = 0.0890833, = 0.2136667, = 0.2345833        (1)
    - u, R8 t8 e& F- C8 e+ h* U2 C3 d& N1 y; A; _$ @# a: ?- D, p
    同样,可以计算股票 A ,B ,C年收益率的协方差矩阵为 % ~' N# X1 n5 C/ ]; y
    , |: ]0 B- O+ V( N) e0 T+ K# c/ z
    $ y/ [8 Y2 E2 }  H' z
    ) H+ t4 l  R" {- g" l0 i* ?, l
    计算的LINGO程序如下:  - l; M. I5 \. M8 R' I
    . ~# w+ L4 f* d- A2 E2 Y
    MODEL: 2 W6 K8 k6 G7 A, S* y4 U
    Title 均值向量Mean与协方差矩阵COV; . ^# l' p8 C! o2 U- c1 Q
    SETS:    ) M, o9 v; J4 P0 x5 x, o) D0 d# \2 H
        YEAR/1..12/;   
    2 [! k# d3 Q- [% n, |# b" s4 u    STOCKS/ A,  B,  C/: Mean;   
    - A0 f1 L9 f" z* i0 k# ~. K! |* W! G" q    link(YEAR, STOCKS): R;   
    9 S% ?, w. w( S% ^0 O    temp/1..5/;       9 f6 [8 v) h+ B" B
        tmatrix(YEAR,temp):tm;   
    " v3 ?0 T6 @* S, M1 C* Z    STST(Stocks,stocks): COV;
    # c" X5 c% T/ f+ h  F4 \ENDSETS 3 X7 ?2 \) L# V. j, B
    DATA:    4 t2 I+ a) c( y# v! h7 f) ]" f( E. P
    tm = + @! Z0 ~' @9 k+ R6 ^/ J5 {
    1943           1.300          1.225           1.149           1.258997 : f+ M( K% v) [/ C
    1944           1.103          1.290           1.260           1.197526
    & F9 }; _3 c6 v% E8 C0 u1945           1.216          1.216           1.419           1.364361
    " l4 e. Z; z( |6 X% ^1946           0.954          0.728           0.922           0.919287
    ' A/ K+ a1 N: m1947           0.929          1.144           1.169           1.057080 ! M% b$ \0 \8 ]$ C
    1948           1.056          1.107           0.965           1.055012
    6 h5 w7 S. D0 r* ~% {% K/ @1949           1.038          1.321           1.133           1.187925 ! v: d0 z9 y# k! _5 q& a& X
    1950           1.089          1.305           1.732           1.317130 / Y; {+ ?1 ]5 ]! j, v2 a3 q1 ~
    1951           1.090          1.195           1.021           1.240164
    # D- ]8 i! B& d. D; W$ A  _1952           1.083          1.390           1.131           1.183675 % [0 l; e* ]0 A) ~
    1953           1.035          0.928           1.006           0.990108   
    8 i3 ?$ D, r9 S/ \. f9 r 1954           1.176          1.715           1.908           1.526236;
    . k! V# R; c! K* P: y6 l" k@text(data1.txt)=R; 8 x' {3 ~# a9 X4 V
    @text(data2.txt)=Mean;
    2 g# d5 k* d1 U( u@text(data3.txt)=COV; , D% L% v1 ?9 x( N6 \
    ENDDATA " a* ?! u, ?# t( q4 R  m" A; n0 T
    CALC: !计算均值向量Mean与协方差矩阵COV; % W' c* Q4 }; ~4 k1 O
    @for(tmatrix(i,j)|j #ge#2 #and# j #le# 4:R(i,j-1)=tm(i,j)-1); 2 ^6 S& R9 ^3 R
    @for(stocks(i): Mean(i) = @sum(year(j): R(j,i)) / @size(year) );  
    " n! Z9 e2 M4 {: \6 @+ L6 [: B- [@for(stst(i,j): COV(i,j) = @sum(year(k): (R(k,i)-mean(i))*(R(k,j)-mean(j))) / (@size(year)-1) );
    8 x6 O3 n/ p& g0 w1 g8 MENDCALC
    7 d' H5 U8 E$ f5 ?" \& ]END
    : o; Q% V0 }( p8 n: i+ `% E      注意模型中计算协方差矩阵COV时,分母是样本数减去1,而不是样本数,这是常 用的计算方法,主要是为了保持这个估计的无偏性(当然,样本数较大时两者差别不太)。4 @, A) ?4 I% g5 I7 t  p
    (2)模型建立
    , C, w- P! t) L+ @  r# \, |' _         用决策变量  分别表示投资人投资股票 A ,B , C 的比例。假设市场上没有其它投资渠道,且手上资金(可以不妨假设只有1个单位的资金)必须全部用于投资这 三种股票,则
    9 p2 m' ^" \! |$ K7 D2 ^9 o1 \  i, n/ x- r+ Q% s  D" y
                 ( 2  )
    4 i$ v% I+ w# V! s& c: h- C# d* c8 ?% m% ^/ L; p
    年投资收益率   也是一个随机变量。根据概率论的知识,投& l5 r" G6 d+ S9 R
    资的总期望收益为          ( 3 )
    # j. G7 h" \# Z5 H$ b  Z& o. g/ I5 B' i/ }1 u& E$ p
       ( 4 )
    5 e, [. ]5 N2 f0 X2 h' \, H6 h+ W  W" j7 b: E) e: J6 G* n
    实际的投资者可能面临许多约束条件,这里只考虑题中要求的年收益率(的数学期 望)不低于15%,即 & E  r! \" u- J. w1 q! g

    ' X$ l# B- c' n: z, S             ( 5 )2 U0 [% q! S8 V" o8 q- y: P
    3 y6 r) T. J' c9 L; w9 r2 S0 n
    所以,最后的优化模型就是在约束(2)和(5)下极小化(4)。由于目标函数V 是决策变量的二次函数,而约束都是线性函数,所以这是一个二次规划问题。
    9 W/ O( m  y6 ~
    9 ^, H1 X- o3 N& @# o1 X(3)用LINGO求解模型 # D4 d: N5 |6 h0 d# A
    $ F( N$ b$ Q4 \
    编写LINGO程序如下
    * ?. w8 }4 E* X& V' b0 `" _  l. p( i. q; h# ~+ b% I1 L
    MODEL:
    ; B: e6 c6 n2 s/ z. E2 A* ~SETS:   
    ( j+ h, ^" K' Y1 u$ V) _    YEAR/1..12/;   
    + o, p/ r) b" E5 l, t  g8 m! x    STOCKS/ A,  B,  C/: Mean,X;    2 G+ p3 Y; }. z+ b
        link(YEAR, STOCKS): R;    * ~+ y6 n# u8 V$ M( }1 x
        STST(Stocks,stocks): COV;
    $ K2 o* g% a9 N6 N: vENDSETS
    - E" a. v! {9 m* H! zDATA:
    $ V, j* I& o) l3 B9 f% N5 E" U1 hTARGET=0.15; $ t4 W2 ^3 }. w" B) K; \
    Mean=@file(data2.txt);
    / e7 {. c& L9 N: b" k& eCOV=@file(data3.txt); 6 t& {% g2 L  }' W$ V3 o4 f2 e! ^1 j
    ENDDATA . x3 O. q0 o- Q9 @& f- {
    [OBJ] MIN = @sum(STST(i,j): COV(i,j)*x(i)*x(j));
    % E( |5 m; a, h* _[ONE] @SUM(STOCKS: X) = 1; - L8 f" I! _9 Y' D( {
    [TWO] @SUM(stocks: mean*x) >= TARGET;
    4 K, u6 ~9 Z& U: H. g- j! \END# I! m( Z, B- k* \' }: g9 i% w
    求得投资三种股票的比例大致是: A占53%,B 占36%,C 占11%。风险(方 差)为0.02241378。; P* L7 ~4 y/ B! b5 f* z
    8 H5 D. @- K1 y" o) \
    (4)用MATLAB软件对模型进行参数分析
    6 f% }' |- J0 C( s; F3 b/ M# }* n+ \0 H# A9 u6 P: e4 J
    对实际投资人来说,可能不仅希望知道指定的期望投资回报率下的风险(回报率的方差),可能更希望知道风险随着不同的投资回报率是如何变化的,然后作出后的投 资决策。这当然可以通过在上面的模型中不断修改约束中的参数(目前为0.15)来实 现,如将0.15改为0.2345,则表示投资回报率希望达到23.45%(这几乎是可能达到的最大值了,因为这几乎是三种股票中最大的投资回报率,即股票C 的回报率)。可 以想到,这时应主要投资在股票C 上。实际求解一下,可以知道最优解中投资股票C 的 份额大约是99.6%(剩余的大约0.4%投资在股票B 上)。 目前LINGO软件还没有二次规划灵敏度分析的功能。下面我们利用MATLAB软件进 行灵敏度分析,回报率的取值区间为[0.09,0.234],变化步长为0.002。
    & J: v4 ^3 R" B$ o/ O! `9 R9 B$ U
    1 A; |! S. M. ?  Z编写程序 如下:   L8 }+ m+ R* y; v7 D! w
    # Q/ x" F' I" W% O" J2 i7 @
    clc,clear & X# b' R% m3 v
    load data2.txt,load data3.txt
    % V9 g6 s; l* l/ c4 q; hh=reshape(data3,[3,3]); ) W- Z/ V$ e' Q8 o/ t5 D( V
    a=data2'; : r/ I* v2 V1 y+ C
    solution=[]; ! }! W: x  j* e$ I3 e# ?
    target=0.09;
    6 n0 a6 n: S5 V% f: Zhold on
    9 `/ |/ d: c  H8 O3 n7 D8 Gwhile target<0.234     
    ' }9 B' n# I# @7 U6 h3 I/ [6 G7 j+ v    [x,y]=quadprog(2*h,[],-a,-target,ones(1,3),1,zeros(3,1));  
    ; ]6 X, y+ u- H$ Z) {. P' A    plot(target,y,'*b');     # Y& m; F7 `& C2 [- A4 }( N, K3 o
        solution=[solution [target;x;y]];     
    5 X! T* U% z4 W0 ?: f% c    target=target+0.002;     
    3 m$ N5 [1 N1 q    end 3 R7 r1 G% R$ g
    solution & `5 V8 o6 U3 K
             得到的投资回报率与风险之间的关系曲线如图3所示 。
    ' P! P5 c, ]8 |8 S' ~& P& v7 t1 _4 N/ O- w/ u
    3 O6 A' Q/ G( r- n3 ~) I6 s

    0 `: E8 h, F( u' x从图3可以看出,投资回报率从22%附近开始,风险迅速增大。 2 O9 t; R, V9 ], F9 L

    / x$ ~8 d6 S7 z
    ' v% y! q7 V8 P( l' s- I/ X! Q) f$ f: L% W% M" M
    存在无风险资产时的投资组合模型 6 Z6 b: p( u; t9 v/ l
    例6  假设除了例5中的三种股票外,投资人还有一种无风险的投资方式,如购买国库券。假设国 库券的年收益率为5%,如何考虑例5中的问题?* Y: z$ a& g/ \! b+ J2 Z

    % F/ J9 G- ^8 v- w  B, a(1)问题分析
      f* t% y7 p/ k! _
    % H: P" Z& r: M5 x其实,无风险的投资方式(如国库券、银行存款等)是有风险的投资方式(如股票)的一种特例, 所以这就意味着例5中的模型仍然是适用的。只不过无风险的投资方式的收益是固定的,所以方差(包 括它与其它投资方式的收益的协方差)都是0。
    9 s- {2 _% n" j, t8 H* t! R- c9 n4 X6 ]
    " k9 q" A  ?  t& C/ q(2)问题求解
    ) W$ s' j( ?! v/ }" v. U+ N
    * F5 y* H& C0 L0 T假设国库券的投资方式记为D,则当希望回报率为15%时,对应的LINGO模型如下: $ Q6 ~& W6 Y* J1 Y7 U& i- d/ {& |" V, d

    ) v9 i& @- H1 K) W- m/ NMODEL: . K! W% k1 T# o6 y9 O
    Title 含有国库券的投资组合模型;
    0 h/ t$ ]- U$ y! h0 RSETS:    8 _, b+ G2 |6 C7 X
        STOCKS1/ A,  B,  C/: mean1;    # j1 f; Z4 S4 o  `' P
        STST1(Stocks1,stocks1): COV1;   
    & D4 o8 Y1 B* ?6 f: |4 v    STOCKS/ A,  B,  C, D/: mean,X;    $ H3 z: M5 R6 K8 x8 \6 q
        STST(Stocks,stocks): COV;
    # ^" F% S4 P/ y5 n3 KENDSETS # P' N0 K1 W7 m" D
    DATA:
    ) a7 B. X5 P# s6 t( }" \& c4 W; Zmean=0.05; 1 Z# t9 s! p: I9 S1 J5 T5 P; S
    COV=0; 9 `/ D( Z: e; I2 T% v
    COV1=@file(data3.txt);
    # |2 U0 N, o- BMean1=@file(data2.txt);   $ o2 d/ g8 Z0 ^( N8 N' X
    TARGET = 0.15; ! 0.10;
    6 P! p% A: a* B4 U! I' aENDDATA
    , _5 V9 L: B4 \/ d; I. R7 hcalc: ) V' a& \  T+ n7 m, k, G3 y
    @for(STOCKS(i)|i #ne# 4: mean(i)=mean1(i)); ( [8 b- l+ u& f. r% r: M
    @for(STST(i,j)| i #ne# 4 #and# j #ne# 4:COV(i,j)=COV1(i,j));  ) q0 c2 ^7 [9 L4 e
    endcalc5 Q8 _, j# l) i0 a6 A6 d" d3 d8 a
    [OBJ] MIN = @sum(STST(i,j): COV(i,j)*x(i)*x(j));
    6 T1 ]% r4 a. K5 t[ONE] @SUM(STOCKS: X) = 1;
    * Y. {2 E( v& K; {: S[TWO] @SUM(stocks: mean*x) >= TARGET;
    " I# g$ `# @! o. O% vEND
    0 ]% V! p) V  q# _" J) g. P计算结果为,投资 A占8.7%,B 占42.9%,C 占14.3%, D (国库券)占34.1%,风险 (方差)为0.02080347。与例5中的风险(方差为0.02241378)比较,无风险资产的存 在可以使得投资风险减少。虽然国库券的收益率只有5%,比希望得到的收益率15%小很多, 但在国库券上的投资要占到34.1%,其原因就是为了减少风险。
    " [$ s; {' }) D4 ?4 [4 I! J* C) g+ ?/ R/ t, ~
    现在,我们把上面模型中的期望收益减少到10%,即把数据段中的语句“TARGET = 0.15”改为“TARGET = 0.10”,重新求解模型。计算结果如下:% T: P* \* k" d# B7 i$ U$ ^
    . G5 e+ r( v* k
    投资 A占4.3%,B 占21.4%,C 占7.2%, D (国库券)占67.1%,此时风险(方差为 0.0052)进一步下降。请特别注意:你能发现这个结果(这里不妨称为结果2)与刚才 “TARGET = 0.15”的结果(这里不妨称为结果1)有什么联系吗? 仔细观察这两个结果,可以发现:结果2中投资在有风险资产(股票 A ,B , C)上的比例大约都是结果1中相应的比例的一半。也就是说,无论你的期望收益和风险偏好如何,你 手上所持有的风险资产本身相互之间的比例居然是不变的!变化的只是投资于风险资产与 无风险资产之间的比例。有趣的是,这一现象在一般情况下也是成立的,一般称为“分离 定理”,即风险资产之间的投资比例与期望收益和风险偏好无关。1981年诺贝尔经济学奖 得主Tobin教授之所以获奖,很大一部分原因就是因为他发现了这个重要的规律。 也正是由于有这样一个重要结果,我们在下面各节的讨论中就不再考虑存在无风险资产 的情形了,而只考虑确定风险资产之间的投资比例。
    , s: ?8 `# S* a4 O% k+ N- g, e1 M7 G: U6 C( u+ U3 s2 A+ t
    * @( [( E3 y- B
    & J, ?# w! L! o" D
    考虑交易成本的投资组合模型/ u, U: ~5 D7 ]' |- E% f; b
    例7  继续考虑例5(期望收益仍定为15%)。假设你目前持有的股票比例为:股票 A占 50%,B 占35%,C 占15%。这个比例与例5中得到的优解有所不同。实际股票市场上 每次股票买卖通常总有交易费,例如按交易额的1%收取交易费,这时你是否仍需要对所持 的股票进行买卖(换手),以便满足优解的要求?& _: @# s6 p! @& {9 Q

    / p0 b1 V# e6 p& s" f(1)建立模型 0 `3 x: c! E+ @2 }4 B
    + l$ F; y! o) k% f# K5 f
    仍用决策变量   分别表示投资人应当投资股票 A 、B 、C 的比例,进一步假设购买股票 A 、B 、C 的比例为 2 ,卖出股票 A 、B 、C 的比例为 。其中,   和 ( i=1,2,3 )中显然多只能有一个严格取正数,且 6 J0 w- l6 T: g& H$ J& A2 Z

    2 |. S% ^. t( |- s$ a3 n                      ( 6 )
    " H7 I5 R+ |& e$ O3 e
    9 w" Q* P4 r/ w' T0 z! ]2 r由于交易费用的存在,这时约束   不一定还成立(只有不进行股票买卖,即  时,这个约束才成立)。其实,这个关系式的本质是:当前持有的总资金是守恒的,在有交易成本(1%)的情况下,应当表示成如下形式:) W# k: @9 D( F0 `1 q( c" [9 x

    / M  a1 j- E1 C& V              (  7 ) ! X7 V( O0 c8 m; C

    $ x+ I  g4 T+ f另外,考虑到当前持有的各只股票的份额    之间也应该满足守恒关系式 6 ^1 X7 A1 U7 n) i1 d

    . h" d+ j4 a/ ^                     (  8 )  # F) q' J9 R- u2 v  }% {+ d% |2 f( }

    # I. d  S" O4 C1 ]; @" W(2)模型求解7 Q  X3 n' V/ A' W( @+ ]  c8 w
    ( C2 s1 D  `. J& t$ R- l7 H% s* W6 g
    问题对应的LINGO程序如下:  " s! Y$ X% S# L
    6 [$ |, W) H* g: M
    MODEL: : U+ X! C; Q, x6 c$ I
    Title 考虑交易费的投资组合模型;
    5 B1 y5 R3 x& ], \0 D1 |2 o6 bSETS:   
    " B$ Y5 P# I" W) C    STOCKS/ A,  B,  C/: C,Mean,X,Y,Z;   
    % T; k/ S0 K1 a* t3 S2 r  O) j; H    STST(Stocks,stocks): COV; + }2 Z8 i# j& m5 u$ [; N4 p
    ENDSETS
    1 X3 @# L' Y- M. Z, e; ZDATA: 2 n1 t% b" H, m$ J! ^6 d, _2 i
    ! 股票的初始份额;
    ! |$ c( b( g, g( W0 S/ ac=0.5  0.35  0.15;   & r& a6 p5 b5 I. E2 c2 Z
    TARGET=0.15;
    3 N0 H* P) g7 f5 EMean=@file(data2.txt);
    9 ?0 c" S. }  v  g+ k$ VCOV=@file(data3.txt);
    3 t" D' m8 A3 c" i9 ~ENDDATA & ?" N2 R. B+ c' b3 T; T
    [OBJ] MIN = @sum(STST(i,j): COV(i,j)*x(i)*x(j)); 5 l, F1 `2 a& I% x# z
    [ONE] @SUM(STOCKS: X+0.01*Y+0.01*Z) = 1;
    / j  w7 ?0 J* ]* b' i9 W8 u9 c$ d[TWO] @SUM(stocks: mean*x) >= TARGET; ' W, v. X" _3 [. W9 Y$ k; c
    @FOR(stocks: [ADD] x = c + y - z); 1 t# e- q( P$ S" U6 B. p
    END
    : e% Z0 M6 D; f1 E0 z在这个LINGO模型中,股票C是原始集合“STOCKS”的一个元素,不会因为与集 合的属性C同名而混淆。这是LINGO新版本比LINGO旧版本的一个改进之处。 1 I8 S' G, L0 y* t  \
    - X) m8 z  _" b3 ^& E

    9 N9 z( y+ g% Y' I+ E
    1 {1 |! g9 k! c! E利用股票指数简化投资组合模型
    : R% `5 u0 l8 }3 e* v例8  继续考虑例5(期望收益率仍定为15%)。在实际的股票市场上,一般存在 成千上万的股票,这时计算两两之间的相关性(协方差矩阵)将是一件非常费事甚至不可能的事情。例如,1000只股票就需要计算  =499500  个协方差。能否通过一定方式避免协方差的计算,对模型进行简化呢?例如,例5中还给出了当时股票指数的信 息,但我们到此为止一直没有利用。我们这一节就考虑利用股票指数对前面的模型进行 修改和简化。
    ! d. g4 o" n" t9 T' V1 x. B- V. `6 l, t) V: m& r, l
    (1)问题分析
    : x8 M! v- I3 Y/ e& J5 Z
      u6 L9 S3 c1 `, M) E可以认为股票指数反映的是股票市场的大势信息,对具体每只股票的涨跌通常是有 显著影响的。我们这里简单地假设每只股票的收益与股票指数成线性关系,从而可以 通过线性回归方法找出这个线性关系。
    4 a1 V" ?; ?6 P& O; {' m- }
    0 H7 {1 R1 z8 n# f(2)线性回归
    * \" _, I9 Y. S0 J1 w/ N2 }- Z$ D2 u$ W+ F' z8 K
    具体地说,用M 表示股票指数(也是一个随机变量),其均值为 ,方差为 。根据上面的线性关系的假定,对某只具体的股票i,其价值 (随机变量)可以表示成 : @' M! V- y; p5 v& L5 w8 T
                                                                 (9)! H% @9 N2 q% u6 ^; z

    : Z" N/ E7 {: D. U其中  和  需要根据所给数据经过回归计算得到,  是一个随机误差项,其均值为 ,方差为 。此外,假设随机误差项 与其它股票 j( j ≠ i )和股票指数M 都是独立的,所以  。 5 o5 z9 c1 k% @5 _1 m' T
      ]) ]3 ~7 F  \4 U' h
    先看看如何根据所给数据经过回归计算得到 和  。记所给的12年的数据为,( k =1,2,..,12),线性回归实际上是要使误差的平方和小,即要解如下优化问题: 6 m; _2 q, B: J9 c8 g) A( w- j
      M5 I4 b4 v2 q
                        (10)5 J+ W( N2 p1 O! `
    3 c9 c% b$ m% h: o5 L' Y/ E) @
    对这里给出的三种股票,可以编写如下LINGO程序求出线性回归的系数  (同时也在计算(CALC)段计算M 的均值 和方差   ,标准差  的值): 3 w9 L4 p' }3 s  f" G

    7 t6 D3 L* P9 i" K0 CMODEL:
    & G7 W6 \, ?7 J# n. uTitle 线性回归模型;
    0 l  A& H, ^. \' ?- U9 g3 w- z( s, iSETS:    YEAR/1..12/:M;   
    ( R$ N8 a7 a" e  |( lSTOCKS/A,  B,  C/: u, b, s2, s;   
    3 |9 S6 X8 [' e) [& jtemp/1..5/;       : ]& D7 w; C8 x$ j$ U2 R
    tmatrix(YEAR,temp):tm;   
    / |  w- p. f( ^8 T  i) H% Ilink(YEAR, STOCKS): R, e; & U- h! G! [1 [( y
    ENDSETS
    6 O" Z4 P. J" i4 x) nDATA:    - q$ P0 R1 _5 O5 k5 _3 j
    num=?;   
    ! ~$ q6 F, N( U6 r/ ztm =
      L4 ]$ D1 G) D5 X0 `1943           1.300          1.225           1.149           1.258997
    - I. w8 g! `" c+ J7 _1944           1.103          1.290           1.260           1.197526 ( H9 v7 D% _& V/ r, w& i/ p
    1945           1.216          1.216           1.419           1.364361
    , |# |5 H- `) ?& p1946           0.954          0.728           0.922           0.919287
    ; Y- u$ `/ M0 Z: v- @- [9 Q0 \1947           0.929          1.144           1.169           1.057080
    * X' [6 v0 S. f/ o1948           1.056          1.107           0.965           1.055012
    9 T- o% Z0 S$ ^  @8 A3 y* h1949           1.038          1.321           1.133           1.1879254 V( w- s% y  N/ G2 x; [) I
    1950           1.089          1.305           1.732           1.317130 : S- n# K3 s9 ^$ g: ?9 `
    1951           1.090          1.195           1.021           1.240164 4 J0 O4 i/ V, P6 r$ [* e
    1952           1.083          1.390           1.131           1.183675
    * N5 }- e. O: T5 S( q$ _1953           1.035          0.928           1.006           0.990108    % \# r1 w4 K. J' p; x) J, U
    1954           1.176          1.715           1.908           1.526236; ) F- C7 R; q4 I$ n" h, g
    ENDDATA ; d3 u2 ]2 N- @7 s7 I& O
    CALC:
    % {. U1 A) S2 N) H! d7 B@for(tmatrix(i,j)|j #ge#2 #and# j #le# 4:R(i,j-1)=tm(i,j)); 0 I& s! @' p# e, w/ e# P, i, K$ l' y9 V
    @for(tmatrix(i,j)| j #eq# 5:M(i)=tm(i,j));
    % I; I. f8 x, e( C+ P. R. Mmean0=@sum(year: M)/@size(year); 0 t6 X7 W- F- s0 O7 M
    s20=@sum(year: @sqr(M-mean0)) / (@size(year)-1);
    " `9 b# l& b7 p: r; e3 Cs0=@sqrt(s20); $ O! ^9 C6 N; d- }  O* C0 a
    ENDCALC 0 n+ N/ H# o  U5 ~" O
    [OBJ] MIN = @sum(stocks(i)|i#eq#num: s2(i));
    1 s! ^; z" c% @2 }5 z@for(link(k,i)|i#eq#num: [ERROR] e(k,i) = R(k,i)-u(i)-b(i)*M(k)); @for(stocks(i)|i#eq#num:[VAR] s2(i)=(@sum(year(k): @sqr(e(k,i))) / (@size(year)-2)); [STD] s(i)=@sqrt(s2(i)) );
    , p4 u1 R: `5 _@for(stocks: @free(u);@free(b) );
    1 p8 j, x$ f* U0 w; Z6 Q+ @@for(link: @free(e) );
    2 C3 C3 H- v- I# z. r( VEND 2 l2 a4 @; k6 Q5 _8 z" S" H2 P: z  Q
    对上面的这个程序,请注意以下几点:) A. w$ Y2 u# G

    7 x. ]* |, ~5 p( E" Ci)在CALC段直接计算了M 的均值  和方差  (为了使这个估计是无偏估计,分母是11而不是12)以及标准差
    ( f, [. e' c: V* y- f; U, n  r% T$ x
    ii)程序中使用了两个常用的数学函数:平方函数@sqr和平方根函数@sqrt。) P; G+ _4 ~9 @0 k" T

    ; [$ S3 E, y4 D9 y! f9 B0 z+ S7 G; [iii)除了计算回归系数外,我们同时估计了回归误差的方差 和标准差   。为了使这个估计是无偏估计,计算 时分母是10而不是11或12,这时因为此时估计了两个参数,自由度少了两个。% N4 `# d% u0 ]
    2 u! L( X! u/ ]! f5 q6 z; Q- t
    iv)@free(u),@free(b),@free(e)三个语句一定不能少,因为这几个变量不一定是非负的。
    0 B6 u/ i1 K8 p1 O8 V0 _9 J* J( m9 v2 K# ?
    v)DATA段定义了一个变量num,并用“num=?”语句表示其具体值需要由使用者在程序 运行时输入。变量num的作用是控制当前对哪只股票进行线性回归(num=1,2,3分别对应于股票 A , B, C)。
    + E0 z7 Q& G7 I5 v
    : T1 d: N' E0 j. O# b" Q% Yvi)其实,这个问题也可以对三只股票的回归不加区分,即放在同一个模型中同时优化(相应地,只需要去掉上面程序中的控制变量num和所有的过滤条件“i#eq#num”),不 过这样就会增加变量的个数,我们不建议大家那样做。也就是说,对于能够分解成小规模 问题的优化问题,好一个一个分开做,这样可以减少问题规模,有助于求到比较好的解。) W& J. ?/ ^% s! l* u1 b

    / A1 {- S, J+ ^- {! J% V, I) m运行上述LINGO程序,得到的计算结果为:股票指数M 的均值 = 1.191458  ,方差为 =0.02873661  ,标准差为 = 0.1695188  ;对股票 A,回归系数 = 0.5639761, = 0.4407264 ,误差的方差  = 0.00574832  ,误差的标准差  = 0.07581767  。
    1 r6 Z( I0 U/ Q1 ]$ R" N' A; e7 p
    2 b' {! V$ a3 L7 W7 @( Q同理(运行时输入num=2或3),可以得到:对股票B ,回归系数 = - 0.2635028 , = 1.239799 ,误差的方差 = 0.01564263  ,误差的标准差 =  0.1250705  。对股票C ,回归系数  = - 0.5809592 3  , = 1.523798  ,误差的方差 = 0.03025165 ,误差的标准差 =0.17393  。
      m. ]8 E8 B/ I# n+ t* |0 d, J# D9 ~1 ?) l0 t
    (3)优化模型   * H. T+ k/ L: B$ i/ m+ v1 n
    % b7 h7 Z" h% t2 X1 M
    现在,仍用决策变量 分别表示投资人应当投资股票 A ,B ,C 的比例,其中
    ; d4 s0 u( B& p$ K2 B2 H& x) {) Q9 R0 [0 p; y2 i
                  (11)0 |9 ^$ R+ }8 x/ z# N7 }: Y
    ) K* @* w( y9 l/ c  N
    此时,与第1节的讨论类似,对应的收益应该表示为 / f/ d/ U9 O! O) [; p
    7 S8 Z% d! R/ E0 p8 ~
    - k+ p0 h# s% v- W6 l# R5 D% X0 P
    9 {" z' {- u+ M6 g/ W
    (4)模型求解# \+ `" Y  x+ P" _7 k

    6 Y, v" n/ q9 v为了数据传递方便,我们把三个回归模型同时计算。LINGO程序如下
    % w* [1 W- ^: d1 T$ v' g
    0 `: G9 Y/ A  f$ W& U7 p$ U8 MMODEL:
    ! Z! C) B# _) c5 BTitle 线性回归模型;
    . u' A" {- b# [( A6 @SETS: # U& ]2 J6 B* A. `
        YEAR/1..12/:M;   
    + ~9 S) l* r$ v6 U/ c) q& y    STOCKS/A,  B,  C/: u, b, s2, s;    ' {% L, K3 T8 ]/ b. ?4 v; i* g
        temp/1..5/;       ) m4 F; w; Y& _1 Z6 ]: K) q/ w# K
        tmatrix(YEAR,temp):tm;   
    9 P8 @1 [" X$ i' p. m    link(YEAR, STOCKS): R, e; $ |$ c- M8 l, ]8 @; X1 S
    ENDSETS
    # X6 t! B" d& D) JDATA:
    0 E0 L$ j* e& ]; M! |tm =
    4 k$ e% B* W# g1943           1.300          1.225           1.149           1.258997
    - Q$ h& H6 b  d$ S" H7 ~% s0 M1944           1.103          1.290           1.260           1.197526
    % ~- w# N- M. N1 K8 x1945           1.216          1.216           1.419           1.364361
    " b% x: u4 S1 q/ T, _1946           0.954          0.728           0.922           0.919287 7 R' A7 r/ D" |
    1947           0.929          1.144           1.169           1.057080
    6 I9 Y  o/ K/ A1948           1.056          1.107           0.965           1.055012
    4 p' x# j4 G* v- ?5 i8 R1949           1.038          1.321           1.133           1.187925
      f& |7 |! y: _/ G1950           1.089          1.305           1.732           1.317130
    " G. j) W6 E6 V7 @) _$ _1951           1.090          1.195           1.021           1.240164
    : [) {/ J% v  w  U( D1952           1.083          1.390           1.131           1.183675 # I$ A8 z+ ^: i
    1953           1.035          0.928           1.006           0.990108    % F) V* m  R1 ?0 X+ |" Z
    1954           1.176          1.715           1.908           1.526236; @text('data4.txt')=s2; 1 Q; F4 K- o0 y' K+ q( ?2 u8 n; A
    @text('data5.txt')=u; - ^1 t3 v* ], a1 j# L
    @text('data6.txt')=b; , S. a3 y) A1 e( I: f1 D
    ENDDATA / B4 v; d. I6 C# {# [
    CALC: + C. R2 k; L8 h  B) V
    @for(tmatrix(i,j)|j #ge#2 #and# j #le# 4:R(i,j-1)=tm(i,j));
    ) J2 z2 Y# W; I) Q- O@for(tmatrix(i,j)| j #eq# 5:M(i)=tm(i,j));
    0 J3 V7 |- T/ h( `/ l; B* Lmean0=@sum(year: M)/@size(year);
    4 F5 ]$ K+ w6 Ns20=@sum(year: @sqr(M-mean0)) / (@size(year)-1); + N0 B: I& [4 Y3 w5 S. q, F- j
    s0=@sqrt(s20); . Y4 s; d; W/ S( R  p
    ENDCALC ! V$ N; [# T* @8 b
    [OBJ] MIN = @sum(stocks(i): s2(i));
    % v% I5 O$ v) t8 h9 t, g@for(link(k,i): [ERROR] e(k,i) = R(k,i)-u(i)-b(i)*M(k)); % d: U( l, [7 |/ c4 Z$ Y' v4 S: Y
    @for(stocks(i):[VAR] s2(i)=(@sum(year(k): @sqr(e(k,i))) / (@size(year)-2));
    # V+ ^; ^2 P5 G, J, h7 A- x9 \[STD] s(i)=@sqrt(s2(i)) );
    # E( @2 R: l, x4 C  b@for(stocks: @free(u);@free(b) );
    0 f+ `( h0 y% d6 i& F. C) g@for(link: @free(e) );
    ! u  m5 p& v  JEND8 _+ g2 `& b  b/ T7 ?5 p. M' a
    上述二次规划的LINGO程序如下:
    # K$ T4 u1 K8 {4 Q
    1 D" w: f. v2 b7 F' ?MODEL:
    0 y5 i% _  ]' |8 w+ B8 F, mTitle 利用股票指数简化投资组合模型; $ b) V" R0 o3 F2 c5 w8 F8 G, t
    SETS:    8 q  l4 H- v$ @, X7 X
        STOCKS/A,  B,  C/: u, b, s2, x;
    6 p, C+ o$ y$ V; H& g2 B. B4 [ENDSETS
    0 e+ r9 B' _3 U' BDATA:
    / i* j1 x$ x& a) V. S7 d! dmean0=1.191458;
    , S' K) }  i" j' c. K3 }! Qs20 = 0.02873661;
    # e) @2 W1 |: W! x; `$ R9 Ds2  =  @file(data4.txt);
    / o* _' R1 _, c& p0 N3 Q1 }u   =  @file(data5.txt); . b# g2 G( Y9 h' y/ `2 N6 x3 P
    b   =  @file(data6.txt); , t. z! K" ?6 G2 D  B! @0 }
    ENDDATA 1 ^  ^( H2 z, ]0 _
    [OBJ] MIN = @sum(stocks: (@sqr(b)*s20+s2)*@sqr(x));  
    7 h) d% ~0 c8 b. ?3 _9 R5 K@sum(stocks: x)=1;
    / U: B5 X3 q% I: q  @. [8 C@sum(stocks: (u+b*mean0)*x)>1.15;
    7 p8 [4 `( O2 C; a9 U9 FEND 3 a3 k# r3 M  h/ b2 O6 [: ~
    计算结果为,最后的持股情况是:A大约占初始时刻总资产的54%,B 占27%,C 占19%。这个结果与例5的结果是不同的。 6 u- R( s8 ]4 S
    5 ~# S+ A7 Z1 W3 o9 ^- G. f( s
    其它目标下的投资组合模型4 L5 x2 u3 ?7 H0 I
    前面介绍的模型中都是在可能获得的收益的数学期望满足一定低要求的前提下, 用可能获得的收益的方差来衡量投资风险,将其作为小化的目标。这种做法的合理性 通常至少需要有两个基本假设:2 R3 u0 c$ T1 [( ~6 P

    2 b7 l$ H1 j9 Z( m1 ?(1)可能获得的收益的分布是对称的(如正态分布)。因为这时未来收益高于设 定的低要求和低于设定的低要求的数量和概率是一样的。可惜的是,实际中这个假设往往难以验证。" {2 \6 J0 E: {; C1 E( y
    , h) S. n  W# u) r3 {
    (2)投资者对风险(或偏好)的效用函数是二次的,否则为什么值选择效益(随 机变量)的二阶矩(方差)来衡量风险使之小化,而不采用其它阶数的矩? 一般来说,投资者实际关心的通常是未来收益低于设定的低要求的数量(即低多 少)和概率,也就是说更关心的是下侧风险(downside risk)。所以,如果分布不 是对称的,则采用收益的方差来衡量投资风险就不一定合适。为了克服这个缺陷,可以 用收益低于低要求的数量的均值(一阶矩)作为下侧风险的衡量依据,即作为最小化 的目标。此外,也可以采用收益低于低要求的数量的二阶矩(即收益的半方差, semivariance)作为衡量投资风险的依据。其实,半方差计算与方差计算类似,只是只有当收益低于低要求的收益率时,才把两者之差的平方记入总风险,而对收益高于低要求的收益率时的数据忽略不计。这方面的具体模型这里就不再详细介绍了。 下面介绍一个与上面这些优化目标完全不同的投资组合模型,这个模型虽然很简单,但却会产生一些非常有趣的现象。) X6 ?, I$ d& l5 N9 k

    % k6 {  D2 g) n( X' c; Q例9  假设市场上只有两只股票 A, B可供某个投资者购买,且该投资者对未来一年的股票市场进行了仔细分析,认为市场只能出现两种可能的情况(1和2)。此外, 该投资者对每种情况出现的概率、每种情况出现时两只股票的增值情况都进行了预测和 分析(见表7,可以看出股票 B A、 的均值和方差都是一样的)。该投资者是一位非常 保守的投资人,其投资目标是使两种情况下小的收益最大化(也就是说,不管未来发生哪种情况,他都能至少获得这个收益)。如何建立模型和求解? - ]8 e  i0 Q: \6 \8 j2 R
    ) O# R6 }' @7 Q' V/ g; c

    " L& ^6 \" A- l" k% f
    . e/ Y5 o$ w" X. |* p! j  (1)优化模型与求解
    : T) z! y6 C8 F" y2 t. Y7 [- f' T" g3 a
    设年初投资股票A、 B 的比例分别为 ,决策变量显然应该满足 ) e  @' J! [+ a/ P/ F
                   (12) / ]- j" t) ~# J! ?$ M9 q/ N
    此外,使最小收益最大的“保守”目标实际上就是希望:
    / g# h# b- |5 ]8 ?; i4 D1 W! O/ ~& k& t. ^* z

    1 q4 e# H( \: R4 ^0 ]9 ?: i& @; R引入一个辅助变量  8 K% @7 y! j; Y2 P3 h. D7 }

    $ c$ _1 Z3 s8 O( ?4 a  4 Y% t- p! x1 s' ]& h

    " N  [9 E8 Y1 K# ~3 K9 J( q这个模型就可以线性化为  
    ! P. S( `( b1 m1 E/ N4 s" C
    ! j5 E( {0 e' t5 X' e
    5 c* j  s4 V. X" m
    1 N9 j. i. K9 t/ `编写LINGO程序如下:  
    # R' O: G) L; @8 O4 v
    2 o1 K% \( F0 s1 Dmodel:
    * `, b7 O: P' @7 a# `sets: 4 S/ J4 p8 G3 P$ J; V# m
    COL/1..2/:x;
    & G) K8 a6 d6 k& hROW/1..2/; ; e6 S2 e" c/ k8 e# d; K9 `
    link(ROW,COL):a;
    1 v" v+ y/ ]3 `# e$ `7 Fendsets 4 p( s5 N4 f, Z( H* i2 `
    data:
    / J& G# V" M5 G* D  v; g) Ja=1 1.2 1.5 0.7; 8 N  k8 n$ x) }+ M& a
    enddata + a! ~% T; G9 R
    max=y;
    3 C; x1 ~. M- M& g. ]( F@sum(COL:x)=1; 5 ?( }* J& e. j& Q# E& |
    @for(ROW(i)sum(COL(j):a(i,j)*x(j))>y);
    ( ~5 a' N' H5 X' a5 ~+ JEnd
    # N: h9 D: P* ~$ V  W可见,此时应该投资 A、B 股票各50%,至少可以增值10%。
    $ E0 Y$ s6 u3 d2 d% S
    8 i8 M8 P" B& u( K. ^7 p# P(2)讨论
    * W& v& S  j: |+ E4 ]0 ~6 q
    * p6 \3 J, G5 U. D' B) V. M现在,假设有一位绝对可靠的朋友告诉该投资者一条重要信息:如果情形1发生, 股票B 的增值将达到30%而不是表7中给出的20%。那么,一般人的想法应该是增加 对股票B 的持有份额。果真如此吗?这个投资人如果将上面模型中的1.2改为1.3计 算,将得到如下结果:  =0.5454545  , =  0.4545455  , y=1.136364 。2 p6 S! [, J0 c& Y6 u

    $ }( ]+ @! k* G/ _; K$ d也就是说,应该减少对股票B 的持有份额,增加对股票 A的持有份额。这真是叫 人大吃一惊!这相当于说:有人告诉你有某只股票涨幅要增加了,你赶紧说:那我马上 把这只股票再卖点吧。之所以出现如此奇怪的现象,就是由于这个例子中的目标的特殊 性引起的:我们可以看到新的解可以保证增值达到13.6364%,确实比原来的10%增 加了。
    ) M+ R0 M4 T' I0 r# f% f: z& N; E& c/ k1 x/ r4 T, L0 Z7 ?. j
    最后需要指出:我们上面所有关于投资组合的这些讨论基本上只是纯技术面的讨 论,只利用历史数据来说话,认为历史数据中包含了引起股票涨跌的所有因素。在实际 股票市场上,影响股票涨跌的因素可能有很多(如政策变化、银行加息、能源短缺、技 术进步等),未来不长时间内可能发生的一些重大事件很可能以前没有发生过,因此也 不可能体现在历史数据中。所以,进行投资选择前,还应该进行基本面分析,需要对未 来的一些重要影响因素、重大事件发生的可能性及其对每种股票涨跌的影响进行预测和分析,最后综合利用历史数据和这些预测数据,决定投资组合。如何将这些预测数据与 历史数据一起使用,建立相应的投资组合模型,这里就不再更多地介绍了。这方面的模 型有很多,有兴趣的可以继续查阅相关的专业书籍和研究文献。' \7 x3 l6 m/ S4 K& r

    ; C% e- A7 j4 I! r4 ~$ K/ V- p+ g2 k; l习题:9 b- E4 e5 a. c% k

    0 X$ B; I( r! D- ^/ ~0 v某银行经理计划用一笔资金进行有价证券的投资,可供购进的证券及其信用等 级、到期年限、收益如表11所示。按照规定,市政证券的收益可以免税,其它证券的 收益需按50%的税率纳税。此外还有以下限制:
    / j9 M6 \/ c- J6 V- R) b
    # h9 M* o% b; Ei)政府及代办机构的证券总共至少要购进400万元;- b( U. N; ?9 i3 A' U; N
    : e8 ~. @  @* }; c
    ii)所购证券的平均信用等级不超过1.4(信用等级数字越小,信用程度越高);9 t, t; `5 g# I5 u/ o
    ' g+ {2 Z" J! A' u
    iii)所购证券的平均到期年限不超过5年。
    & F2 |9 @  h0 D' _* h% _
    * |2 X2 P& Y- `& C# _7 `# Q: n! R8 K7 a. V

    * U: n5 v% I" I5 j (1)若该经理有1000万元资金,应如何投资?
    ' R7 t$ K2 l8 m; }7 u) z5 V7 Z
    ! z* i7 `: c8 K: e: i: i, \(2)如果能够以2.75%的利率借到不超过100万元资金,该经理应如何操作?7 w* ^0 Q7 W4 g. A2 E% `" B

    3 z- e' |5 U% _7 i  N$ }(3)在1000万元资金情况下,若证券 A的税前收益增加为4.5%,投资应否改变? 若证券C 的税前收益减少为4.8%,投资应否改变?   |! N/ o' w5 T2 s, q
    ————————————————
    5 [5 y) |# t9 O1 I版权声明:本文为CSDN博主「wamg潇潇」的原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。6 F% U! }+ O- [' h: w! M
    原文链接:https://blog.csdn.net/qq_29831163/java/article/details/89406997
    3 s: w) a5 k5 R5 y  c* F, T1 ?
    * m5 r# O% p' A6 j* W5 T+ K/ Z. {8 `% g: _
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-8-1 05:38 , Processed in 0.438942 second(s), 51 queries .

    回顶部