久久久久久久999_99精品久久精品一区二区爱城_成人欧美一区二区三区在线播放_国产精品日本一区二区不卡视频_国产午夜视频_欧美精品在线观看免费

標(biāo)題: Matlab盲源分離的一個(gè)小程序 [打印本頁(yè)]
作者: 壁虎歸來(lái)    時(shí)間: 2018-3-13 16:06
標(biāo)題: Matlab盲源分離的一個(gè)小程序
matlab盲源分離的一個(gè)小程序。。。

源程序如下:
  1. %以下程序調(diào)用ICA，輸入觀察信號(hào)，輸出為解混合信號(hào)
  2. function Z=ICA(X)
  3. %-------------去均值------------
  4. [M,T] = size(X); %獲取輸入矩陣的行/列數(shù)，行數(shù)為觀測(cè)數(shù)據(jù)的個(gè)數(shù)，列數(shù)為采樣點(diǎn)點(diǎn)數(shù)
  5. average= mean(X')';  %按行取均值
  6. for i=1:M
  7.     X(i,:)=X(i,:)-average(i)*ones(1,T);
  8. end
  9. %------------白化--------------
  10. Cx = cov(X',1);    %計(jì)算協(xié)方差矩陣Cx
  11. [eigvector,eigvalue] = eig(Cx); %計(jì)算Cx的特征值和特征向量
  12. W=eigvalue^(-1/2)*eigvector';   %白化矩陣
  13. Z=W*X;   %正交矩陣

  15. %-------------迭代-------------
  16. Maxcount=10000;        %最大迭代次數(shù)
  17. Critical=0.00001;   %判斷是否收斂
  18. m=M;                %需要估計(jì)的分量的個(gè)數(shù)
  19. W=rand(m);
  20. for n=1:m  
  21.     WP=W(:,n);  %初始權(quán)矢量（任意）
  22. %     Y=WP'*Z;
  23. %     G=Y.^3;%G為非線性函數(shù)，可取y^3等
  24. %     GG=3*Y.^2;  %G的導(dǎo)數(shù)
  25.     count=0;
  26.     LastWP=zeros(m,1);
  27.     W(:,n)=W(:,n)/norm(W(:,n));%單位化一列向量
  28.     while abs(WP-LastWP)&abs(WP+LastWP)>Critical %兩個(gè)絕對(duì)值同時(shí)大于收斂條件
  29.         count=count+1;   %迭代次數(shù)
  30.         LastWP=WP;       %上次迭代的值
  31.        % WP=1/T*Z*((LastWP'*Z).^3)'-3*LastWP;
  32.         for i=1:m     
  33. WP(i)=mean(Z(i,:).*(tanh((LastWP)'*Z)))-(mean(1-(tanh((LastWP))'*Z).^2)).*LastWP(i);%更新
  34.         end
  35.         WPP=zeros(m,1);%施密特正交化
  36.         for j=1:n-1
  37.             WPP=WPP+(WP'*W(:,j))*W(:,j);
  38.         end
  39.         WP=WP-WPP;
  40.         WP=WP/(norm(WP));

  42.         if count==Maxcount
  43.             fprintf('未找到相應(yīng)的信號(hào)');
  44.             return;
  45.         end
  46.     end
  47.     W(:,n)=WP;
  48. end
  49. Z=W'*Z;
復(fù)制代碼

所有資料51hei提供下載:
mangyuanfenli.rar (972 Bytes, 下載次數(shù): 16)






歡迎光臨 (http://www.zg4o1577.cn/bbs/) Powered by Discuz! X3.1
主站蜘蛛池模板： 伊人一区| 国产一区不卡在线观看 | 日韩精品一区二区三区在线观看 | 久久男人 | 欧美日韩国产一区二区三区不卡 | 国产成人精品一区二区三区 | 天天色官网 | 91精品一区二区三区久久久久 | 国产日韩一区二区三免费高清 | 成人精品鲁一区一区二区 | 日韩精品av一区二区三区 | 国产成人高清 | 亚洲高清一区二区三区 | 丝袜 亚洲 另类 欧美 综合 | 国产精品九九视频 | 欧美一级久久精品 | 成人精品国产一区二区4080 | 国产精品免费在线 | 中国黄色毛片视频 | 精品国产视频 | 久热免费 | 精品成人佐山爱一区二区 | 综合久久亚洲 | 99精品国产一区二区青青牛奶 | 精品国产91亚洲一区二区三区www | 亚洲国产成人精品女人久久久 | 亚洲成人免费观看 | 久久精品99 | 中文字幕一区二区三区四区五区 | 久久综合伊人 | 国产伦精品一区二区三区四区视频 | 欧美精品v国产精品v日韩精品 | 日韩中文一区 | 国产欧美一区二区三区在线播放 | 国产精品1区2区 | 黄色成人亚洲 | 羞羞视频免费观看入口 | 久干网| 亚洲中午字幕 | 日韩成人精品一区二区三区 | 99国产精品久久久久老师 |