У меня есть n-мерный массив ppi с целыми значениями (с дубликатом), второй массив n на 1 'gid', содержащий уникальные целочисленные значения из 'ppi' и третий массив «X0» с 0 и 1 значениями. Ниже приведен код MATLAB, который я хочу реализовать в Java-Экспоненциальная функция, дающая разные результаты в Java и Matlab
MATLAB версию-
gamma0 = 10; gamma1 = 10; beta = 20;
delta = 0;weight=10;
for i=1:length(gid)
gconn = union(ppi(find(ppi(:,1)==gid(i)),2), ppi(find(ppi(:,2)==gid(i)),1));
[a,b] = intersect(gid, gconn);
u1 = (X0(i)*weight + sum(X0(b)==1))/(weight+length(gconn));
u0 = ((1-X0(i))*weight + sum(X0(b)==0))/(weight+length(gconn));
mrfpdf(i) = exp((1-X0(i))*(gamma0-beta*u1)+X0(i)*(gamma1-beta*u0))/(exp(gamma0- beta*u1)+exp(gamma1-beta*u0));
end
Java код для той же реализации
double[] mrfpdf = new double[X0.length];
Integer[] X = new Integer[X0.length];
// find out unique class labels in ppi
Set<Double> gid_set = new LinkedHashSet<Double>();
for (int i = 0; i < ppi.length; i++) {
for (int j = 0; j < ppi[0].length; j++) {
gid_set.add(ppi[i][j]);
}
}
System.out.println("unique gid_set: " + gid_set);
Double[] gid = gid_set.toArray(new Double[gid_set.size()]);
Arrays.sort(gid);
int gamma0 = 10, gamma1 = 10, beta = 20, delta = 0;weight=10;
for (int i = 0; i < gid.length; i++) {
Set<Double> ppi1 = new HashSet<Double>();
Set<Double> ppi2 = new HashSet<Double>();
for (int k = 0; k < ppi.length; k++) {
if (ppi[k][0] == gid[i])
ppi1.add(ppi[k][1]);
if (ppi[k][1] == gid[i])
ppi2.add(ppi[k][0]);
}
Set<Double> gconn = new HashSet<Double>(ppi1);
gconn.addAll(ppi2);
Set<Double> a = new HashSet<Double>(gid_set);
a.retainAll(gconn);
java.util.List<Double> buff = new LinkedList<Double>(gid_set);
int[] b = new int[a.size()];// get index of intersection elements in this array
int x = 0;
for (double j : a) {
for (int c = 0; c < gid.length; c++) {
if (j == gid[c]) {
b[x] = c;
x++;
}
}
}
java.util.List<Double> X0_list = new ArrayList<Double>();
for (double z : X0) {
X0_list.add(z);
}
int sum1 = 0;
for (int y : b) {
if (X0_list.get(y) == 1)
sum1++;
}
int sum2 = 0;
for (int y : b) {
if (X0_list.get(y) == 0)
sum2++;
}
u1 = (int) ((X0[i] * weight + sum1)/(weight + gconn.size()));
u0 = (int) (((1 - X0[i]) * weight + sum2)/(weight + gconn.size()));
mrfpdf[i] = (Math.exp((1 - X0[i]) * (gamma0 - beta * u1) + X0[i]* (gamma1 - beta * u0))/(Math.exp(gamma0 - beta * u1) + Math.exp(gamma1 - beta * u0)));
}// end for(int i : gid)
}
Может кто-то пожалуйста, скажите мне, почему значения в mrfpdf массив отличается в Java и Matlab-коде? Результаты кода Java в 0,99 и 0,5 значения только тогда как Matlab не имеют 0,5 значения
Благодаря
вы не возражаете форматирование? Никто не прочитает это! – Fildor
И я уверен, что вы знаете, что [каждый компьютерный ученый должен знать о арифметике с плавающей точкой] (http://docs.oracle.com/cd/E19957-01/806-3568/ncg_goldberg.html)? – Fildor