%%?1.常規(guī)?CMA?算法

%----?16QAM信號(hào)的CMA算法-----------
%1.?----?參數(shù)設(shè)定-----------
SNR?=?30;????????????????????????%信噪比
L?=?13;???????????????????????????????%均衡器階數(shù)（階數(shù)不要太高，太高容易不穩(wěn)定，得到全局最小誤差需要的迭代次數(shù)也增加）
u_cma?=?1e-4;??????????????????%CMA迭代時(shí)候的步長(zhǎng)
u_lms?=?2e-3;???????????????????%LMS迭代時(shí)候的步長(zhǎng)
W?=?3.1;???????????????????????????%特征值分散度
N?=?10000;???????????????????????%處理的信號(hào)序列的長(zhǎng)度
R2?=?13.2;?????????????????????????%Godard常數(shù)?xd=[-3?-3;-3?-1;-3?1;-3?3;-1?-3;-1?-1;-1?1;-1?3;1?-3;1?-1;1?1;1?3;?3?-3;3?-1;3?1;3?3;]
????????????????????????????????????????%?xs=xd（:1）+xd（:2）*sqrt（-1）;??R2=mean（abs（xs）.^4）./mean（abs（xs）.^2）=13.2;?
a?=?zeros（100001）;
x?=?zeros（100001）;
y?=?zeros（N-（L-1）1）;
y1?=?zeros（N-（L-1）1）;

%產(chǎn)生隨機(jī)信道信號(hào)
b?=?rand（[10000?1]）;
c?=?rand（[10000?1]）;

for?n?=?1:10000
????if?b（n）?>?0.5
????????b（n）?=?3;
????else
????????b（n）?=?1;
????end
????if?c（n）?>?0.5
????????c（n）?=?3;
????else
????????c（n）?=?1;
????end
end
?
for?n?=?1:10000
????a（n）?=?sign（randn（1））*b（n）?+?1i*sign（randn（1））*c（n）;
end

%產(chǎn)生通過LTI信道后的響應(yīng)
h?=?zeros（13）;
for?n?=1:3
????h（n）?=?0.5*（1?+?cos（2*pi*（n-2）/W））;
end

for?n?=?3:10000
????x（n）?=?h（1）*a（n）?+?h（2）*a（n-1）?+?h（3）*a（n-2）?;?????????%通過LTI信道后的響應(yīng)
end

h1=[0.005??0.009?-0.0024?0.854?-0.218?0.049?-0.016];?%典型電話信道
h2=[1?0.5?0.25?0.125];?????????????????????????????????????????????????????%普通信道
h3=[0.0410+0.01091i?0.0495+0.01231i?0.0672+0.01701i?0.0919+0.0235?0.7920+0.12811i?0.3960+0.08711i?0.2715+0.04981i?0.2291+0.04141i?0.1287+0.01541i?0.1032+0.01191i];
h4=[0.3132?-0.104?0.8908?0.3134];???%[0.3132??-?o.??104??0.8908??0.3134JC
hn=h4;
x=conv（hna）;
x=conv（[0.005??0.009?-0.0024?0.854?-0.218?0.049?-0.016]a）;
x?=?awgn（xSNR）;?????????????????????????????????????????????????%加入白噪聲，SNR在文件開始出設(shè)定

%盲均衡算法
for?k=1:N-L+1
????X（k:）?=?x（k:k+L-1）?;
end

w（11:L）?=?0;
w（1（L+1）/2）?=?1;
aL=0.99996;
u_cma?=?5e-5;???
for?k?=?1:N-（L-1）
????y（k）?=?X（k:）*w（k:）‘;???
????e（k）?=?y（k）*（R2-abs（y（k））^2）;
????w（k+1:）?=??w（k:）+u_cma*e（k）‘?*?X（k:）;
????%{
%對(duì)CMA算法計(jì)算得到的點(diǎn)進(jìn)行判定
????if?（abs（y（k）?-?（3+3i））?????????y1（k）?=?3?+?3i;
????elseif?（abs（y（k）?-?（3-3i））?????????y1（k）?=?3?-?3i;
????elseif?（abs（y（k）?-?（-3+3i））?????????y1（k）?=?-3?+?3i;
????elseif?（abs（y（k）?-?（-3-3i））?????????y1（k）?=?-3?-?3i;
????????
????elseif?（abs（y（k）?-?（1+3i））?????????y1（k）?=?1+3i;
????elseif?（abs（y（k）?-?（1-3i））?????????y1（k）?=?1-3i;
????elseif?（abs（y（k）?-?（-1+3i））?????????y1（k）?=?-1+3i;
????elseif?（abs（y（k）?-?（-1-3i））?????????y1（k）?=?-1-3i;
????????
????elseif?（abs（y（k）?-?（3+i））?????????y1（k）?=?3+i;
????elseif?（abs（y（k）?-?（3-i））?????????y1（k）?=?3-i;
????elseif?（abs（y（k）?-?（-3+i））?????????y1（k）?=?-3+i;
????elseif?（abs（y（k）?-?（-3-i））?????????y1（k）?=?-3-i;
????????
????elseif?（abs（y（k）?-?（1+i））?????????y1（k）?=?1+i;
????elseif?（abs（y（k）?-?（1-i））?????????y1（k）?=?1-i;
????elseif?（abs（y（k）?-?（-1+i））?????????y1（k）?=?-1+i;
????elseif?（abs（y（k）?-?（-1-i））?????????y1（k）?=?-1-i;
????????
????else
??????

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????13642193??2017-12-20?11:13??（全文）?自適應(yīng)盲均衡技術(shù)?郭業(yè)才著.pdf
?????目錄???????????0??2018-01-23?09:27??Matlab代碼\
?????文件?????????255??2017-12-25?10:07??Matlab代碼\averageC2.mat
?????文件?????????341??2017-12-27?09:41??Matlab代碼\averageC4.mat
?????文件???????12546??2017-12-22?14:59??Matlab代碼\CMA_DDLMS.m
?????文件????????4411??2017-12-20?15:51??Matlab代碼\Rayleigh_16QAM_CMA.m
?????文件????????2025??2017-12-20?15:48??Matlab代碼\Rayleigh_4QAM_CMA.m
?????目錄???????????0??2018-05-10?19:21??新建文件夾\
?????目錄???????????0??2018-09-11?10:59??新建文件夾（1）\