模式识别中运用ICA进行特征提取的matlab代码

FastICA算法代码
%对数据进行标准化
close all,clear
load np_all
normsr=zscore(npp);
features=normsr";
[r,c]= size(features);
%对数据进行白化，使之为0均值，单位协方差
features = features - mean(features")"*ones(1,c);
[E, D]  = eig(cov(features",1));
V      = E*inv(sqrtm(D))*E";
z = V*features;
%使用FastICA算法对数据进行特征提取.非线性函数选择tanh(a1y)
epsilon=0.0001;
m=r;
p=1;
W=zeros(r,1);
for p=1:m
W(:,p)=rand(r,1);
W(:,p)=W(:,p)/norm(W(:,p));
exit=0;
count=0;
iter=1;
while exit==0;
count=count+1;
temp=W(:,p);
for i=1:m
W(i,p)=mean(z(i,.*(tanh((temp)"*z)))-(mean(1-(tanh((temp))"*z).^2)).*temp(i,1);
end
ssum=zeros(r,1);
for counter=1:p-1
ssum=ssum+(W(:,p)"*W(:,counter))*W(:,counter);
end
W(:,p)=W(:,p)-ssum;
W(:,p)=W(:,p)/norm(W(:,p)); 数据挖掘研究院
if(abs((dot(W(:,p),temp)))<1+epsilon)&(abs((dot(W(:,p),temp)))>1-epsilon)
exit=1;
end
iter=iter+1;
end
end
W=W";
%找到最独立的方向
power=sum(abs(W)");
[m, in]  = sort(power);
W= W(in(r-2+1:r),;
%计算新的特征数据.
features = W*z;
features=features";
X=features(find(kk==1|kk==2),;
group=kk(find(kk==1|kk==2));
%画散点图.
gscatter(X(:,1),X(:,2),group,"mr","o+")
legend("health","npc-initial");
xlabel("1st component");
ylabel("2nd component");

怎么还不行呢?

非常好！请问能不能发给我啊！nj_szm@yahoo.com.cn

楼主，积分太低，现在急需啊，帮忙发给我一份吧，感激不尽！！
邮箱：jinselanling@163.com

技术分不够，可不可以发到我邮箱里？andybeyond@sina.com

请搂主给我发一份好吗？万分谢谢

数据挖掘交友

楼上的几位,程序已经群发,注意接受.

老大可以发一份给我吗？我也想看看，但技术分太低～
dinga@yeah.net

楼主发一份给 我吧。zqm20021009@sina.com
谢谢了！

好，顶！！！！！！！！！！！！！

楼主可以发一份给我吗？azangly@sina.com
谢谢！

我也想要　发份给我吧　zheshimajia@yahoo.com.cn 谢谢

哪个群发啊 ？
怎么找得到啊

楼主，发给我可以吗 ？急用啊！！！谢谢了
ilywenwenzi@163.com

楼主，积分太低，看不到，急需，帮忙发给我一份吧，感激不尽！！
邮箱：ywb-dj@163.com

楼主发一份给 我吧。liuliyz@126.com 谢谢了！ 万分感谢!!!!

楼主发一份给 我吧。liuliyz@126.com 谢谢了！ 万分感谢!!!!

技术分不够，可不可以发到我邮箱里？

技术分不够，可不可以发到我邮箱里~happy12338@yahoo.com.tw

能给我一份吗?非常的需要：）
candyfeifei@sjtu.edu.cn

请给我一份,技术分不够啊

aiwu818@163.com

我想要，积分太低，

发给我吧，谢谢chenfengjun87602@yahoo.com

楼主，积分太低，看不到,帮忙发给我一份吧，感激不尽！！
邮箱：gfengyang@163.com

技术分底，可是等着急用，能不能发到evenyw@21cn.com
不胜感激！

楼主，积分太低，看不到,帮忙发给我一份吧，感激不尽！！
邮箱：hlg5141979@yahoo.com.cn

我也想要，可是看不到
chenfengjun87602@yahoo.com

楼主,可否给我发一份,winny127@163.com,谢谢!

我的积分也太底　　　　给发一个吧！！
gjkkk@163.com

楼主，急需你的程序，但是技术分太低，能不能麻烦您发一封到我的邮箱，gigi1981◎sina.com
万分感谢！

我现在急需这个完成论文的仿真部分，真诚的希望得到您的帮助，ismexiaoli@yahoo.com.cn 劳请您发给我一份，谢谢！

楼主，能传一份吗？
EMAIL：moon_star_1982@163.com

技术分太低，看不到呀，
作毕业设计急用，楼主能不能给我发一份呀
355388@163.com
谢谢

希望您能给我发一分,你就是我的恩人!

希望您能给我发一分,你就是我的恩人! 数据挖掘交友
erduizhang@163.com

给我发一份吧,谢谢 tangjingyuan0001@163.com

楼主，请给我发一份吧,我现在正为特征筛选判决头痛呢，谢谢 ！yinjiajian_100@163.com

楼主，积分太低，现在急需啊，帮忙发给我一份吧，感激不尽！！
邮箱：wangrui3724@163.com

请问能不能发给我啊！ ren.jijun@gmail.com

谢谢,我先运行看看

楼主 能否发给我一份
mengqingyu18@sina.com
非常感谢

哦，谢谢共享！

LZ.可以也麻烦给我发一份吗？我现在正在做特征提取，但是效果很不理想，可以参考下你的吗？谢谢了。我的技术分太低，看不了。麻烦你了
ltma1821@sina.com

好心的楼主，也给我发一份吧
yanhai_cen@yahoo.com.cn

楼主,我现在急需,帮忙发一份给我吧,非常谢谢!!
邮箱:tangjingyuan0001@163.com

楼主，积分太低，现在急需啊，帮忙发给我一份吧，感激不尽！！
邮箱yxmyxm2001◎163.com

Create By Any-Extract(WL-AE)

数据挖掘论坛