91av视频/亚洲h视频/操亚洲美女/外国一级黄色毛片 - 国产三级三级三级三级

  • 大小: 2KB
    文件類型: .rar
    金幣: 2
    下載: 1 次
    發(fā)布日期: 2021-12-01
  • 語(yǔ)言: 其他
  • 標(biāo)簽: 432??

資源簡(jiǎn)介

完美分離各種混合音頻,并且有輸出功能,并且算法在分離性能上比傳統(tǒng)算法高效,

資源截圖

代碼片段和文件信息

function?ICAedS?=?fastica_3_1(MixedS)
%%%%%%%%%%%%%%%%%%%%%%%%%%??標(biāo)準(zhǔn)化??%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
MixedS_mean=zeros(31);
for?i=1:3
????MixedS_mean(i)=mean(MixedS(i:));
end????????????????????????????????????????%?計(jì)算MixedS的均值

for?i=1:3
????for?j=1:size(MixedS2)
????????MixedS(ij)=MixedS(ij)-MixedS_mean(i);
????end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%??白化??%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

MixedS_cov=cov(MixedS‘);????????????????????%?cov為求協(xié)方差的函數(shù)
[ED]=eig(MixedS_cov);??????????????????????%?對(duì)信號(hào)矩陣的協(xié)方差函數(shù)進(jìn)行特征值分解
Q=inv(sqrt(D))*(E)‘;????????????????????????%?Q為白化矩陣
MixedS_white=Q*MixedS;??????????????????????%?MixedS_white為白化后的信號(hào)矩陣
IsI=cov(MixedS_white‘);?????????????????????%?IsI應(yīng)為單位陣????????????

%%%%%%%%%%%%%%%%%%%%%%%% FASTICA算法??%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
X=MixedS_white;????????????????????????????%?以下算法將對(duì)X進(jìn)行操作
[VariableNumSampleNum]=size(X);
numofIC=VariableNum;???????????????????????%?在此應(yīng)用中,獨(dú)立元個(gè)數(shù)等于變量個(gè)數(shù)
B=zeros(numofICVariableNum);??????????????%?初始化列向量w的寄存矩陣B=[b1??b2??...???bd]
for?r=1:numofIC
????i=1;maxIterationsNum=100;???????????????%?設(shè)置最大迭代次數(shù)(即對(duì)于每個(gè)獨(dú)立分量而言迭代均不超過(guò)此次數(shù))
????k=0;
????b=rand(numofIC1);??????????????????%?隨機(jī)設(shè)置b初值
????b=b/norm(b);
????%HHH=b;
????%?對(duì)b標(biāo)準(zhǔn)化?norm(b):向量元素平方和開(kāi)根號(hào)
????while?i<=maxIterationsNum+1
????????if?i?==?maxIterationsNum???????????%?循環(huán)結(jié)束處理
????????????fprintf(‘\n第%d分量在%d次迭代內(nèi)并不收斂?!?rmaxIterationsNum);
????????????break;
????????end
????????bOld=b;??????????????????????????
????????a2=1;
????????t=X‘*b;
????????g=t.*exp(-a2*t.^2/2);
????????dg=(1-a2*t.^2).*exp(-a2*t.^2/2);
%?????????b=((1-u)*t‘*g*b+u*X*g)/SampleNum-mean(dg)*b;
??????????????????????????%?核心公式,參見(jiàn)理論部分公式2.52
????????ava=(X*g)/SampleNum;
????????ave=mean(dg)*b;
????????b=ava-ave;
????????b=b-B*B‘*b;?
????????b=b/norm(b);
????????b0=b;
????????t2=X‘*b;
????????dg2=(1-a2*t2.^2).*exp(-a2*t2.^2/2);
????????ave2=mean(dg2)*bOld;
????????b=2*ava-ave-ave2;
????????b=b-B*B‘*b;????????????????????????%?對(duì)b正交化
????????b=b/norm(b);?
????????k=k+1;
%?????????abs(abs(b‘*bOld)-1)<1e-9
????????if?abs(abs(b‘*b0)-1)<1e-9
???????????%?jjjj=b;?????????????????????%?如果收斂,則
?????????????B(:r)=b;?????????????????????%?保存所得向量b
?????????????fprintf(‘\n第%d分量在%d次迭代?!?rk)
?????????????break;
?????????????
?????????end
????????i=i+1;????????
????end
%????B(:r)=b;????????????????????????????????%?保存所得向量b
end
%%%%%%%%%%%%%%%%%%%%%%%%%%??ICA計(jì)算的數(shù)據(jù)復(fù)原并構(gòu)圖??%%%%%%%%%%%%%%%%%%%%%%%%%
ICAedS=B‘*MixedS_white;???
end%?計(jì)算ICA后的矩陣


?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????2753??2018-09-14?09:01??f_3_1\fastica_3_1.m

?????文件???????1782??2018-09-14?09:01??f_3_1\f_3_1.m

?????目錄??????????0??2018-12-26?10:22??f_3_1

-----------?---------??----------?-----??----

?????????????????4535????????????????????3


評(píng)論

共有 條評(píng)論