亚洲人成77777在线观看网,亚洲国产日韩欧美在线99,国产欧美在线播放

How to read *.data in Matlab and Python

杰哥 — Fri, 12 Oct 2018 20:33:00 GMT

See P2 of "HW2_545_2018_"

In Matlab:

z = dlmread('spambase.data',',');

In Python:

import numpy as np

z = np.genfromtxt('spambase.data', dtype=float, delimiter=',')

杰哥 2018-10-13 04:33 发表评论

matlab legend字号改变

杰哥 — Mon, 18 Dec 2017 20:50:00 GMT

hleg1 = legend('SBQ','DBQ','IBQ',4);
set(hleg1,'Interpreter','none');
set(hleg1,'FontSize',15);

set(gca,'linewidth',2);
set(gca,'FontName','Times New Roman','FontSize',12);
set(gca,'xtick',[16 32 64 128]);

This is with Yuan Cao's help.

杰哥 2017-12-19 04:50 发表评论

[转蝲]matlab中text 函数在显�C�字�W�串时的使用�Ҏ��

杰哥 — Sun, 17 Dec 2017 04:53:00 GMT

原文地址�Q?/strong>matlab中text 函数在显�C�字�W�串时的使用�Ҏ��作者：gypsy

matlab中text 函数在显�C�字�W�串时的使用�Ҏ��

功能在当前��u中创建text对象。函数text是创建text囑�Ş句柄的低�U�函数。可用该函数在图形中指定的位�|�上昄��字符丌Ӏ?/p>

用法 text(x,y,'string')在图形中指定的位�|?x,y)上显�C�字�W�串string

text(x,y,z,'string') 在三�l�图形空间中的指定位�|?x,y,z)上显�C�字�W�串string

