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