clear?all
clc
%產生非線性卷積盲信號分離

%1.產生源信號
T=0.1;????????%取時間區間為[0s??0.1s]
Num=1000;?????%10k采樣率下，采了1000個點
M=3;??????????%信源數=傳感器數=3
Lum=50;???????%卷積的移位值
Kum=100;??????%STFT的窗長度
a=-0.1;b=0.1;
t=linspace（0TNum）;
S1=sign（cos（2*pi*155*t））;???%符號信號
S2=sin（2*pi*800*t）;?????????%高頻正弦信號
S3=sin（2*pi*9*t）.*sin（2*pi*300*t）;??????%幅度調制信號
S=[S1;S2;S3];

%2.產生混合矩陣A，且A為元胞數組
A=cell（MM）;??%分別取不同時刻的混合矩陣完成卷積
for?i=1:M
????for?j=1:M
????????A{ij}=a+（b-a）*rand（1Lum）;
????end
end

X=zeros（size（S））;
ans1=0;ans2=0;ans3=0;
for?i=1:M
????for?k=1:Num
????????for?j=1:M
????????????for?p=0:Lum-1
????????????????if?（k-p>0）&（k-p<=Num）
????????????????????ans1=A{ij}（p+1）*S（jk-p）;
????????????????else?ans1=0;
????????????????end
????????????ans2=ans2+ans1;
????????????end
?????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1667??2008-10-08?10:54??feixianxinghunhe.m

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

?????????????????1689????????????????????2