Ica盲源分离Matlab程序.docx

下载文档

38
0
约1.48千字
约 3页
2017-04-11 发布于重庆
举报
版权申诉
保障服务

Ica盲源分离Matlab程序.docx

1、本文档共3页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

Ica盲源分离Matlab程序

Ica盲源分离Matlab程序close all;clear all;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;for r=1:numOFIC i=1;maxNumIterations=100; w=rand(vectorSize,1)-.5; w=w/norm(w); while i=maxNumIterations+1 w=w-B*B*w; w=w/norm(w); w=(X*((X*w).^3))/numSamples-3*w; w=w/norm(w); i=i+1; end A(:,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))