极值理论(一):极大值极限分布

1 为什么引入广义极值分布:

        考虑随机变量序列\left \{ X_{n} \right \}极大值分布:P(M_{n}\leq x)=F^{n}(x),当n \to \infty时,F^{n}(x)不一定是依分布收敛的(一般我们总是期望它是收敛的)。因此,我们引入广义极值分布来描述标准化M_{n}

2 广义极值分布(GEV):

对于标准化M_{n}

1+\xi x> 0

 \xi=0:Gumbel分布;\xi>0:Frechet分布;\xi<0:weibull分布

何谓标准化:对于实际M_{n},若存在实常数序列\left \{ c_{n} \right \}\left \{ d_{n} \right \},使得:

则称\left ( M_{n}-d_{n} \right )/c_{n}为标准化

由此引出三参数分布族H_{\xi ,\mu ,\sigma }(x)=H((x-\mu )/\sigma ):

P(M_{n}\leq x)=P((M_{n}-d_{n})/c_{n}\leq (x-d_{n})/c_{n})=H((x-d_{n})/c_{n})n \to \infty

 实际上,\mu =d_{n}\sigma =d_{n}

GEV分布图像:

定义与定理:

举个例子: 

 3 极小值收敛:

 4 极大值吸引域:

本节介绍不同\xi下分布的特征

1)Frechet情况:

先引入两个定义:

 一个定理:

 

由此可以推断,对于k> 1/\xi ,E(X^{k})=\infty  

Frechet情况本身即为厚尾分布,并且由上述定理可知,这些服从Frechet情况的均为无限高阶矩厚尾分布,因此也成为了EVT中研究最多的分布。

常见分布:

 

举个例子:

 2)Gumbel情况:

尾部指数衰减,分布函数服从MDA(H_{0})的正值分布的正阶矩是有限的,此类分布的尾部有较多变化,如薄尾的正态分布与后尾的对数正态分布,在区分其尾部行为时(与Frechet情况),需要搜集大量数据,而金融模型常常将本应是Gumbel情况误当作Frechet情况处理

常见分布:

3)weibull情况:

最不重要的情况,主要是因为其有限的分布右端点。但是信用风险模型会有例外

一个定理:

5 严平稳时间序列极大值:

 其中\theta为聚类规模的倒数

 

 注:H^{\theta }(x)\subseteq H_{\xi }(x)

6 如何实证:

实际上通过上述的讨论,我们可以得知无论是独立同分布序列还是严平稳时序的非标准化极大值的极限分布都是 H_{\xi ,\mu ,\sigma }(x)=H((x-\mu )/\sigma ) 类型,考虑用极大似然法估计其参数,首先假设数据分为m个大小为n的组,分别取m个组中的最大值M_{nj}( j=1,2,...,m),得到如下对数似然:

 满足1+\xi (M_{ni}-\mu )/\sigma > 0 

模型检验与置信区间:

  • 6
    点赞
  • 50
    收藏
    觉得还不错? 一键收藏
  • 6
    评论
在MATLAB中,寻找矩阵中的局部极大值通常涉及到信号处理或图像处理中的梯度计算和峰值检测。以下是一个简要的步骤概述: 1. **定义矩阵**:首先,你需要一个二维数组(矩阵)作为输入数据。 2. **计算梯度**:使用MATLAB的`gradient`函数或者`imgradient`函数计算矩阵的x和y方向的梯度,这将给出每个像素变化的方向和幅度。 3. **平滑处理**:为了减少噪声影响,对梯度图进行平滑(如使用`imfilter`函数中的高斯滤波器)。 4. **峰值检测**:利用`maxFilter`或者`imfindpeaks`函数查找局部极大值点。这些函数可以找到梯度图中的局部最大值及其坐标。 5. **验证极值**:确认检测到的点确实是局部极大值,而不是边缘或其他噪声造成的局部最大。 6. **设置邻域大小**:在实际应用中,可能需要指定一个邻域大小,只考虑这个大小内的像素值来确定是否为真正的极大值。 下面是一个示例代码片段: ```matlab % 假设你有一个名为data的矩阵 data = randn(100, 100); % 示例数据 % 计算梯度 [gradX, gradY] = gradient(data); % 平滑梯度 gradX_smoothed = imfilter(gradX, fspecial('gaussian', [5 5], 2)); gradY_smoothed = imfilter(gradY, fspecial('gaussian', [5 5], 2)); % 检测峰值 [maxVal, locs] = imfindpeaks(gradX_smoothed + gradY_smoothed, 'MinPeakHeight', 0.1, 'MaxPeakDistance', 10); % 邻域大小设为10 % 查看局部极大值点 figure; imagesc(data); hold on; scatter(locs(:,2), locs(:,1), 'r', 'filled'); ```

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 6
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值