一个简单的例子解释随机矩阵里的MP law,

#############################
##
## A example to show the MP law in random matrix theory.
##
#############################
library(ggplot2)## The empirical distribution
p <- 1000
n <- 2000
dat <- matrix(rnorm(p*n), nrow=p)
dat.cov <- cov(t(dat))
dat.eig <- eigen(dat.cov)
dat.val <- dat.eig$values
dat.dsty <- density(dat.val)## plot
df.dsty <- data.frame(x=dat.dsty$x, y=dat.dsty$y)
p <- ggplot(df.dsty, aes(x = x, y = y))
p + geom_line()  + ylab('density') + ggtitle("MP law")## The largest eigenvalue
p <- 1000
n <- 2000
K0 <- 100
lam <- list()
for (k in 1:K0){dat <- matrix(rnorm(p*n), nrow=p)dat.cov <- cov(t(dat))dat.eig <- eigen(dat.cov)dat.val <- dat.eig$valueslam[[k]] <- dat.val
}lam.max <- unlist(lapply(lam,max))
lam.min <- unlist(lapply(lam,min))######################
## plot the max and min eigenvalue
######################
plot(lam.max, ylim=c(0,3), main='extreme eigenvalues of sample covariance matrix', xlab='Id', ylab='extreme eigenvalue', col='red', type='p')
lines(lam.min,type='p', col='blue')######################
## Save the data
######################
lam.MP <- lam
save(lam.MP, file="D:/BinbinChen/Literature/Orcal Defence/Simulation.MP.Rdata")

转载于:https://www.cnblogs.com/BinbinChen/articles/3354881.html

MP law simulation相关推荐

  1. matlab产生均匀白噪声,各种分布白噪声的产生matlab.pdf

    各种分布白噪声的产生matlab 3.各种分布白噪声的产生 均匀分布白噪声的产生 1.物理方法 2.数学方法 线性同余法.联合法.反馈位移寄存器法 非均匀分布白噪声的产生 1.理论方法 反变换法. ...

  2. Poland’s ruling Law and Justice party is doing lasting damage

    扫码进入阅时即查模式 Once hollowed out, the rule of law is hard to restore FOR a glimpse of Poland under the p ...

  3. SOLIDWORKS Simulation攻略丨赫兹接触应力分析

    赫兹接触应力是弹性力学中的经典力学问题.其适用范围很广,最典型的比如:滚动轴承中滚珠与坐圈的接触,两个啮合齿轮在齿面上的接触,车轮与轨道的接触等. 以机械设计中直齿轮齿面接触在接触面上产生接触应力为例 ...

  4. sumo添加车辆_sumo demo 我们通过使用交通仿真软件SUMO(Simulation of 联合开发网 - pudn.com...

    sumo demo 所属分类:Windows编程 开发工具:C/C++ 文件大小:3KB 下载次数:15 上传日期:2019-03-07 16:12:12 上 传 者:tonythetiger1110 ...

  5. HP小型机MP简易使用手册

    第二章 访问MP 可以通过MP serial port 或者MP LAN port 登陆MP,默认的用户名/口令为:Admin/Admin 登陆过程: ************************ ...

  6. Lotus.Engine.Simulation.v5.06f 1CD(PC平台上用于模拟发动机运行表现的一款软件)

    Lotus Engineering Software产品: Lotus.Base.Engine.Analysis.Tools.v4.02g 1CD(发动机制造前的理论性研究分析工具) Lotus.Co ...

  7. Paschen law(帕邢/巴申定律)

    Paschen law:百度百科有二种解释帕邢定律和巴申定律,在国内文献中引用也存在二种称呼混用.个人习惯称呼为帕邢定律,可能是第一次的印象较深也可能是这个名字更霸气. 对产品设计基带工程师,帕邢定律 ...

  8. MP(Maxpressure)算法的信号控制仿真实例

    系列文章目录 利用SUMO进行交叉口速度引导 SUMO使用E2检测器获取信号交叉口车道信息和信号控制方案 SUMO轨迹图的绘制和traci接口的简单使用 SUMO利用转弯率构建车流rou文件 前言 博 ...

  9. MP实战系列(二)之集成swagger

    其实与spring+springmvc+mybatis集成swagger没什么区别,只是之前写的太不好了,所以这次决定详细写. 提到swagger不得不提rest,rest是一种架构风格,里面有对不同 ...

最新文章

  1. 索引块与数据块的区别
  2. [Leetcode总结] 98.验证二叉搜索树
  3. 7-11 堆栈操作合法性 (20 分)以及堆栈讲解(c语言实现)
  4. Unionid-微信开发学习
  5. 阿里云学生轻量级应用服务器安装MySQL
  6. My97DatePicker日历的平面显示,不是文本框点击事件后显示
  7. 数学建模算法与应用_《数学建模算法与应用》笔记【1】
  8. 警告: 程序集绑定日志记录被关闭。解决方法
  9. 学习 shell —— 相对路径转换为绝对路径
  10. UniqueIdentifier 数据类型
  11. 开始时间和结束时间,用例设计
  12. PAT 1055 集体照 python
  13. 解释一下积分变上限函数
  14. js做gif动画【MovieClip插件】
  15. 拆解玩具电池充电器:充久了可能会爆,廉价电路方案让人震惊!
  16. 文本文件与二进制文件区别 r 与 rb 方式 w 与 wb方式(windows)—————— 开开开山怪
  17. LXR( Linux超文本交叉代码检索工具)
  18. 【已解决】PEP 8: W292 no newline at end of file等相关PEP 8错误与警告
  19. 科学计算库Pandas的使用(附思维导图)
  20. chrome浏览器中自带input样式input:-internal-autofill-selected(修改input背景色)

热门文章

  1. 第16 17章节-Python3.5-Django知识点整理 15
  2. Linux内核【链表】整理笔记(1)
  3. 谈谈物联网产业4G转型中的隐忧与应对
  4. [译] Node.js 流: 你需要知道的一切
  5. Objective-C原理系列(一)
  6. 国内首个基于liftweb开发的项目上线运营
  7. Android多媒体学习一:Android中Image的简单实例。
  8. Linux 系统文件及结构
  9. Map And Reduce Join的使用案例
  10. VUE-CLI可选的配置文件vue.config.js