基本信息

---

下载

   友情提示：（无需注册或充值，赞助后即可获取资源下载链接）

---

     嘿，亲！知识可是无价之宝呢，但咱这精心整理的资料也耗费了不少心血呀。小小地破费一下，绝对物超所值哦！如有下载和支付问题，请联系我们QQ(微信同号)：813200300

本次赞助数额为： 2 元　

支付宝支付

微信支付

   源码介绍

---

clear all
clc
N=5000;
M=16;
m=2000;
L=11;
n=500;
d=0.65;
SNR=25;
uc=10^-4;
ud=6*10^-5;
w=[0 0 0 0 0 1 0 0 0 0 0 ];
h=[0.005 0.009 -0.024 0.854 -0.218 0.049 -0.016 ];
mmse=zeros(1,N-L 1);
for j=1:m
a=randint(1,N,M);
s=qammod(a,M);

a1=abs(s).^4;
a2=abs(s).^2;
R=mean(a1./a2);
s1=filter(h,1,s);
x=awgn(s1,SNR,'measured');
for i=1:N-L 1

  y=x(i L-1:-1:i);
  ds(i)=w*conj(y');
  e(i)=abs(ds(i))^2-R;
 w=w-uc*e(i)*ds(i)*conj(y);
 if (i>n) & (abs(ds(i)-x(i))<d)
     RD=abs(ds(i))^2;
      e(i)=abs(ds(i))^2-RD;
      w=w-ud*e(i)*ds(i)*conj(y);
 end
 mse(i)=e(i)^2;
end
mmse=mmse mse;
end
mmse=mmse/m;

plot([1:N-L 1],mmse)
xlabel('iteration numbers')
ylabel('MQAM CMA MMA mse')
hold on
scatterplot(s)
title('original signal');
grid on
hold on
scatterplot(x)
title('received signal');
grid on
hold on
scatterplot(ds)
title('CMA MMAdetected signal');
grid on
hold on