2013-10-15 60 views
0

我有一個n乘2的數組'ppi',它具有整數值(帶有重複),第二個數組n'1'gid'包含來自'ppi'的唯一整數值​​,數組'X0'有0和1個值。下面是我想在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代碼相同的實施MATLAB代碼

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代碼中有所不同?在0.99 Java代碼的結果和0.5的值只,而Matlab的沒有任何值0.5

感謝

+1

你不介意格式化?沒有人會閱讀! – Fildor

+0

我相信你知道[每個計算機科學家應該知道什麼是浮點運算](http://docs.oracle.com/cd/E19957-01/806-3568/ncg_goldberg.html)? – Fildor

回答

1

您的MATLAB代碼使用雙打,而Java的代碼使用整數(例如U0和U1)

+1

謝謝!我改變了u1和u2,讓它變成了雙倍效果! – user2006189

相關問題