Ica盲源分离Matlab程序closeall;clearall;i4=imread('1.jpg');i5=imread('2.jpg');i6=imread('2.png');i1=rgb2gray(i4);i2=rgb2gray(i5);i3=rgb2gray(i6);s1=reshape(i1,[1,256*256]);s2=reshape(i2,[1,256*256]);s3=reshape(i3,[1,256*256]);s=[s1;s2;s3];sig=double(s);aorig=rand(size(sig,1));mixedsig=aorig*sig;ms1=reshape(mixedsig(1,:),[256,256]);ms2=reshape(mixedsig(2,:),[256,256]);ms3=reshape(mixedsig(3,:),[256,256]);figure;subplot(331),imshow(i1),subplot(332),imshow(i2),subplot(333),imshow(i3);subplot(334),imshow(uint8(ms1)),subplot(335),imshow(uint8(ms2)),subplot(336),imshow(uint8(ms3));%mixedsig=zeros(size(mixedsig));meanValue=mean(mixedsig')';mixedsig=mixedsig-meanValue*ones(1,size(mixedsig,2));covarianceMatrix=cov(mixedsig',1);[E,D]=eig(covarianceMatrix);eigenvalues=flipud(sort(diag(D)));whiteningMatrix=inv(sqrt(D))*E';dewhiteningMatrix=E*sqrt(D);whitesig=whiteningMatrix*mixedsig;X=whitesig;[vectorSize,numSamples]=size(X);B=zeros(vectorSize);numOFIC=vectorSize;forr=1:numOFICi=1;maxNumIterations=100;w=rand(vectorSize,1)-.5;w=w/norm(w);whilei<=maxNumIterations1w=w-B*B'*w;w=w/norm(w);w=(X*((X'*w).^3))/numSamples-3*w;w=w/norm(w);i=i1;endA(:,r)=dewhiteningMatrix*w;W(r,:)=w'*whiteningMatrix;B(:,r)=w;endicasig=W*mixedsig(W*meanValue)*ones(1,numSamples);icasig=abs(50*icasig);a=reshape(icasig(1,:),[256,256]);b=reshape(icasig(2,:),[256,256]);c=reshape(icasig(3,:),[256,256]);subplot(337),imshow(uint8(a)),subplot(338),imshow(uint8(b)),subplot(339),imshow(uint8(c))
本文档为【Ica盲源分离Matlab程序(2)】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。
该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。
[版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
浙公网安备 33021202002483