二项分布均值和方差的简单推导

   前一篇文章《二项分布》中说过,伯努利分布(也称为两点分布或0-1分布)是二项分布在n=1时的特例。我们先看伯努利分布的均值和方差的推导。

   根据离散型随机变量均值和方差的定义,若离散型随机变量X的分布列为:

Xx1x2...xi...xn
Pp1p2...pi...pn

   则称E(X)=x1*p1+x2*p2+...+...xi*pi+...+xn*pn为随机变量X的均值或数学期望,称为随机变量X的方差。

   伯努利分布的分布列为:

X01
P1-pp

   则根据离散型随机变量的均值和方差定义:
E(X)=0*(1-p)+1*p=p
   

D(X)=(0-E(X))2(1-p)+(1-E(X))2p=p2(1-p)+(1-p)2p=p2-p3+p3-2p2+p=p-p2=p(1-p)

   对于二项分布X~B(n,p)X表示的是n次伯努利试验中事件发生次数的随机变量。用Xi表示第i次伯努利试验中的随机变量,那么n次伯努利试验总的随机变量X可以表示成:

X=X1+X2+...+Xi+...+Xn

   一直没有找到满意的随机变量和、差、积、商的物理/几何/现实意义,如果有了解的朋友不妨留言,不甚感激。

   根据均值和方差的性质,如果两个随机变量X,Y相互独立,那么:

E(X+Y)=E(X)+E(Y)

D(X+Y)=D(X)+D(Y)

   对于二项分布X~B(n,p),每一次伯努利试验都相互独立,因此:

E(X)=E(X1)+E(X2)+...+E(Xi)+...+E(Xn)=p+p+...+p+...p=np

D(X)=D(X1)+D(X2)+...+D(Xi)+...+D(Xn)=p(1-p)+p(1-p)+...+p(1-p)+...+p(1-p)=np(1-p)

### 马科维茨均值-方差模型的数学推导 马科维茨均值-方差模型的核心在于通过优化投资组合中的权重 $ w_i $ 来实现风险最小化的同时满足一定的收益目标。以下是该模型的具体数学推导过程: #### 1. 投资组合的风险表示 对于由 $ N $ 种资产组成的投资组合,其总风险可以用协方差矩阵来描述。设每种资产的权重为 $ w_i $ ($ i = 1, 2, ..., N $),则投资组合的总体风险(方差)可写成如下形式: $$ \sigma^2(R_p) = \sum_{i=1}^{N} w_i^2 \sigma^2(R_i) + \sum_{i \neq j} w_i w_j \sigma(R_i, R_j) $$ 其中: - $ \sigma^2(R_i) $ 是第 $ i $ 种资产的方差; - $ \sigma(R_i, R_j) $ 是第 $ i $ 第 $ j $ 资产之间的协方差。 上述表达式可以通过向量矩阵的形式进一步简化为: $$ \sigma^2(R_p) = w^\top \Sigma w $$ 这里: - $ w = [w_1, w_2, ..., w_N]^\top $ 表示权重向量; - $ \Sigma $ 是 $ N \times N $ 的协方差矩阵[^4]。 #### 2. 收益约束条件 为了使投资组合达到特定的目标收益率 $ \overline{R}_p $,需引入以下约束条件: $$ \overline{R}_p = \sum_{i=1}^{N} w_i E(R_i) $$ 或者用向量形式表示为: $$ \overline{R}_p = w^\top \mu $$ 其中 $ \mu = [E(R_1), E(R_2), ..., E(R_N)]^\top $ 是期望收益率向量。 此外,还需满足权重之等于1的归一化条件: $$ \sum_{i=1}^{N} w_i = 1 \quad 或者 \quad \mathbf{1}^\top w = 1 $$ 这里的 $ \mathbf{1} $ 是一个长度为 $ N $ 的全1列向量。 #### 3. 最优投资组合问题建模 综合以上两点,最优投资组合问题可以表述为以下二次规划问题: $$ \min_w w^\top \Sigma w $$ 受制于两个主要约束条件: 1. $ w^\top \mu = \overline{R}_p $ 2. $ \mathbf{1}^\top w = 1 $ 这是一个典型的带等式约束的凸优化问题,可通过拉格朗日乘子法解决。 #### 4. 使用拉格朗日乘子法求解 定义拉格朗日函数为: $$ L(w, \lambda_1, \lambda_2) = w^\top \Sigma w - \lambda_1 (w^\top \mu - \overline{R}_p) - \lambda_2 (\mathbf{1}^\top w - 1) $$ 对 $ L $ 关于 $ w $、$ \lambda_1 $ $ \lambda_2 $ 分别取偏导并令其为零,得到一组必要条件: $$ \frac{\partial L}{\partial w} = 2\Sigma w - \lambda_1 \mu - \lambda_2 \mathbf{1} = 0 $$ $$ \frac{\partial L}{\partial \lambda_1} = -(w^\top \mu - \overline{R}_p) = 0 $$ $$ \frac{\partial L}{\partial \lambda_2} = -(\mathbf{1}^\top w - 1) = 0 $$ 由此得出三个关键方程: 1. $ \Sigma w = \frac{\lambda_1}{2} \mu + \frac{\lambda_2}{2} \mathbf{1} $ 2. $ w^\top \mu = \overline{R}_p $ 3. $ \mathbf{1}^\top w = 1 $ 这些方程联立起来即可求得最优权重 $ w $ 及对应的参数 $ \lambda_1 $ $ \lambda_2 $。 --- ### MATLAB 实现示例 MATLAB 提供了一个名为 `portcons` 的工具箱用于构建投资组合约束,并结合其他优化算法完成实际应用。下面是一个简单的代码示例展示如何利用 MATLAB 解决这一问题: ```matlab % 参数设置 Mu = [0.12; 0.18]; % 各资产的期望收益率 Sigma = [0.2 0.1; 0.1 0.3]; % 协方差矩阵 TargetReturn = 0.15; % 目标收益率 % 构造约束矩阵 Aeq beq Aeq = [Mu', ones(1, length(Mu))]; % 均值与单位向量构成的约束 beq = [TargetReturn, 1]; % 初始化优化器 f = zeros(length(Mu), 1); % 定义目标函数系数 H = Sigma; % 求解二次规划问题 [w, ~] = quadprog(H, f, [], [], Aeq, beq); disp('最优权重:'); disp(w); ``` 此代码实现了基于给定目标收益率寻找最优权重的过程。 ---
评论 8
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值