盲源分离算法仿真分析系统.docVIP

盲源分离算法仿真分析系统 程序界面： 使用帮助： 信号种类： MASK：2ASK，4ASK，8ASK MPSK：2PSK、4PSK、8PSK MFSK：2FSK、4FSK、8FSK 分连续相位CPFSK和离散相位DPFSK两种 MQAM：16QAM、32QAM、64QAM、128QAM OFDM：子载波数可任意选定，映射方式有BPSK、QPSK、4QAM、16QAM、32QAM、64QAM Sign：符号函数 Noise：噪声 算法种类： 文件说明： runme.m 主程序，可运行 coef.m 计算相似系数 CrossTalkError.m 计算串音误差 FastICA.m FastICA算法 MASK.m 产生ASK信号 MFSK.m 产生FSK信号 MPSK.m 产生PSK信号 MQAM.m 产生QAM信号 NG.m 自然梯度算法 NJD.m 非正交联合对角化算法 OFDM.m 产生OFDM信号 proj_NG.m 投影自然梯度算法 runme.fig 界面图像文件 sir.m 计算信干比 SOBI.m SOBI算法 sub_njd.m 非正交联合对角化子函数，由NJD.m调用 作者： HYPERLINK mailto  西安电子科技大学 MFSK.m function y=MFSK(x,fc,df,fs,rc,M,type) %----MFSK调制 %参数说明 % y 输出序列 % x 输入序列,为0，1序列，可以用函数randint（1，N，2）产生 % fc 载波频率 % df 频偏 % fs 采样频率 % M 调制方式，M=2、4、8 % type type=1为相位连续MFSK调制（CPFSK），type=2为相位不连续MFSK调制（DPFSK） N=length(x); Ns=fs/fc; ini_detfy=rand()*2*pi; %随机初始相位，根据需要可以手动设置为0 switch M case 2, n=1:floor(N*Ns*fc/rc)-1; t=n/(Ns*fc); m=x(ceil(n*rc/(Ns*fc))); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% if type==1 %相位连续调制 detfy=0; %频率变化时相位改变值 temp=0; %累计相位改变值 s(1)=cos(2*pi*(fc+m(1)*df)*t(1)); %初始化第一个点 for k=2:length(n) if m(k)-m(k-1)~=0 temp=temp+detfy; %计算累计相位改变值 detfy=2*pi*(m(k)-m(k-1))*df*t(k-1); %计算瞬时相位改变值 s(k)=cos(2*pi*(fc+m(k)*df)*t(k)+ini_detfy-temp-detfy);%计算信号值 else s(k)=cos(2*pi*(fc+m(k)*df)*t(k)+ini_detfy-temp-detfy); %频率不改变时 end end else s=cos(2*pi*(fc+m*df).*t+ini_detfy); %相位不连续调制 end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% case 4, if mod(N,2)==1 N=