项目中遇到了用java计算二维标准正态分布累计函数。网上了查了好久才找到解决的方法,特此记录下来。
问题描述:求解二维标准正态分布累计函数M(a,b:ρ),其中,a,b分别为两个变量的最大取值。ρ为a和b的相关系数。 运用Matlabl里面的mvncdf进行计算,得到结果。
系统以及程序版本
Win 7(64 bits)
MATLAB 2016a(64 bits)
Eclipse LUNA(64 bits)
Java 1.7(64 bits)
1. 环境配置
主要是Java的环境配置。这里就不详细阐述了。
2. 利用MATLAB生成Jar包
这里以matlab中标准二维正态分布函数为例:
function [f] = calmvncdf(a,b,p)
f = mvncdf([a,b],[0,0],[1,p;p,1]);
end
(1) 查看MATLAB的JDK版本
MATLAB 2016a生成Jar包,一定要用JDK1.7,用JDK1.8会报错。
(2) 新建calmvncdf.m脚本
点击新建脚本
编辑器里面编辑自定义函数
保存为.m文件
(3) 在命令行窗口输入deploytool
(4) 弹出MATLAB Compiler窗口,选择library Compiler
(5) 在application type窗口选择Java package
(6) 在exported functions中点击加号按钮,在弹出的添加文件窗口中选择calmvncdf.m文件
(7) 在main file窗口选择runtime downloaded from web
(8) 选择Package按钮,开始打包
(9) 点击保存之后,显示正在编译
(10) 如图所示,则表示编译完成
如果失败,检查下,java环境是否为JDK1.7
(11) 编译完成之后,会形成四个文件,三个文件夹和一个html文件
(12) 打开for_redistribution_files_only文件夹,文件夹中的calmvncdf.jar就是打包成的jar包
(13) 其中for_redistribution_files_only\doc\html中的index.html是关于jar包的使用说明
(14) 拷贝出calmvncdf.jar和Matlab安装目录\toolbox\javabuilder\jar文件夹中的javabuilder.jar文件
3. 将Jar包添加到Eclipse工程中
将...\MATLAB\R2016a\toolbox\javabuilder\jar中的javabuilder.jar和生成的calmvncdf.ja包一同加入到工程中,并通过build path添加到Referenced libraries 中
4. 现在就可以在java里使用MATLAB中的函数计算了
public class TestMvncdf {
public static void main(String[] args) {
double a = 4.81;
double b = 2.78;
double p = 0.38;
Class1 MyClass;
try {
MyClass = new Class1();
Object ans[] = MyClass.calmvncdf(1,a,b,p);
System.out.println("ans[0]:"+ans[0]);
} catch (MWException e) {
e.printStackTrace();
}
}
}
结果:
ans[0]:0.9973
Matlab中的结果:
Java和Matlab中计算结果一致。
该方法适用于Java调用Matlab中的自带函数和自定义函数。
如果你要在没有安装matlab的机器上用java调用这些jar包,怎么办?也是可以的,不过,不过你要先在机器上安装matlab提供的MATLAB Compiler Runtime (MCR) ,
它的安装文件在D:\Program Files\MATLAB\R2016a\toolbox\compiler\deploy\win64下(你的安装目录下)文件名MCRInstaller.exe
参考资料:
1、http://www.ilovematlab.cn/forum.php?mod=viewthread&tid=265570