function?outputimage=RX（A）
[mnq]=size（A）;
newimage=?reshape（?Am*nq）;
lastimage=newimage‘;
meanRow=mean（lastimage2）;%求每一行的均值
%?meanRow=round（meanRow）;%取整
%?C=zeros（qq）;
%?for?i=1:m*n?
%??????lastimage（:i）=?lastimage（:i）-meanRow;
%?????
%?C=?lastimage*lastimage‘+C;?
%?Csum=C;
%?end
%?Wcov=Csum/（m*n）;%協方差矩陣

for?i=1:m*n?
?????lastimage（:i）=?lastimage（:i）-meanRow;
end
Wcov=cov（lastimage‘）;
invWcov=inv（Wcov）/（m*n-1）;


rx=zeros（1m*n）;
for?i=1:m*n
????rx（i）=lastimage（:i）‘*?invWcov?*lastimage（:i）;??
end
outputimage=reshape（rxmn）;
outputimage=mat2gray（outputimage）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????627??2014-01-02?09:26??RX.m

-----------?---------??----------?-----??----

??????????????????627????????????????????1