text(x,y,z,’string’.'PropertyName',PropertyValue…) 对引号中的文字string定位于用坐标轴指定的位置�Q�且�Ҏ��定的属性进行设�|�。表7-6�l�出文字属性名、含义及属性倹{�?/p>

�?-6

属性名	属性说�?/p>	属性�?/p>
定义字符�?/p>
Editing	能否�Ҏ��字进行编�?/p>	有效��\|��on、off �~�省��\|��off
Interpretation	TeX字符是否可用	有效��\|��tex、none �~�省��\|��tex
String	字符�?包括TeX字符�?	有效��\|��可见字符�?/p>
攄��字符�?/p>
Extent	text对象的范��_��位置与大��）	有效��\|��[left, bottom, width, height]
HorizontalAlignment	文字水��^方向的对齐方�?/p>	有效��\|��left(文本外框左边寚w��Q�缺省对齐方�?、center(文本外框中间寚w��)、right(文本外框双��寚w��) �~�省��\|��left
Position	文字范围的位�\|?/p>	有效��\|��[x,y,z]直角坐标�p?/p> �~�省��\|��[]�Q�空矩阵�Q?/p>
Rotation	文字对象的方位角�?/p>	有效��\|��标量�Q�单位�ؓ度） �~�省��\|��0
Units	文字范围与位�\|�的单位	有效��\|��pixels (屏幕上的像素�?、normalized (把屏�q�看成一个长、宽�?的矩�?、inches(英寸)、centimeters(厘米)、points (图象�?、data �~�省��\|��data
VerticalAlignment	文字垂直方向的对齐方�?/p>	有效��\|��top (文本外框��上寚w��)、cap(文本字符��上寚w��)、middle(文本外框中间寚w��)、baseline(文本字符底线�?、bottom(文本外框底线寚w��) �~�省��\|��middle
指定文字字体
FontAngle	讄��斜体文字模式	有效��\|��normal(正常字体)、italic(斜体�?、oblique(斜角�? �~�省��\|��normal
FontName	讄��文字字体名称	有效��\|��用户�pȝ��支持的字体名或者字�W�串FixedWidth�?/p> �~�省��gؓ Helvetica
FontSize	文字字体大小	有效��\|��l�合字体单位的数�?/p> �~�省��gؓ�Q?0 points
FontUnits	讄��属性FontSize的单�?/p>	有效��\|��points (1�?1/72英寸)、normalized(把父对象坐标轴作��Z��单位长的一个整体；当改变坐标��u的尺寸时�Q�系�l�会自动改变字体的大��?、inches (英寸)、Centimeters(厘米)、Pixels(像素) �~�省��\|��points
FontWeight	讄��文字字体的粗�l?/p>	有效��\|��light(�l�字�?、normal(正常字体)、demi(黑体�?、Bold(黑体�? �~�省��\|��normal
控制文字外观
Clipping	讄��坐标轴中矩�Ş的剪辑模�?/p>	有效��\|��on、off on�Q�当文本��出坐标轴的矩�Ş�Ӟ��出的部分不昄��Q?/p> off�Q�当文本��出坐标轴的矩�Ş�Ӟ��出的部分显�C��?/p> �~�省��\|��off
EraseMode	讄��昄��与擦除文字的模式。这些模式对生成动画�p�d��与改�q�文字的昄��效果很有好处�?/p>	有效��\|��normal、none�?xor�?background �~�省��\|��normal
SelectionHighlight	讄��选中文字是否�H�出昄��	有效��\|��on、off �~�省��\|��on
Visible	讄��文字是否可见	有效��\|��on、off �~�省��\|��on
Color	讄��文字颜色	有效的颜色��\|��ColorSpec
控制�Ҏ��字对象的讉K��
HandleVisibility	讄��文字对象句柄对其他函数是否可�?/p>	有效��\|��on、callback、off �~�省��\|��on
HitTest	讄��文字对象能否成�ؓ当前对象�Q�见囑�ŞCurrentObject属性）	有效��\|��on、off �~�省��\|��on
文字对象的一般信�?/p>
Children	文字对象的子对象�Q�文字对象没有子对象�Q?/p>	有效��\|��[]�Q�即�I�矩阵）
Parent	文字对象的父对象�Q�通常为axes对象�Q?/p>	有效��\|��axes的句�?/p>
Seleted	讄��文字是否昄��?#8220;选中”状�?/p>	有效��\|��on、off �~�省��\|��off
Tag	讄��用户指定的标�{?/p>	有效��\|��M��字符�?/p> �~�省��\|��’’�Q�即�I�字�W�串�Q?/p>
Type	讄��囑�Ş对象的类型（只读�c�d��Q?/p>	有效��\|��字符�?#8217;text’
UserData	讄��用户指定数据	有效��\|��M��矩阵 �~�省��\|��[]�Q�即�I�矩阵）
控制回调例行执行�E�序
BusyAction	讄��如何处理�Ҏ��字回调过�E�中断的句柄	有效��\|��cancel、queue �~�省��\|��queue
ButtonDownFcn	讄��当鼠标在文字上单��L��Q�程序做出的反应�Q�即执行回调�E�序�Q?/p>	有效��\|��字符�?/p> �~�省��\|��' '�Q�空字符�Ԍ��
CreateFcn	讄��当文字被创徏�Ӟ��E�序做出的反应（��x��行的回调�E�序�Q?/p>	有效��\|��字符�?/p> �~�省��\|��' '�Q�空字符�Ԍ��
DeleteFcn	讄��当文字被删除�Q�通过关闭或删除操作）�Ӟ��E�序做出的反应（��x��行的回调�E�序�Q?/p>	有效��\|��字符�?/p> �~�省��\|��' '�Q�空字符�Ԍ��
Interruptible	讄��回调�q�程是否可中�?/p>	有效��\|��on、off �~�省��\|��on�Q�能中断�Q?/p>
UIContextMenu	讄��与文字相关的菜单��?/p>	有效��\|��用户相关菜单句柄

MATLAB资源�|�：http://www.ymlib.net/article/sort010/info-276.html
http://blog.sciencenet.cn/blog-711035-566055.html

杰哥 2017-12-17 12:53 发表评论

matlab ColorSpec (Color Specification)

杰哥 — Sat, 09 Dec 2017 18:27:00 GMT
https://www.mathworks.com/help/matlab/ref/colorspec.html

RGB Triplet
Short Name
Long Name
[1 1 0]
y
yellow
[1 0 1]
m
magenta
[0 1 1]
c
cyan
[1 0 0]
r
red
[0 1 0]
g
green
[0 0 1]
b
blue
[1 1 1]
w
white
[0 0 0]
k
black

杰哥 2017-12-10 02:27 发表评论

Professor Deng Cai's code, very good code

杰哥 — Sun, 26 Nov 2017 05:10:00 GMT
20171125 rutgers email: Please check https://github.com/dengcai78/MatlabFunc for the most up-to-date codes.

杰哥 2017-11-26 13:10 发表评论

[zz]MATLAB如何实现十进制数与二�q�制数的转换

杰哥 — Wed, 08 Nov 2017 04:49:00 GMT
在MATLAB 中，实现十进制与二进制�{化的函数有dec2bin()�Q�bin2dec()�?br />
Reference
http://jingyan.baidu.com/article/59703552c5699a8fc00740d2.html

杰哥 2017-11-08 12:49 发表评论

matlab dir

杰哥 — Mon, 25 Sep 2017 16:28:00 GMT

Matlab帮助: list folder contents.
针对根目录，列出所有子文�g夹，包括隐藏文�g。注意针寚w��根目录，�q�有多出两个文�g�?."�?.."�Q�分别对应本目录和上一层目录。故TrainDataLoading.m的第八行和第九行如下�Q?br />
templist = dir(session_dir_name);
folderlist = templist(3:end); % delete the first two null directory
见本机dir_.m和TrainDataLoading.m�?br />用windows菜单�Q�开始，所有程序，附�g�Q�命令提�C�符�Q�输入dir��׃��列出所有目录，会看到有两个文�g�?."�?.."。cd C:\�Q�定位到根目录，再dir�Q�没有两个文件夹"."�?.."�?br />This is with Nannan and Zhongze's help.

杰哥 2017-09-26 00:28 发表评论

杰哥 — Fri, 29 Jul 2016 07:13:00 GMT

首先�Q�P文�g��Z��保护知识产权设计的一�U�加密文�Ӟ��是不能查看的�?/div>
�q�行比较��单，和调用m文�g�Ҏ��一栗��?/div>
下面是一些解释�?/div>
P文�g是对应M文�g的一�U�预解析版本�Q�preparsed version�Q�。因为当你第一�ơ执行M文�g�Ӟ��Matlab需要将其解析（parse�Q�一�ơ（�W�一�ơ执行后的已解析内容会放入内存作�W�二�ơ执行时使用�Q�即�W�二�ơ执行时无需再解析）�Q�这无�Ş中增加了执行旉��。所以我们就预先作解释，那么以后再��用该M文�g�Ӟ��便会直接执行对应的已解析版本�Q�即P文�g。但又因为Matlab的解析速度非常快，一般不用自�׃��预解析。只有当一些程序要调用到非常多的M文�g�Ӟ��如GUI应用�E�序�Ӟ��才会作预解析�Q�以增加以后的调用速度�?nbsp;
P文�g可以用来作保密代码之用，如果你给别�h一个M文�g�Q�别人可以打开来看��C��所有的代码和算法�?span style="color: red;">如果你的代码不想被别人看刎ͼ�那可以给他P文�g�?/span>
如何使用.p文�g
1.首先�Q�将当前工作目录切换�?p文�g所在的目录�Q�然后就可以在左侧的工作�I�间�H�口看见该目录所包含的所有文件了�Q�当然包�?p文�g啦~
2.�׃��.p文�g�?m文�g的预解析版本�Q�不能查看、编辑）�Q�显�?p文�g�?m文�g一样都是一个函��C��Q�只不过此时该函数是一个黑盒子�|�了。�D例如test.p文�g便意味着存在一个已�l�编好的函数test�Q�如何实现我们不知道�Q�，��用这个函数test�Q�我们需要知道它的参数�Ş�?此时在命令窗口下键入help test�Q�即可查看到test函数的参�?然后我们��可以调用啦
如何自己制作.p文�g
在Matlab中，我们�l�常把语句或函数写在一个M文�g里，便于模块化和代码重用�Q�比如，写一个test.m的文�Ӟ��卌��定义一个test函数�Q�：
% ———————-
% test.m
a = 1;
b = 2;
% ———————-
�q�个M语句文�g�Q�或�Q?/div>
% ———————-
% test2.m
function a = test2(b)
a = b * 2;
% ———————-
�q�个M函数�?/div>
我们在Command Window里输入：
>> test
��x��{�于输入�?/div>
>> a = 1;
>> b = 2;
而输入：
>> a = test2(2)
a =
4
便是执行了一个函数�?/div>
以上的内�Ҏ��惛_��部䆾人都知道是怎么一回事了，以下说一说P代码文�g�?/div>
如果在Command Window里输入：
>> pcode test
便会在相应test.m的文件夹里��生了test.p。如果这时在Command Window里输入：
>> test
其实是运行了test.p�Q�而不是test.m�?/span>
P文�g是对应M文�g的一�U�预解析版本�Q�preparsed version�Q�。因为当你第一�ơ执行M文�g�Ӟ��Matlab需要将其解析（parse�Q�一�ơ（�W�一�ơ执行后的已解析内容会放入内存作�W�二�ơ执行时使用�Q�即�W�二�ơ执行时无需再解析）�Q�这无�Ş中增加了执行旉��。所以我们就预先作解释，那么以后再��用该M文�g�Ӟ��便会直接执行对应的已解析版本�Q�即P文�g。但又因为Matlab的解析速度非常快，一般不用自�׃��预解析。只有当一些程序要调用到非常多的M文�g�Ӟ��如GUI应用�E�序�Ӟ��才会作预解析�Q�以增加以后的调用速度�?/div>
如Matlab的当前目录（Current Directory�Q�有test.m文�g�Q�作预解析后�Q�又有test.p文�g。因为P文�g的调用优先��比M文�g要高�Q�所以当你调用test�Ӟ��会作优先选择而调用了test.p�?/div>
可以修改test.m的代码�ؓ�Q?/div>
% ———————-
% test.m
a = 3;
b = 4;
% ———————-
再在Command Window里调用test�Q?/div>
>> test
Warning: P-file C:\Program Files\MATLAB71\work\test.p is older than M-file C:\Program Files\MATLAB71\work\test.m.
C:\Program Files\MATLAB71\work\test.p may be obsolete and may need to be regenerated.
Type "help pcode" for information about generating P-files from M-files.
会出��C��个Warning�Q�警告你所调用的P文�g比同名的M文�g要旧�Q�即表示M文�g已被修改了。尽��如此，调用的还是旧的P文�g�Q�即得出a = 1�?b = 2�?/div>
P文�g可以用来作保密代码之用，如果你给别�h一个M文�g�Q�别人可以打开来看��C��所有的代码和算法。如果你的代码不惌��别�h看到�Q�那可以�l�他P文�g�?/div>
pcode函数也可以应用在M函数文�g�?/div>
参考：
http://zmryyj.blog.163.com/blog/static/98450751201132081451490/
http://frenselx.spaces.live.com/blog/cns!3D7342371BEE8602!290.entry
阅读记录: read twice, understand completely

杰哥 2016-07-29 15:13 发表评论

The matlab code of Locality sensitive hashing (LSH)

杰哥 — Fri, 27 Nov 2015 14:24:00 GMT
%see the matlab code test.m of CCA-ITQ. It is very simple. Only one sentence using random projection
XX = XX * randn(size(XX,2),bit);
Y = zeros(size(XX));
Y(XX>=0)=1;
Y = compactbit(Y);

杰哥 2015-11-27 22:24 发表评论

matlab函数meshgrid

杰哥 — Fri, 11 Sep 2015 04:42:00 GMT
问题�Q?/span>
使用matlab怎么能够产生�q�样的一�l�数�l�呢�Q?br />x = 1,2,3,4
y = 5,6

产生一�l�点�Q?/p>
(1,5) (1,6)
(2,5) (2,6)
(3,5) (3,6)
(4,5) (4,6)

�{�案�Q?/span>

x = [1,2,3,4];
y = [5,6];
[X,Y] = meshgrid(x,y);
[X(:), Y(:)]

杰哥 2015-09-11 12:42 发表评论

matlab函数fminunc:求无�U�束优化的最优解

杰哥 — Thu, 13 Aug 2015 04:59:00 GMT
例：K-means Hashing提供的代码optimize_one_iter函数内部用到fminunc�Q�代价函数是eval_one_center函数�Q�看eval_one_center函数的注释是this function computes the cost of Eqn. (8) given a center c_j

杰哥 2015-08-13 12:59 发表评论

matlab加�? bsxfun

杰哥 — Wed, 05 Aug 2015 03:18:00 GMT
%��训�l�数据和��试数据都去中心�?/span>
X = traindata;
label = traingnd;
m = mean(X);
X_zm = bsxfun(@minus, X, m);
traindata_zm = bsxfun(@minus, traindata, m);
testdata_zm = bsxfun(@minus, testdata, m);

matlab函数 bsxfun��谈�Q��{载）
http://blog.sina.com.cn/s/blog_9e67285801010ttn.html
�|�上关于bsxfun的东西不多，今天需要看��C��个，�׃��原博文插入的囄��昄��不出来，于是�W�者大发善心进行了contrl+V 以及alt+ctrl+A的操作，供大家交��学习�?/span>

bsxfun是一�?/span>matlab自版�?/span>R2007a来就提供的一个函敎ͼ�作用�?/span>”applies an element-by-element binary operation to arrays a and b, with singleton expansion enabled.”

举个例子。假设我们有一列向量和一行向量�?/span>

a = randn(3,1), b = randn(1,3) a = -0.2453 -0.2766 -0.1913 b = 0.6062 0.5655 0.9057

我们可以很简单的使用matlab的外�?/span>c=a*b来得刎ͼ�如图 bsxfun��谈�Q��{载）" o:button="t" target='"_blank"' o:spid="_x0000_i1025">bsxfun��谈�Q��{载）" src="file:///C:\Users\jie\AppData\Local\Temp\msohtmlclip1\01\clip_image001.jpg">
但如果我们想�?/span>”外加”呢？也就是说把上式求解过�E�中的乘��h��做加��P��

�q�时我们可以�?/span>c=bsxfun(@plus,a,b)来实现�?/span>

bsxfun的执行是�q�样的，如果a�?/span>b的大��相同，那么c=a+b. 但如果有某维不同�Q�且a�?/span>b必须有一个在�q�一�l�的�l�数�?/span>1, 那么bsxfun��将��的�q�个虚拟的复制一些来使与多的�l�数一栗��在我们�q�里�Q?/span>b的第一�l�只�?/span>1�Q�只一行）�Q�所�?/span>bsxfun��?/span>b复制3�ơ�Ş成一�?/span>3×3的矩阵，同样也将a复制�?/span>3×3的矩��c��这个等价于c=repmat(a,1,3)+repmat(b,3,1)。这�?/span>

repmat(a,1,3) ans = -0.2453 -0.2453 -0.2453 -0.2766 -0.2766 -0.2766 -0.1913 -0.1913 -0.1913

repmat是显式的复制�Q�当然带来内存的消耗。�?/span>bsxfun是虚拟的复制�Q�实际上通过for来实玎ͼ��{�效�?/span>for(i=1:3),for(j=1:3),c(i,j)=a(i)+b(j);end,end。但bsxfun不会有��?/span>matlab�?/span>for所带来额外旉��。实际验证下�q�三�U�方�?/span>

>> c = bsxfun(@plus,a,b) c = 0.3609 0.3202 0.6604 0.3296 0.2889 0.6291 0.4149 0.3742 0.7144 >> c = repmat(a,1,3)+repmat(b,3,1) c = 0.3609 0.3202 0.6604 0.3296 0.2889 0.6291 0.4149 0.3742 0.7144 >> for(i=1:3),for(j=1:3),c(i,j)=a(i)+b(j);end,end,c c = 0.3609 0.3202 0.6604 0.3296 0.2889 0.6291 0.4149 0.3742 0.7144

从计��时间上来说前两�U�实现差不多�Q�远高于for的实现。但如果数据很大�Q�第二种实现可能会有内存上的问题�?/span>所以bsxfun最好�?br />

下面看一个更为实际的情况。假设我们有数据A和B, 每行是一个样本，每列是一个特征。我们要计算高斯核，既：

�q�里@plus是加法的函数数柄�Q�相应的有减法@minus, 乘法@times, 左右除等�Q�具体可�?doc bsxfun.

k(||x-xc||)=exp{- ||x-xc||^2/(2*σ)^2) } 其中xc为核函数中心,σ为函数的宽度参数 , 控制了函数的径向作用范围�?/span>

当然可以用双重for实现�Q�如果第一直觉是用三重for的话…�Q��?br />

K1 = zeros(size(A,1),size(B,1)); for i = 1 : size(A,1) for j = 1 : size(B,1) K1(i,j) = exp(-sum((A(i,:)-B(j,:)).^2)/beta); end end

使用2,000×1,000大小的A和B, �q�行旉��?8�U��?/span>
考虑下面向量化后的版本：

sA = (sum(A.^2, 2)); sB = (sum(B.^2, 2)); K2 = exp(bsxfun(@minus,bsxfun(@minus,2*A*B', sA), sB')/beta);

使用同样数据�Q�运行时间仅0.85�U�，加速超�q?00倍�?/span>
如要判断两者结果是不是一��P��可以如下

assert(all(all(abs(K1-K2)<1e-12)))

杰哥 2015-08-05 11:18 发表评论

matlab函数cat

杰哥 — Thu, 23 Jul 2015 00:36:00 GMT
问题�Q�两个图像矩阵A=M*N*K1�Q�B=M*N*K2�Q�第三维K1�Q�K2表示数量�Q�matlab怎么快速实��C��们在�W�三�l�数量上的叠加�?br />�{�：cat(3,A,B)

杰哥 2015-07-23 08:36 发表评论

How to solve AX + XB = C for X using matlab?

杰哥 — Mon, 06 Jul 2015 07:28:00 GMT
X = sylvester(A,B,C)
http://cn.mathworks.com/help/matlab/ref/sylvester.html

杰哥 2015-07-06 15:28 发表评论

f-measure matlab

杰哥 — Tue, 30 Jun 2015 13:07:00 GMT
supervised discrete hashing的evaluate_macro代码理解�Q�在�ȝ��序中调用�Ӟ��evaluate_macro的前一个输入参数是cateTrainTest�Q�是训练��h��?��试��h��数的矩阵。evaluate_macro函数体中�Q�对每个��试��h��Q�retrieved_relevant_num代表TP�Q�relevant_num代表TP+FN�Q�retrieved_num代表TP+FP。所有样本分别的�q�_��precision��x��最�l�的输出precision�Q�所有样本分别的�q�_��recall��x��最�l�的输出recall�?br />cat_apcal函数是计��MAP的，代码意义很好理解�Q�比如按距离排序�Q�前七个只有1�?�?�?是和query相同的类�Q�则MAP=(1+2/3+3/5+4/7)/4 (This is with discussing with Shu Zhang). 但Deep hashing for compact (CVPR 2015) mean average
precision (mAP): which computes the area under the precision-recall curve.
matlab曲线下如何求面积?
如果知道函数表达式的话，调用quad函数��可以了。如果不知道函数表达式只知道�q�一�p�d��L��点，x和y�Q�则trapz(x,y)卛_��
�?�Q�本机Hashing\code\ITQ\delete\test_cifar_PCA_ITQ_V01.m�Q�对��d��的precison和recall曲线�Q�计��曲�U�下的面�U�，用trapz(recall,precision)卛_��
Shu Zhang说不知以上两�U�是否等��P��如果按照Deep hashing for compact (CVPR 2015)��面�U�，则直接就能看��Z��同的�Ҏ��谁的MAP大，因�ؓ谁的曲线在上方就谁大

�w�䆾证照片和真�h的比�Ҏ��两个指标
误识率：他明明不是张三，但你错误地把他判断成了张三。误识率1%�Q�意味着�?00个�h来冒充其他�h�Q�有1个�h会冒充成功。应该就是FAR
识别率：张三�?0%左右的概率被�pȝ��说你不是张三�Q�识别率90%。应该就�?-FRR
二代证比真�h�Q�识别率90% (误识�?%)
即��?%的认假率的情况下�Q�识别率目前可能大概做到90%左右。当然这个前提条件是��_��我们用的是二代证卡内的那张低质量�Q�压�~�得非常狠，照片的分辨率不够高的那张��照片。如果我们用的是一张清晰的�Q�近期的证�g照，那么�q�个�l�论可以做到大概�Q�误识率可以做到万分之一的情况下�Q�识别率可以做到90%以上的正��率�?/div>
Reference�Q?/div>
造就Talk | �׃��光：天眼�pȝ��l�将开启，你准备好了吗�Q�http://mp.weixin.qq.com/s?__biz=MzAwNTcyNDU5MQ==&mid=509662068&idx=1&sn=f14840fb7927497ad0d49e1d18eb6b5a&scene=1&srcid=0529xR8z0qC240GLkQZuUAG9&from=groupmessage&isappinstalled=0#wechat_redirect 已看两次

http://www.mathworks.com/matlabcentral/fileexchange/37758-performance-measures-for-classification/content/Evaluate.m
function EVAL = Evaluate(ACTUAL,PREDICTED) % This fucntion evaluates the performance of a classification model by % calculating the common performance measures: Accuracy, Sensitivity, % Specificity, Precision, Recall, F-Measure, G-mean. % Input: ACTUAL = Column matrix with actual class labels of the training % examples % PREDICTED = Column matrix with predicted class labels by the % classification model % Output: EVAL = Row matrix with all the performance measures idx = (ACTUAL()==1); p = length(ACTUAL(idx)); n = length(ACTUAL(~idx)); N = p+n; tp = sum(ACTUAL(idx)==PREDICTED(idx)); tn = sum(ACTUAL(~idx)==PREDICTED(~idx)); fp = n-tn; fn = p-tp; tp_rate = tp/p; tn_rate = tn/n; accuracy = (tp+tn)/N; sensitivity = tp_rate; specificity = tn_rate; precision = tp/(tp+fp); recall = sensitivity; f_measure = 2*((precision*recall)/(precision + recall)); gmean = sqrt(tp_rate*tn_rate); EVAL = [accuracy sensitivity specificity precision recall f_measure gmean];

杰哥 2015-06-30 21:07 发表评论

The matlab code of canonical correlation analysis (CCA)

杰哥 — Tue, 02 Jun 2015 14:10:00 GMT

%Reference: the matlab code of "CVPR 2011 Iterative Quantization". I have checked the correctness of the corresponding derivation and matlab code. No need to check again. Just use it

function [Wx, r] = cca(X,Y,reg)

%
% X is input data for the 1st view (image), n*d, n images, d dim
% Y is input data for the 2nd view (Tag), n*D, D dim
% reg is regularization parameter, usually set by validation
% in our work, reg = 0.0001 which works well
%
% Wx is the embedding function for image
% r is the eigenvalue
%
%
% to run the code, you need visual data X and tag data Y
%
% bit = 32, 64, 128 ...
% [eigenvector,r] = cca(X, Y, 0.0001); % this computes CCA projections
% eigenvector = eigenvector(:,1:bit)*diag(r(1:bit)); % this performs a scaling using eigenvalues
% E = X*eigenvector; % final projection to obtain embedding E
%

z = [X, Y];
C = cov(z);
sx = size(X,2);
sy = size(Y,2);
Cxx = C(1:sx, 1:sx) + reg*eye(sx);
Cxy = C(1:sx, sx+1:sx+sy);
Cyx = Cxy';
Cyy = C(sx+1:sx+sy,sx+1:sx+sy) + reg*eye(sy);

Rx = chol(Cxx);
invRx = inv(Rx);
Z = invRx'*Cxy*(Cyy\Cyx)*invRx;
Z = 0.5*(Z' + Z);

[Wx,r] = eig(Z);   % basis in h (X)
r = sqrt(real(r)); % as the original r we get is lamda^2
Wx = invRx * Wx;   % actual Wx values
r = diag(r);

[r index] = sort(r,'descend');
Wx = Wx(:,index);

Other references:
Information Fusion 2017 Multi-view learning overview的公�?�?
Section 2.1 of "NCA 2013 A survey of multi-view machine learning"

杰哥 2015-06-02 22:10 发表评论

[zz]matlab囑փ�滤�L

杰哥 — Mon, 11 May 2015 12:33:00 GMT
clc;
clear;
I=imread('a.jpg');             %��d��囑փ�
K1=medfilt2(I,[10 10]);          %中值��o�?br />h1=ones(5,5)/25;
K2=imfilter(K1,h1);       %均值��o�?br />figure, imshow(K2), title('滤�L后的囑փ�')

杰哥 2015-05-11 20:33 发表评论

The matlab code to normalize each feature of training data and test data to unit mean and zero variance

杰哥 — Fri, 05 Dec 2014 04:32:00 GMT
    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    fea_Train = fea(:,trainIdx);
    gnd_Train = gnd(trainIdx);
    fea_Test = fea(:,testIdx);
    gnd_Test = gnd(testIdx);
    train_size = length(gnd_Train);
    test_size = length(gnd_Test);

    %%%%%%%%%%%%%%%%%%%%%%%Normalization begins%%%%%%%%%%%%%%%%%%%%%%%
    train_data_std = std(fea_Train');
    train_data_mean = mean(fea_Train');
    train_data_std_matrix = ones(train_size,dim)*diag(train_data_std);
    train_data_std_matrix(train_data_std_matrix==0)=1;
    train_data_mean_matrix = ones(train_size,dim)*diag(train_data_mean);
    train_scale_data = (fea_Train' - train_data_mean_matrix)./train_data_std_matrix;
    clear train_data_mean_matrix train_data_std_matrix;
    fea_Train = train_scale_data';
    clear train_scale_data;

    test_data_std_matrix = ones(test_size,dim)*diag(train_data_std);
    test_data_mean_matrix = ones(test_size,dim)*diag(train_data_mean);
    clear train_data_mean train_data_std;
    test_data = (fea_Test' - test_data_mean_matrix)./test_data_std_matrix;
    clear test_data_std_matrix test_data_mean_matrix;
    fea_Test = test_data';
    clear test_data;
    %%%%%%%%%%%%%%%%%%%%%%%Normalization ends%%%%%%%%%%%%%%%%%%%%%%%

杰哥 2014-12-05 12:32 发表评论

杰哥 — Mon, 14 Jul 2014 07:47:00 GMT

http://featureselection.asu.edu/software.php中的Information Gain�Q�运行会提示�Q�错误：Undefined variable "weka" or class "weka.attributeSelection.InfoGainAttributeEval".

�q�个软�g包是*.jar文�g�Q�原以�ؓ解压��p��扑ֈ�安装文�g�Q�但没有�?/span>Weiqiang讲要安装java软�g�?/span>Feature Selection Package软�g安装�q�程如下

�W�一�?/span>�Q�安�?/span>java(使用我电脑中�?/span>“JAVA 虚拟�?/span>”)�Q��ƈ配置环境变量�Q�计��机->属�?/span>->高��pȝ��讄��->环境变量->lenovo的用户变�?/span>->Path(如果没有PATH�Q�新��Z��个变量PATH卛_��)->�~�辑->(加入”; C:\Program Files (x86)\Java\jre1.6.0_02\bin”�?#8220;变量�?#8221;,注意此处路径名可能不�?/span>)�Q?/span>

�W�二�?/span>�Q�安�?/span>FeatureSelection-Installer.jar�?/span>E: \Feature Selection Package�Q?/span> �q�行cmd->E:(定位到E�?;cmd->cd E:\Kingsoft\other\matlab 2007a\work\Feature selection\说明(目的�Q�定位到FeatureSelection-Installer.jar所在目�?/span>) ->java -jar FeatureSelection-Installer.jar

在安装过�E�中会有提示�Q?/span>

what to do once you've installed Feature Selection Package:

1.   Open MATLAB

2.   Use MATLAB's 'cd' function to navigate to the folder you

       installed Feature Selection Package into.

3.   Run the load_fspackage function with the command: 'load_fspackage'

4.   Select all of the commands you have run in this session, right click,

       and select 'Create Shortcut'. This will give you a shortcut which

       you can use to open Feature Selection Package at your choosing.

只要到这里的2�Q?/span>3�?/span>4不要��，Weiqiang也不知此�?/span>4步代表啥含义�Q�没��这一步�?/span>

�W�三�?/span>�Q�将load_fspackage.m的最开始一�?/span>”curPath = pwd;”(pwd的含义查matlab帮助卛_��Q�就是当前目�?/span>)中将pwd换成load_fspackage.m所在的目录'E:\Feature Selection Package\fspackage'卛_��?/span>在要使用该工具包�?/span>�E�序的最开始加上加�?/span>run('E:\Feature Selection Package\fspackage\load_fspackage.m')�Q�这样就不会出现如下错了�Q?/span>Undefined variable "weka" or class "weka.attributeSelection.InfoGainAttributeEval".

已经按照上述步骤�Q�在我台式机成功安装FeatureSelection-Installer.jar软�g(台式机原来没安装java软�g)�?/span>This is with Weiqiang Ren's help.

杰哥 2014-07-14 15:47 发表评论

杰哥 — Sat, 02 Nov 2013 12:00:00 GMT

http://blog.sina.com.cn/s/blog_401aa0540100ufl1.html

Matlab �l�图�Ӟ�� 可以通过讄�� "interpreter" 属�?昄��数学公式�?/div>
E.g.

xlabel('$n$','interpreter','latex');

ylabel('${\gamma}$','interpreter','latex');

title('${\gamma}(n)$','interpreter','latex');

text(0.5,0.5, '$\int_{a}^{b}f(x)dx$', 'interpreter','latex');

�?�Q?/span>我的�E�序Flicker.m (所有字体都是Latex)
h = legend('ours','$l_p$-norm MKL',1);

set(h,'Interpreter','latex');

�?�Q?/span>我的�E�序AR_Sparse.m (仅仅公式字体是latex)

h = legend('l1','DLSR-FS',4);

h1 = findobj(get(h,'Children'),'String','l1');

set(h1,'String','$l_1$','Interpreter','LaTex');

---------------------------如何在matlab中的xlabel�Q�ylabel�Q�legend和text函数中��用latex(以下关于legend的不必看�Q�已到我的例�?/span>)---------------------------
http://blog.sina.com.cn/s/blog_6e0693f70100nj22.html
以下例子中展�C�Z��如何用在matlab函数中��用latex
t = 0:0.1:2*pi;
x = sin(t);

hold on;

plot(t,x,'-*');

plot(t,2*x,'-.');

%% Using Latex in xlabel and ylabel

xlabel('$Time$','Interpreter','LaTex');

ylabel('$Value$','Interpreter','LaTex');

%% Using Latex in legend

h = legend('sin(x)__','2*sin(x)__');

h1 = findobj(get(h,'Children'),'String','sin(x)__');

set(h1,'String','$sin(\hat{x})$','Interpreter','LaTex');

h2 = findobj(get(h,'Children'),'String','2*sin(x)__');

set(h2,'String','$2*sin(\hat{x})$','Interpreter','LaTex');

%% Using Latex in text

text('Interpreter','latex',...

'String','$$\int_0^x\!\int_y dF(u,v)$$',...

'Position',[3 1],...

'FontSize',16);

%% Using Latex in title

title('$How \ to \ use \ latex \ in \ figure$','Interpreter','LaTex');

杰哥 2013-11-02 20:00 发表评论

元胞(cell)与结构体(struct)

杰哥 — Sat, 21 Sep 2013 03:31:00 GMT
元胞数组cell array

Element:cell;以下标index讉K��cell�Q�以元胞内编址content addressing讉K��元胞内容�Q�cell中可以存放�Q何类型，��M��大小的数�l?/div>

cell的创�?/span>

cell()%创徏元胞数组

c=cell(2);%创徏2×2�?/div>
c=cell(m,n);%创徏m×n�?用cell函数创徏元胞数组�Q�创建的数组为空元胞�?span style="color: red">cell函数创徏�I�元胞数�l�的主要目的是�ؓ数组预先分配�q�箋的存储空��_��节约内存占用�Q�提高执行效率。同cppblog/MATLAB�E�序的优化、预先分配存储空�?/span>

c_str=char('This is cell array');c_R=rand(3,3);c_comp=1+2i;c_sym=sym('sin(-3*t)*exp(-t)');

两种赋值方式：index赋值和content addressing赋�?/span>

index赋��|��A(1,1)={c_str};A(1,2)={c_R};A(2,1)={c_comp};A(2,2)={c_sym};

content addressing赋�?B{1,1}=c_str;B{1,2}=c_R;B{2,1}=c_comp;B{2,2}=c_sym;% The class({c_str}) and class(c_str) are cell and string, respectively.

cell array的访�?br />B{1,2}
B{1,2}(1,2)

问题�Q?/span>我有一个cell变量。比方说3*1的size  每一个里面存了一�?�?的矩��c��matlab中有没有快速的语句把每一�?*9的矩阵，reshape�?*3的矩��c��也��是说对一个cell变量中的每一个元�?应该是一个矩�?�q�行reshape操作。当然不想用循环做�?br />�{�案�Q?/span>A = {rand(1,9),rand(1,9),rand(1,9)}; cellfun(@(x) reshape(x, 3, 3).', A, 'UniformOutput', false)

cell转化成矩�늚�函数cell2mat
例：A = {rand(1,9),rand(1,9),rand(1,9)};cell2mat(A)

�l�构体数�l?structure array

Element:structure

域访�?/div>
域中可以存放��M��c�d��、�Q何大��的数组

�c�C

cell和struct的�{换cell2struct.

Matlab提供了两�U�定义结构的方式�Q�直接应用和使用struct函数�?/span>
1. 使用直接引用方式定义�l�构
与徏立数值型数组一��P��建立新struct对象不需要事先申明，可以直接引用�Q�而且可以动态扩充。比如徏立一个复数变量x�Q?br />x.real = 0; % 创徏字段名�ؓreal�Q��ƈ��字段赋��gؓ0
x.imag = 0 % 为x创徏一个新的字�D�imag�Q��ƈ��字段赋��gؓ0
x =
real: 0
imag: 0
然后可以��旗动态扩充�ؓ数组�Q?br />x(2).real = 0; % ��x扩充�?×2的结构数�l?br />x(2).imag = 0;
在�Q何需要的时候，也可以�ؓ数组动态扩充字�D�，如增加字�D�scale�Q?br />x(1).scale = 0;
�q�样�Q�所有x都增加了一个scale字段�Q�而x(1)之外的其他变量的scale字段为空�Q?br />x(1) % 查看�l�构数组的第一个元素的各个字段的内�?
ans =
real: 0
imag: 0
scale: 0
2. 使用struct函数创徏�l�构
使用struct函数也可以创建结构，该函��C�生或吧其他�Ş式的数据转换为结构数�l��?br />struct的��用格式�ؓ�Q?br />s = sturct('field1',values1,'field2',values2,…);
该函数将生成一个具有指定字�D�名和相应数据的�l�构数组�Q�其包含的数据values1、valuese2�{�必��Mؓ��h��相同�l�数的数据，数据的存放位�|�域其他�l�构位置一一对应的。对于struct的赋值用��C��元胞数组。数�l�values1、values2�{�可以是元胞数组、标量元胞单元或者单个数倹{��每个values的数据被赋值给相应的field字段�?br />当valuesx为元胞数�l�的时候，生成的结构数�l�的�l�数与元胞数�l�的�l�数相同。而在数据中不包含元胞的时候，得到的结构数�l�的�l�数�?×1的。例如：
s = struct('type',{'big','little'},'color',{'blue','red'},'x',{3,4})
s =
1x2 struct array with fields:
type
color
x

参考文献：http://blog.sciencenet.cn/blog-436588-320694.html

杰哥 2013-09-21 11:31 发表评论

matlab fscanf

杰哥 — Fri, 13 Sep 2013 06:34:00 GMT
�l�合我的教材P31
matlab文�g操作�Q�见电脑目录Blessing of Dimensionality\Features\code\readFea.m �?findImsFeaIdx.m
http://blog.sina.com.cn/s/blog_4cfb5a6201015i8q.html

�?�Q?br />路径+文�g名：d:\moon.txt
内容�Q?br />
13,1,3.4
3,2.1,23
1,12,2
4,5.4,6
现在��Z��d��moon中的数据存在一个数�l�里�Q�可以用如下�Ҏ��

fid=fopen('d:\moon.txt');

data=fscanf(fid,'%f,%f,%f',[3,inf]) ;%�q�里得用单引�?/span>

fclose(fid);

�q�时data中的数据如下�Q?/span>
13 3 1 4
1 2.1 12 5.4
3.4 23 2 6

�?�Q?br />数据在d:\test.txt
0.00 good 2
0.10 bot 3
1.02 yes 4
1.00 yes 5
1.00 yes 6
1.00 yes 3
1.00 yes 5
1.00 yes 6
1.00 yes 1
1.00 yes 3
1.00 yes 7
1.00 yes 3
1.00 yes 2

�E�序�Q?/span>
fid = fopen('d:\test.txt', 'r');
a = fscanf(fid, '%f %*s %d ', [2 inf]) % It has two rows now.
fclose(fid);
a
解释�?/span>�Q�第一列和�W�二列之间有四个�I�格�Q?span style="color: red;">format也要四空格哦�Q�有三列即三�U�类型，要有三种format�Q?span style="color: red;">%*s即�ؓ不输出字�W�串型�?br />

fid = fopen('E:\temp\test.txt', 'r');
a = fscanf(fid, '%f %*s %*f ', 5) % It has five rows and one column now. %*s %*f 卌��两个不输�?/span>
fclose(fid);
a

�E�序�l�果为：
a =
0
0.1000
1.0200
1.0000
1.0000

杰哥 2013-09-13 14:34 发表评论

matlab中对角阵的存储，diag vs spdiags

杰哥 — Sun, 07 Jul 2013 03:39:00 GMT

Feiping Nie的NIPS 2010公布代码L21R21�W?4行用�?/span>D = spdiags(d,0,n,n);可用D =diag(d)实现�c�M��功能�Q�但可能效率更低�Q�在54服务器FRGC 30Train�Q�前者耗时22�U�，后者耗时53�U�。因为前者采用稀疏存储，更节省内存。Nie老师�Ԍ��功能相同�Q�但�׃��是sparse的矩阵，矩阵很大时会快很多�?/span>
��单例子：
p=[1 2 3]; D1=diag(p);D2 = spdiags(p',0,3,3);

杰哥 2013-07-07 11:39 发表评论

杰哥 — Thu, 04 Jul 2013 06:28:00 GMT

[dim,num]=size(Data);weight是一�?×num的行向量�Q�够成一个对角阵P=diag(weight);目标�Q�实现Data*P
方式一�Q�Data*diag(weight)�Q�这�U�需要dim×num×(2num-1)个加法或乘法�Q?br />方式�?/span>�Q�Data.*repmat(weight,dim,1)�Q�这�U�只需要dim×num个乘法�?br />�q�就是Ran He的CVPR 12代码CRFS�W?1行用:X1 = Data.*repmat(weight,dim,1)实现Data*P的原因�?/div>

杰哥 2013-07-04 14:28 发表评论

能不用inv��量不用

杰哥 — Wed, 03 Jul 2013 14:47:00 GMT
Ran He老师讲他CVPR 2012公式20不是用inv求的�Q�他和学生说�q�能不用inv�Q�尽量不用，用左除或者右除，解线性方�E�组。不能受我matlab教材影响�Q�左除或者右除和求inv是不�{��h的。Feiping Nie的NIPS 2010公布代码L21R21�W?5�?span style="color: red;">实现一个矩�늚�逆乘以另外一个矩阵，也用右除。类似实��C��个矩阵乘以另外一个矩�늚�逆，用左除，不用inv

http://zhidao.baidu.com/question/308277515.html
左除要不要求矩阵可逆？��Z��么一个不可逆的矩阵a在matlaB中用inv�Q�a�Q?y不可计算�Q�用a\y��可以计��了�Q?/div>
�{�：左除不要求矩阵可逆，在matlab中inv�Q�a�Q?y与a\y�q�不�{�同�Q�左除命令有更深的内涵，当a为方阵且可逆时�Q�二者都能得到正��结果，�?span style="color: red;">左除�_�ֺ�更高�Q�运��数度快(我实验下来，��实是这�?�Q�特别当A为病态矩阉|��Q�inv(A)*y的计��结果不可靠�Q�不�q�这�U�矩阵出现的概率很低�Q�，但用左除计算出的�l�果却很可靠�Q�因为左除是matlab设计的一个智能算�W�，他可以智能地选择��法�Q�当方程为矛盾方�E�（也就是独立方�E�个数多余未知数个数�Q�时�Q�此时显然没有精��解�Q�但很多时候我们需要x在最��二乘意义下的解�Q�此时用A\y��p��得到x的最��二乘解�?br />
Haiqing Li说他从来不看�q�些中文的，有些人可能回�{�得不专业或误解�Q�到matlab帮助搜烦arithmetic operators有左除或者右除详�l�解�?br />
��试例子�Q?br />�W�一个matlab的m文�g
% 生成数据
rand('state',0);
n = 10000;% 也可讄��?000
A = randn(n,n);
B = randn(n,1);
save MyFile1 A B;

�W�二个matlab的m文�g
% ��试求逆时�?/div>
tic;
load MyFile1;
n = size(A,1);
C1 = inv(A'*A+0.01*eye(n))*B;
save MyFile2 C1;
toc;

�W�三个matlab的m文�g
% ��试左除旉��
tic;
load MyFile1;
n = size(A,1);
C2 = (A'*A+0.01*eye(n))\B;
save MyFile3 C2;
toc;

杰哥 2013-07-03 22:47 发表评论

如何用matlab求矩�늚�最大奇异��|��(Matlab code to solve the largest singular value of a matrix)

杰哥 — Fri, 03 May 2013 13:44:00 GMT

n = norm(A) returns the largest singular value of A, max(svd(A)).

杰哥 2013-05-03 21:44 发表评论

杰哥 — Thu, 13 Dec 2012 00:30:00 GMT

�~�写NFL和NFP分类器应用于ASLAN数据库发现combntns函数所引�v的内存不��I��到网上搜索有人提问，但没有解��x��案�?br />问题�Q�我惛_��一个序列的��L��l�合�Q�用combntns�Q�在数目较小�Ӟ��可以辑ֈ�目的�Q�但是数目变大时�Q�就不行了，求助各位�?/p>
例如�?:5,里面3个数字的�l�合�Q�很�Ҏ��得到�?br />combntns(1:5,3) ��可以了
但是�?:100里面10个数字的�l�合�Q�就不行了�?br />我的解决�Ҏ��Q?/span>如果你未必要取所有组合，可以�?�Q?00中�Q�?0个数�Q�再取这20个数字里�?0个数字的�l�合�?br />%matlab code:
rand('state',0);
temp1 = randperm(100);
temp2 = temp1(1:20);
combntns(temp2,10)

杰哥 2012-12-13 08:30 发表评论

杰哥 — Wed, 05 Dec 2012 11:32:00 GMT

http://www.ilovematlab.cn/thread-68860-1-1.html
在win7�Q?2位家庭版�?有两�U�方法解��x��常运行问题，一是将配色�Ҏ��改�ؓ�l�典(在桌面空白处叛_��鼠标�Q�选择个性化�Q�在弹出的对话框中选择windows�l�典主题卛_��)�Q�二是兼�Ҏ��改为vista sp1或sp2.卸蝲完了后，再把主题改回��C��原来的主题�?br />
��Z��么matlab2010a只能在断�|�时才能打开�Q?/span>
20121207发现matlab打不开�Q�上�|�搜索解��x��案。有��断网��可以。我实验两次�Q�注销后断�|�能打开matlab�Q�注销后连�|�，不打开其他��M��E�序�Q�matlab�q�是打不开�Q�等了五分钟。能看到matlab�q�程。��l�搜�?#8220;Windows 7 matlab断网才能打开”。查到解��x��?http://wenwen.soso.com/z/q294348187.htm)�Q?span style="font-family: Arial, Helvetica, sans-serif, SimSun; line-height: 22px;">你可能用的是window7吧，�q�种现象的��生是�׃��MatLab与Windows 7 的兼�Ҏ��问题，解决的方法如下：从纯净的Windows XP�pȝ��中system32目录中拷贝一份iphlpapi.dll�?Matlab 2010a安装目录\bin\win32中问题即可解冻I��当然�Q�X64的系�l�要�?Windows XP X64 中去拗��我问Lianyang Ma要了�q�个文�g�Q�重启下电脑��搞好了(即��在同时打开qq,gtalk,��h��览�?FileZilla的情况下�Q�也能打开matlab)�? �?/span>配色�Ҏ��和兼�Ҏ��没有关�p�R�?/div>

杰哥 2012-12-05 19:32 发表评论

How to use matlab solve optimization quadratic?

杰哥 — Wed, 21 Nov 2012 10:31:00 GMT
Nannan gives me a fold named "Matlab Help". On page 46 of "Optimization Toolbox User Guide", it lists the constrain and objective type, and the matlab function. For example, if the constrain is linear and the objective is quadratic, we can use quadprog. Note that it can not slove $D_1$ in Section 4.1 of "Smooth minimization of non-smooth functions". Problem: max ((X^T)HX) and H is positive semi definite. The matlab function "quadratic" can not solve this kind of problem. It can only solve the problem: min ((X^T)HX) and H is positive semi definite.

杰哥 2012-11-21 18:31 发表评论

How to compile the mexfile?

杰哥 — Thu, 18 Oct 2012 00:27:00 GMT

An Efficient Projection for l1,∞ Regularization该文提供代码�Q�在Readme文�g中说To run the projection in Matlab, you first need to compile the mexfile:

$ mex -output projL1Inf mex-projL1Inf.c

Then, from Matlab, test the projection by running the script "runprojection",

(1)在我�W�记本电脑，打开matlab�Q�打开该程序所在的文�g夹，在matlab命��o行输入mex -setup�Q?span style="color: red">matlab输出�Q?/span>Would you like mex to locate installed compilers [y]/n? 手动输入�Q?/span>y �Q?span style="color: red">matlab输出�Q?nbsp;
Select a compiler:  手动输入�Q?, matlab输出�Q?/span> Are these correct [y]/n?   手动输入�Q?/span>y �?span style="color: #ff0000">手动输入�Q?/span> mex -output projL1Inf mex-projL1Inf.c �Q�可能是因�ؓ我的�W�记本电脑没装VC�Q�有错误。在我本��能运行，关键��是生成projL1Inf.mexw32�q�个文�g��可以了�Q�在Mingming Gong电脑上生成，拯��到我的电脑就可以了�?/span> Mingming电脑装VC了，但还是不能编译，他讲可能原因是这个c语言代码本来��是在linux下写的，在windows上跑不了�Q�naiyang也这么说�?span style="color: red">20121023, after Mingming Gong saw my blog, he found the solution to produce the file projL1Inf.mexw32. 只要��mex-projL1Inf.c后缀改�ؓmex-projL1Inf.cpp�Q�projL1Inf.c后缀改�ؓprojL1Inf.cpp�Q�mex-projL1Inf.cpp�W?9行改�?include "projL1Inf.cpp"�Q�只要修改这三点��可以在windows下编译，我已�l�从他那里拷贝了projL1Inf.mexw32 。从我本��是不能编译，因�ؓ本机没装VC�?span style="color: #ff0000">�?/span>因：*.cpp是支持C和C++的；*.c只支持C�Q�故//不支持。在Linux下C�?/span>C++是两个编译器(在linux下c是用gcc�~�译�Q�C++使用g++�~�译)�Q�在windows下，C和C++是可以同时用VC�~�译的�?br />
20130629在NLPR�?4服务器用mex -output projL1Inf mex-projL1Inf.cpp成功�~�译生成文�gprojL1Inf.mexw64�Q�不知�ؓ什么用mex -setup一步一步来不行�?br />
(2)在uts cluster�Q�按照上面输入的命��o�Q�有错误和警�?(Mingming said that the warning does not matter)

Warning: You are using gcc version "4.4.6". The version

currently supported with MEX is "4.3.4".
错误�Q�projL1Inf.c:248: error: expected expression before '/' token
projL1Inf.c:250: error: expected expression before '/' token

projL1Inf.c:252: error: expected expression before '/' token
打开projL1Inf.c �Q�将�q�三行的注释//�Q�均换成/* * /
�q�样的作用，��是生成projL1Inf.mexa64�q�个文�g.如果用命令mex mex-projL1Inf.c,��是生成mex-projL1Inf.mexa64 (已经实践�Q?/span>��原�?/span> projL1Inf.mexa64 删除�Q�用命��omex mex-projL1Inf.c,��是生成mex-projL1Inf.mexa64�Q?/span> 再将�q�个文�g改名为projL1Inf.mexa64�Q�再�q�行runprojection �q�是可以�?

20121208在uts服务器编译libsvm-mat-2.9-1�Q�编译有问题�Q�根据提�C�，��所有文�?.c都变�?.cpp,�q�有问题g++: svm.obj: No such file or directory�Q�上�|�搜索解��x��案http://blog.csdn.net/kit_147/article/details/7417734�Q�将文�gsvm.obj改�ؓsvm.o�Q�编译就没有问题了，但用�E�序��试�q�有问题�Q�提�C?"-largeArrayDims"的问题，make.m开始注释：% add -largeArrayDims on 64-bit machines�Q�问mingming gong怎么add�Q�将makefile中的mex...修改为mex...-largeArrayDims卛_��。例�?br />
mex -O -c svm.cpp

mex -O -c svm_model_matlab.c
修改�?/span>

mex -O -c svm.cpp -largeArrayDims

mex -O -c svm_model_matlab.c -largeArrayDims
(This is with Mingming Gong's help)

20150710在跑FastHash提示错误�Q�binaryTreeTrain1.mexw64不是有效的w32应用�E�序,感到费解�Q�这�?4位机�Q�怎么会提�C��个错误。Jian Liang讲重新编译下��p��了，他就在matlab�H�口敲入mex binaryTreeTrain1.cpp�Q�再敲入mex forestInds.cpp��可以了

杰哥 2012-10-18 08:27 发表评论

RGB Triplet	Short Name	Long Name
`[1 1 0]`	`y`	`yellow`
`[1 0 1]`	`m`	`magenta`
`[0 1 1]`	`c`	`cyan`
`[1 0 0]`	`r`	`red`
`[0 1 0]`	`g`	`green`
`[0 0 1]`	`b`	`blue`
`[1 1 1]`	`w`	`white`
`[0 0 0]`	`k`	`black`

亚洲人成77777在线观看网,亚洲国产日韩欧美在线99,国产欧美在线播放

How to read *.data in Matlab and Python

matlab legend字号改变

[转蝲]matlab中text 函数在显�C�字�W�串时的使用�Ҏ��

matlab ColorSpec (Color Specification)

Professor Deng Cai's code, very good code

[zz]MATLAB如何实现十进制数与二�q�制数的转换

matlab dir

The matlab code of Locality sensitive hashing (LSH)

matlab函数meshgrid

matlab函数fminunc:求无�U�束优化的最优解

matlab加�? bsxfun

matlab函数cat

How to solve AX + XB = C for X using matlab?

f-measure matlab

The matlab code of canonical correlation analysis (CCA)

[zz]matlab囑փ�滤�L

The matlab code to normalize each feature of training data and test data to unit mean and zero variance

元胞(cell)与结构体(struct)

matlab fscanf

matlab中对角阵的存储，diag vs spdiags

能不用inv���量不用

如何用matlab求矩�늚�最大奇异��|��(Matlab code to solve the largest singular value of a matrix)

How to use matlab solve optimization quadratic?

How to compile the mexfile?

能不用inv��量不用