mcma.m

基于matlab的MCMA算法仿真程序

% 盲信道均衡MCMA算法
clear all;
clc
N=12000;%数据点总数
SNRdB=25;%信噪比
Lf=11;%滤波器阶数=Lf + 1
% Lh=5;% 信道阶数
Dy=round((Lf+1)/2) ;
% h = [0.9344-1.0311i 2.2483-1.6682i -1.0780 -0.2138i -0.6742+0.8835i 0.3755-0.2771i 0.0886+0.0628i];
h=[1 0 0.5];
% 信道冲激响应
M=256; % 产生8PSK信号
x=randint(N,1,M);
s=qammod(x,M)';
s=s/std(s);
sx=s;
% sx=s.*exp(-j*pi/6);
% sx=awgn(sx,SNRdB,'measured');
% 
% x=filter(h,1,sx) ; % 信号通过信道
x=sx;
x=x/std(x) ; % 信号功率归一化


Lp=10000;%迭代次数

% Lp=N-Lf;%迭代次数
X=zeros(Lf+1,1);
e=zeros(1,Lp);
w=zeros(Lf+1,1);w(Dy)=1; % 初始化

Rr=mean(abs(real(s)).^4)/mean(abs(real(s)).^2)
Ri=mean(abs(imag(s)).^4)/mean(abs(imag(s)).^2)
mu=0.002; % 学习速率
for i=1:Lp
    X=x(i+Lf:-1:i).';
    yk=w'*X;% 滤波器输出
    er(i)=real(yk)*(real(yk)^2-Rr); %产生误差信号
    ei(i)=imag(yk)*(imag(yk)^2-Ri); %产生误差信号
    e(i)=er(i)+j*ei(i);
    w =w-mu*e(i)'*X; %更新
end
for i=1:Lp
    sb(i)=w'*x(i+Lf:-1:i).'; % 均衡接收信号
end ;
% 
% % c=subclust([real(sb)' imag(sb)'],0.10)
subplot(1,2,1)
plot (sx/std(sx),'.') ;
grid ,title ('存在相位偏转信号星座图，SNR=25dB') ; xlabel ('Real') ,ylabel ('Image')
subplot(1,2,2)
plot (x ,'.') ;
grid , title ('接收信号星座图') ; xlabel ('Real') , ylabel ('Image')
figure(3),
% plot (real(sb),imag(sb) ,'.',c(:,1),c(:,2),'rs') ;
plot (real(sb),imag(sb) ,'.') ;
grid , title ( '均衡后的接收信号图') , xlabel ( 'Re2al') , ylabel ('Image')