%連續波信號用基于影響因子的頻域盲源分離排序算法，效果較好
clear;
clc;
%源信號，50000點采樣，連續信號同時到達
Fs=100e6;%采樣率
N=50000;%總共仿真500us，采樣點數50000
t=[0:1/（100e6）:500e-6-0.01e-6];

%s1（t）產生LFM信號，中心頻率0.5M
s1=chirp（t0500e-6-0.01e-61e6‘linear‘）;%產生脈沖周期序列的基礎

%s2（t）BPSK信號，10位，1110010101，中心頻率為0.1M，每個碼元持續50us
%?b1=[1?1?1?0?0?1?0?1?0?1];
%?b=（b1-1）*（-pi）;
%?tt2=[0:1/（100e6）:50e-6-0.01e-6;
%?????50e-6:1/（100e6）:100e-6-0.01e-6;
%?????100e-6:1/（100e6）:150e-6-0.01e-6;
%?????150e-6:1/（100e6）:200e-6-0.01e-6;
%?????200e-6:1/（100e6）:250e-6-0.01e-6;
%?????250e-6:1/（100e6）:300e-6-0.01e-6;
%?????300e-6:1/（100e6）:350e-6-0.01e-6;
%?????350e-6:1/（100e6）:400e-6-0.01e-6;
%?????400e-6:1/（100e6）:450e-6-0.01e-6;
%?????450e-6:1/（100e6）:500e-6-0.01e-6];
%?????%10*5000點
%?ss2=zeros（105000）;
%?for?i=1:10
%?????ss2（i:）=cos（2*pi*0.1e6*tt2（i:）+b（i））;
%?end
%?s2=zeros（150000）;
%?for?i=1:10
%?????s2（5000*（i-1）+1:5000*i）=ss2（i:）;
%?end
%?
%s3產生2FSK信號，10位1110010101，1頻率0.5M，0頻率0.1M
c1=[1?0?1?0?1?0?0?1?0?0];
tt3=[0:1/（100e6）:50e-6-0.01e-6;
????50e-6:1/（100e6）:100e-6-0.01e-6;
????100e-6:1/（100e6）:150e-6-0.01e-6;
????150e-6:1/（100e6）:200e-6-0.01e-6;
????200e-6:1/（100e6）:250e-6-0.01e-6;
????250e-6:1/（100e6）:300e-6-0.01e-6;
????300e-6:1/（100e6）:350e-6-0.01e-6;
????350e-6:1/（100e6）:400e-6-0.01e-6;
????400e-6:1/（100e6）:450e-6-0.01e-6;
????450e-6:1/（100e6）:500e-6-0.01e-6];
????%10*5000點
ss3=zeros（105000）;
for?i=1:10
????if?c1（i）==1
????????ss3（i:）=cos（2*pi*0.5e6*tt3（i:））;
????else
????????ss3（i:）=cos（2*pi*0.1e6*tt3（i:））;
????end
end
s3=zeros（150000）;
for?i=1:10
????s3（5000*（i-1）+1:5000*i）=ss3（i:）;
end

originalsig2=[s1;s3];
figure（1）;
subplot（‘211‘）;plot（ts1）;axis（[0500e-6-22]）;title（‘s1（t）LFM信號‘）;xlabel（‘t（單位s）‘）;
%?subplot（‘312‘）;plot（ts2）;axis（[0500e-6-22]）;title（‘s2（t）BPSK信號‘）;xlabel（‘t（單位s）‘）;
subplot（‘212‘）;plot（ts3）;axis（[0500e-6-22]）;title（‘s3（t）2FSK信號‘）;xlabel（‘t（單位s）‘）;

%混合過程，要求混合濾波器可逆，怎么做到？
Nobserv=size（originalsig21）;%求行數
Nsource=Nobserv;
FilterLength=10;%10階濾波器
H?=?rand（NobservFilterLengthNsource）;???????%random?generation?of?the?mixing?filter

%卷積混合信號
mixedsig2=zeros（Nsourcesize（originalsig22））;%2行采樣點數
for?j=1:size（originalsig21）%H行數
????for?i=1:size（originalsig21）
????????mixedsig2（j:）=mixedsig2（j:）+filter（H（j:i）1originalsig2（i:））;
????end
end
figure（2）;
subplot（‘211‘）;plot（tmixedsig2（1:））;axis（[0500e-6-2020]）;title（‘觀測信號1‘）;xlabel（‘t（單位s）‘）;
subplot（‘212‘）;plot（tmixedsig2（2:））;axis（[0500e-6-2020]）;title（‘觀測信號2‘）;xlabel（‘t（單位s）‘）;
%?subplot（‘313‘）;plot（tmixedsig2（3:））;axis（[0500e-6-2020]）;title（‘觀測信號3‘）;xlabel（‘t（單位s）‘）;????????
%STFT變換
%window=hamming（512）;
nfft=512;%nfft是做FFT變換的點數，這里取nfft=窗長
%noverlap=256;%滑動步長默認是幀長的一半
%[S1F1T1]=spectrogram（mixedsig2（1:）windownoverlapnfftFs）;
%[tfr1tTF1]?=?tfrstft（mixedsig2（1:）‘tNhamming（127）0）;
%[tfr1tTF1]?=?tfrstft（mixedsig2（1:）‘tNhamming（127））;
%?X_w1=specgram（x1N_frequencyFswindowoverlap）
[S1?F1?T1]?=?s

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2019-04-12?18:58??卷積盲源分離（2?source?基于能量相關排序）\
?????文件????????5763??2016-06-30?15:38??卷積盲源分離（2?source?基于能量相關排序）\demo_2source.m
?????文件????????1886??2016-06-15?20:57??卷積盲源分離（2?source?基于能量相關排序）\istft.m
?????文件????????9218??2016-06-15?20:57??卷積盲源分離（2?source?基于能量相關排序）\jade.m
?????文件?????????370??2016-09-13?23:07??卷積盲源分離（2?source?基于能量相關排序）\Rxiang.m
?????文件??????????58??2016-06-29?10:19??卷積盲源分離（2?source?基于能量相關排序）\S.m
?????文件????????1625??2016-06-15?20:57??卷積盲源分離（2?source?基于能量相關排序）\stft.m
?????文件?????????344??2016-09-11?21:25??卷積盲源分離（2?source?基于能量相關排序）\W.m
?????文件?????????290??2016-06-15?20:57??卷積盲源分離（2?source?基于能量相關排序）\xiangsi.m