24小时热门版块排行榜

>论坛更新日志 (3396)
>虫友互识 (477)
>导师招生 (171)
>休闲灌水 (166)
>文献求助 (162)
>论文投稿 (114)
>考博 (110)
>博后之家 (83)
>论文道贺祈福 (65)
>招聘信息布告栏 (62)
>基金申请 (52)
>硕博家园 (49)
>教师之家 (47)
>找工作 (45)
>考研 (42)
>公派出国 (20)

返回列表

H_LG

银虫 (正式写手)

应助: 2 (幼儿园)
金币: 157.8
散金: 228
帖子: 334
在线: 82.3小时
虫号: 579678
注册: 2008-07-11
性别: GG
专业: 计算数学与科学工程计算

[求助] 求高手将一个M函数转换成mex文件或mex语法编写的c代码

如题，将以下程序进行转换
因为这个程序在我的matlab程序中多次调用，很耗时。
求高手解答或指点。
谢过！

function [A]=OMPerr(D,X,errorGoal)%#codegen
% ========================================================
% Sparse coding of a group of signals based on a given dictionary and specified representation
% error to get.
% input arguments: D - the dictionary
% X - the signals to represent
% errorGoal - the maximal allowed representation error
% output arguments: A - sparse coefficient matrix.
% ========================================================
[n,P]=size(X);
[n,K]=size(D);
E2 = errorGoal^2*n;
maxNumCoef = n/2;
A = zeros(size(D,2),size(X,2));
% h=waitbar(0,'OMP on each example ...');
for k=1:1

,
%    waitbar(k/P);
a=[];
x=X(:,k);
residual=x;
indx = [];
a = [];
currResNorm2 = sum(residual.^2);
j = 0;
while currResNorm2>E2 && j < maxNumCoef,
      j = j+1;
      proj=D'*residual;
      pos=find(abs(proj)==max(abs(proj)));
      pos=pos(1);
      indx(j)=pos;
      a=pinv(D(:,indx(1:j)))*x;
      residual=x-D(:,indx(1:j))*a;
      currResNorm2 = sum(residual.^2);
end;
if (~isempty(indx))
      A(indx,k)=a;
end
end;
% close(h);
return;[ 来自科研家族控制&图像家族 ]

回复此楼