| 查看: 662 | 回复: 0 | ||
H_LG银虫 (正式写手)
|
[求助]
求高手将一个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;[ 来自科研家族 控制&图像家族 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : OMPerr.m
2011-12-23 17:13:13, 1.17 K
» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
听说今天filecode变了
已经有34人回复
-
基金中了
已经有6人回复
-
filecode
已经有14人回复
-
固定端突然变了，今天
已经有7人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
找到一些相关的精华帖子，希望有用哦~
【JPCC】-太阳能光转换：TiO2@Graphene核壳结构 VS TiO2/Graphene负载结构
已经有118人回复
- 帮忙看一下几句话的英文翻译。谢谢！水处理类。帮忙改一下语法和结构。BB
已经有5人回复
- 有关量化计算中稀有气体He Ne Ar Kr Xe的赝势基组函数
已经有3人回复
- 请问如果知道了导热系数，如何转换为热阻，
已经有5人回复
- 数据是.CSV格式的，我想把它转换成红外光谱图，请问该怎么转换？
已经有4人回复
- 【求助】帮忙读一个mex程序，100金币相送
已经有5人回复
- 【求助】急求高手帮忙！只有cif如何晶体画图（绝对没造假！）
已经有10人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
100