计算beta分布并画图(1)
2018-06-18 02:52:32来源:未知 阅读 ()
贝塔函数与伽马函数的关系
详细推导过程见LDA漫游指南
java的org.apache.commons.math3.special.Gamma封装了Gamma函数可以直接使用
本文首先计算出B(3.9,2.9)与B(3.9,5.3)的100个散点,然后绘制出Beta分布图像
import java.io.BufferedWriter; import java.io.FileWriter; import java.io.IOException; import org.apache.commons.math3.special.Gamma; public class betadistr { public static double Gammafun(double alpha){ return Gamma.gamma(alpha); } public static double Bfunction(double alpha,double beta) { double totgamm = Gammafun(alpha+beta); double addgamm = Gammafun(alpha)*Gammafun(beta); return addgamm/totgamm; } public static void valued() throws IOException { BufferedWriter writer = new BufferedWriter(new FileWriter(filename)); int K=2; double[][] m=new double[K][100]; double b[][]= new double[2][2]; b[0][0]= 3.9; b[0][1]= 2.9; b[1][0]= 3.9; b[1][1]= 5.3; for (int i = 0; i < K; i++) { for (int x = 0; x < 100; x++) { double y = (double)x/100; double f=Math.pow(y,b[i][0]-1)*Math.pow(1-y, b[i][1]-1); double z=Bfunction(b[i][0],b[i][1]); m[i][x]=f/z; writer.write(m[i][x] + "\t"); System.out.println(m[i][x]); } writer.write("\n"); } writer.close(); } public static void main(String[] args) throws IOException { valued(); } }
标签:
版权申明:本站文章部分自网络,如有侵权,请联系:west999com@outlook.com
特别注意:本站所有转载文章言论不代表本站观点,本站所提供的摄影照片,插画,设计作品,如需使用,请与原作者联系,版权归原作者所有
下一篇:JWT 简介
- 聊聊微服务架构及分布式事务解决方案! 2020-06-10
- 作为一个面试官,我想问问你Redis分布式锁怎么搞? 2020-06-10
- Java实现圆柱体表面积和体积的计算 2020-06-08
- 计算机基础到底是哪些基础?为什么很重要! 2020-06-08
- 分布式锁没那么难,手把手教你实现 Redis 分布锁!|保姆级教 2020-06-08
IDC资讯: 主机资讯 注册资讯 托管资讯 vps资讯 网站建设
网站运营: 建站经验 策划盈利 搜索优化 网站推广 免费资源
网络编程: Asp.Net编程 Asp编程 Php编程 Xml编程 Access Mssql Mysql 其它
服务器技术: Web服务器 Ftp服务器 Mail服务器 Dns服务器 安全防护
软件技巧: 其它软件 Word Excel Powerpoint Ghost Vista QQ空间 QQ FlashGet 迅雷
网页制作: FrontPages Dreamweaver Javascript css photoshop fireworks Flash