compare_marginals.m.svn-base

来自「Probabilistic graphical models in matlab」· SVN-BASE 代码 · 共 58 行

function [error] = compare_marginals(filename_1, filename_2)

fid = fopen(filename_1, 'r');

a = fscanf(fid, '<marginals> ');

fscanf(fid, '( ');

t = fscanf(fid, '%d');

while size(t) > 0

fscanf(fid, ' ) (');

temp =  fscanf(fid, ' %g');

marginals_1(t) = { temp' };

fscanf(fid, ') ( ');

t = fscanf(fid, '%d');

end

fid = fopen(filename_2, 'r');

a = fscanf(fid, '<marginals> ');

fscanf(fid, '( ');

t = fscanf(fid, '%d');

while size(t) > 0

fscanf(fid, ' ) (');

temp =  fscanf(fid, ' %g');

marginals_2(t) = { temp' };

fscanf(fid, ') ( ');

t = fscanf(fid, '%d');

end

error = 0;

for i =1:size(marginals_1,2)
    
    for j =1:size(marginals_1{i},2)
    
   error = error + (marginals_1{i}(j) - marginals_2{i}(j) )^2;
    
    end
end