matalb——生成离散非均匀随机数

今日感冒,状态非常差,连锁故障风险评估的程序仍然没有调试出bug的问题所在。在和Yu XingBin的交流中发现了程序设计中另外一个问题,会导致程序计算的指标出现较大偏差。

一、问题

连锁故障仿真程序中,第一步就是需要随机选择一条故障线路。起先我认为按照均匀分布产生[1 n]之间的任一随机数(n为总线路数),结果发现仿真结果出现母线分离的概率非常大,远远超过了Yu XingBin论文中的结论。疑惑不解,遂于之邮件交流,得知两点信息:

1.输电线路的故障概率远远大于变压器线路故障概率。因为输电线路很长,并且暴露在外部环境中,易受天气等因素的影响,其故障概率显然大于变压器。

2.输电线路的故障概率近视的认为与线路的长度成正比。

根据以上信息分析可以得知,对于固定拓扑结构的电力网络,故障线路的概率是随机的、离散的,而不是均匀分布的。故障线路的选择应该按照每条线路的故障概率进行抽样。所以需要实现离散非均匀分布的随机整数抽样。

二、解决方案

网上现成有一个解决方案如下: function y=dis_rand(x,p,n)

%dis_rand 产生离散分布随机数

% x:可能取值

% p:取值概率

% n:拟生成随机数的数目

cp=cumsum(p);

y=zeros(1,n);

for i=1:n

y(i)=x(sum(cp<=rand(1))+1);

end

算法说明:利用matlab中的rand函数获取区间[0,1]上均匀分布的随机数,根据希望的分布律把[0,1]区间分段,各段的长度依次为随机数取各个可能值的概率。当那个均匀分布的随机数落在其中第i个小段上时,就把第i个可能取值输出。

例:y=dis_rand([1 :3],[0.4 0.3 0.4],1)  即按照0.4 0.3 0.4的概率从1,2,3三个数字中产生一个随机数。

matlab坐标轴非均匀分布,matalb——生成离散非均匀随机数相关推荐

  1. 使用SQL生成非均匀随机数

    参考: http://www.cnblogs.com/CareySon/archive/2012/07/11/GenerateNURNsUsingSQLServer.html 正如"随机数的 ...

  2. 扩增子分析解读4去嵌合体,非细菌序列,生成代表性序列和OTU表

    写在前面 之前发布的<扩增子图表解读>系列,相信关注过我的朋友大部分都看过了(链接直达7月文章目录).这些内容的最初是写本实验室的学生们学习的材料,加速大家对同行文章的解读能力. < ...

  3. 用力和应变片计算弹性模量_一种沿深度非均匀分布的残余应力测试计算方法与流程...

    本发明涉及一种沿深度非均匀分布的残余应力测试计算方法,属于残余应力测试技术领域. 背景技术: 机械零件在加工制造.装配过程中,不可避免地产生残余应力,对零件的强度.刚度.疲劳等机械性能产生严重的影响. ...

  4. 成功解决pillow模块内生成的非UTF-8格式(比如label中有汉字时)

    成功解决pillow模块内生成的非UTF-8格式(比如label中有汉字时) 目录 解决问题 解决方法 解决问题 解决pillow模块内生成的非UTF-8格式(比如label中有汉字时),需要进行格式 ...

  5. matlab—赋值具有的非单一 rhs 维度多于非单一下标数

    刚开始学习接触matlab,有许多地方有疑问,想在学习的过程中做一些简单的记录与总结. 下面是想创建一个元胞数组,2*2的. B(1,1)={'This is the first cell'}; B( ...

  6. MATLAB 赋值具有的非单一 rhs 维度多于非单一下标数

    Tstr=char(Policy); for i=1:nfor j=1:nif Policy(i,j)==1Tstr(i,j)='\uparrow';elseif Policy(i,j)==2Tstr ...

  7. PHP多文字,二维码(动态、非动态)生成海报方式

    PHP多文字,二维码(动态(支持带logo).非动态)生成海报方式 1.下载二维码插件Phpqrcode,地址 [https://sourceforge.net/projects/phpqrcode/ ...

  8. springfox源码_【开源项目】springfox-bridge:随心所欲地为非restful接口生成API文档...

    一.引言 目前,利用swagger框架为restful接口编写API文档非常流行,在spring web项目中,利用springfox+swagger更是可以通过注解的方式直接进行API文档的生成,这 ...

  9. matlab离散方波信号,MATLAB - 在ODE15s函数中生成方波信号?(MATLAB - Generate square signal in ODE15s function?)...

    MATLAB - 在ODE15s函数中生成方波信号?(MATLAB - Generate square signal in ODE15s function?) 我试图从函数内部生成一个方形信号,通过O ...

  10. Matlab学习记录-矩阵的生成

    Matlab学习记录-矩阵的生成 matlab中生成矩阵有两种方式: 1.中括号加分号. 中括号表示矩阵,分号表示分行:每一行之间的元素可以用逗号分开也可以用空格分开: 例如 a=[1 2 3; 4 ...

最新文章

  1. “Account is trusted for delegation is not shown || cannot see the Delegation tab
  2. 【Axure原型图】—— 4. Tab Control(选项卡)
  3. c语言图形库怎样画一条弧线,自己写的一个图形库,有点问题.来请教大家了
  4. 安装 Cacti 监控
  5. 列了一些自己会但是不怎么精通的编程语言和知识,做个记录,空余时间加强学习...
  6. java 点餐界面_Java小项目点餐系统(二)之服务端 | 学步园
  7. C#多线程与UI响应 防止界面假死不响应(子线程创建的窗体获取消息响应用Application.DoEvent )
  8. vue中使用cookies和crypto-js实现记住密码和加密
  9. UVa 129 - Krypton Factor(回溯法)
  10. 【从0到1,搭建Spring Boot+RESTful API+Shiro+Mybatis+SQLServer权限系统】05、Shiro集成
  11. Go 如何利用 Linux 内核的负载均衡能力
  12. Slog57_玩转NPM之NPM_package的制作发布和使用
  13. SQL Server查询执行计划–基础
  14. Python路飞学城老男孩内部书籍,Python全栈开发实战pdf
  15. vue 登录页qq快速登录功能
  16. 树莓派有线网络设置_树莓派通过网线实现与电脑共享网络
  17. 关于最佳观影距离分享 投影仪测量的那些事
  18. macbook pro window蓝牙驱动 下载
  19. Servlet+JSP一文完结
  20. MLDLRL:MLDL相关概念的原始英文解释——对理解最初的ML和DL的相关概念的定义非常有用

热门文章

  1. STM32的ucgui移植
  2. python获取浏览器cookie_python3实现读取chrome浏览器cookie
  3. 射频电路设计中的热量分析
  4. NMAKE makefile内容和功能
  5. matlab 编写雷达波形,雷达信号处理+Matlab程序
  6. nodejs 游戏框架_Pomelo:网易开源基于 Node.js 的游戏服务端框架
  7. android抓取微信朋友圈,一种快速提取Android微信朋友圈数据的方法
  8. 易观分析联合中小银行联盟发布海南数字经济指数,敬请期待!
  9. Java基础语法结构
  10. [unity3d插件]插件效果以及下载