%对数据进行标准化
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,
!["Smile"]("images/smiles/smile.gif")
.*(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 请搂主给我发一份好吗？万分谢谢
ren.jijun@gmail.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 发给我吧，谢谢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
