请问各位高手：
（1）以下是本人编写的CMA算法的程序，从仿真的结果得到，该算法滤波后的星座图（眼图）可以张开，但是MSE曲线却不收敛，而是形成一条直线的样子。这是怎么回事？？该如何解决？？
（2）恒模算法的权值迭代中，有的量需要取共轭复数，而再实际仿真的时候如果取共轭复数则算法就发散了，这是为什么？？？对于LMS算法也是这种情况。如何解释这种理论与仿真的差异呢？？？
（3）请大家针对我下面的程序给予指教。
（4）最后，如果那位专家有CMA仿真的程序，可以参考一下吗 ？？？？
谢谢！！！！！！！！！！！！！！！！！！！！！！！！！！！！！
N=6000; %信息序列的长度
K=11; %11个抽头系数
%实数信道
%actual_isi=[0.05 -0.063 0.088 -0.126 -0.25 0.9047 0.25 0 0.126 0.038 0.088];%有ISI的信道响应
%复数信道
%actual_isi=[0.0410+0.0109j 0.0495+0.0123j 0.0672+0.0170j 0.0919+0.0235j 0.7920+0.1281j 0.3960+0.0871j 0.2715+0.0498j 0.2291+0.0414j 0.1287+0.0154j 0.1032+0.0119j];
actual_isi=[0.005 0.009 -0.024 0.854 -0.218 0.049 -0.016];
Num_of_realizations=50; %学习次数
mse_av=zeros(1,N-K+1);
snr=30;
delta_k=0.0001;
sigpower=1;
d=1; %QAM的欧式空间距离就是1了
M=16; %QAM的进制数
mapping=[-3*d 3*d;-d 3*d; d 3*d; 3*d 3*d;-3*d d;-d d; d d; 3*d d;
-3*d -d;-d -d; d -d; 3*d -d;-3*d -3*d;-d -3*d; d -3*d; 3*d -3*d];
A=mapping;
A=A(:,1)+A(:,2)*j;
R2=mean((abs(A)).^4)/mean((abs(A)).^2); %输入信号的自相关
for jt=1:Num_of_realizations
%信源产生
for i=1:N
temp=rand; %rand是均匀随机数产生命令
dsource(i)=1+floor(M*temp);
end;
%信息到信号点的映射
for i=1:N
qam_sig(i,:)=mapping(dsource(i),:);
end;
%发送复数信号
r=qam_sig;
plot(r(:,1),r(:,2),'*'); %看星座图
rr=r(:,1)+j*r(:,2);
if(jt==1) end;
%序列通过有ISI的信道滤波器
ya1=filter(actual_isi,1,rr);
ya2=ya1';
yout= awgn(ya2,snr,sigpower);
estimated_c(1:K)=0;
estimated_c((K+1)/2)=1; %均衡器抽头系数的初值
for k=1:N-K+1 %迭代次数，490次。
y_k=yout(k:k+K-1);%;*exp(-j*an_w);
z_k=estimated_c*y_k'; %均衡器的输出信号
yt(k)=z_k;
gi_k=real(z_k)*(1+R2-real(z_k)^2); %误差信号，即期望（检测器输出）－均衡器的输出
gq_k=imag(z_k)*(1+R2-imag(z_k)^2);
g_k=gi_k+j*gq_k;
e_k=g_k-z_k;
mse(k)=abs(e_k)^2/2; %误差的平方
estimated_c=estimated_c+delta_k*y_k*e_k; %LMS算法的核心表达式。
end;
if(jt==1) end;
mse_av=mse_av+mse; %均方误差的累加。
end;
mse_av=mse_av/Num_of_realizations; %取平均
figure(1);
plot(mse_av);
figure(3);
plot(real(yt),imag(yt),'*'); %保留500个点数
