find_db.m

cootes提出的ASM算法的matlab实现。在cootes用于课堂使用的源代码的基础上改写。改正了原来代码中的一些错误

function db=find_db(dx,P);
%function db=find_db(dx,P);


% find db by transforming dx into the model parameter space
%        db  = P\dx
% dx=P*db

db=P\dx;


% 以下部分在调试C程序时使用
fid=fopen('find_db_record.txt','a');
if fid==0
    msgbox('fopen(''find_db_record.txt'',''a'') error');
    return;
end

fprintf(fid,'\r\n\r\n----------------------------------------\r\n\r\n');

fprintf(fid,'P=\r\n');
for i=1:1:size(P,1)
    for j=1:1:size(P,2)
        fprintf(fid,'%f',P(i,j));
        if j==size(P,2)
            fprintf(fid,'\r\n');
        else
            fprintf(fid,'\t');
        end
    end
end
fprintf(fid,'\r\ndx=\r\n');
for i=1:1:length(dx)
    fprintf(fid,'%f\r\n',dx(i));
end
fprintf(fid,'\r\ndb=\r\n');
for i=1:1:length(db)
    fprintf(fid,'%f\r\n',db(i));
end

fclose(fid);