function TEST_Convergence_LMS(L1,delta1,L2,delta2,L3,delta3,interations,trials)
% Chuong trinh mo phong khao sat su hoi tu cua giai thuat LMS
d = interations;
er1 =[];er2=[];er3=[];
for n = 1:trials
M = 4;
SigLen = 5000;
EbNodB = 20;
% Dieu che M_PSK
Ii = randi([0 M-1],1,SigLen);
theta = (2*pi/M)*Ii+pi/M;
Si = exp(1i*theta);
% Tinh toan cac thong so lien quan
EbNo = 10.^(EbNodB./10);
Es = norm(Si)^2/SigLen;
Eb = Es/log2(M);
No = Eb./EbNo;
% delta =(2/((2*L + 1)*(No + 2*Es))/10);
% Bo can bang thich nghi su dung giai thuat LMS
% L1 = lengthEq1;
% delta1 = stepsize1;
eq1 = lineareq(L1,lms(delta1));
eq1.SigConst = exp(1i*((2*pi/M)*(0:M-1)+pi/M));
eq1.ResetBeforeFiltering=1;
% L2 = lengthEq2;
% delta2 = stepsize2;
eq2 = lineareq(L2,lms(delta2));
eq2.SigConst = exp(1i*((2*pi/M)*(0:M-1)+pi/M));
eq2.ResetBeforeFiltering=1;
% L3 = lengthEq3;
% delta3 = stepsize3;
eq3 = lineareq(L3,lms(delta3));
eq3.SigConst = exp(1i*((2*pi/M)*(0:M-1)+pi/M));
eq3.ResetBeforeFiltering=1;
% Mo hinh tuong duong dap ung cua kenh truyen theo he so cua bo loc FIR
chan = [.986; .845; .237; .123+.31i];
% AWGN
noise = sqrt(No./2).*(randn(1,SigLen)+1i*randn(1,SigLen));
% Tin hieu tai ngo vao cua bo can bang
ychan = filter(chan,1,Si);yin = ychan + noise;
% Tin hieu tai ngo ra cua bo can bang
[Sest1,Sd1,err1] = equalize(eq1,yin,Si);
[Sest2,Sd2,err2] = equalize(eq2,yin,Si);
[Sest3,Sd3,err3] = equalize(eq3,yin,Si);
er1 = [er1;err1];er2 = [er2;err2];er3 = [er3;err3];
end
% Tinh toan Mean Square Error
er1 = mean(er1);er2 = mean(er2);er3 = mean(er3);
e1 = abs(er1);MSE1 = 1/2*(e1.^2);MSE1 = MSE1(1:d);
e2 = abs(er2);MSE2 = 1/2*(e2.^2);MSE2 = MSE2(1:d);
e3 = abs(er3);MSE3 = 1/2*(e3.^2);MSE3 = MSE3(1:d);
maxEr = roundn(max([MSE1 MSE2 MSE3]),-1);
eravg1 = mean(MSE1);eravg2 = mean(MSE2);eravg3 = mean(MSE3);
N = 1:d;
% Plot
figure();
subplot(3,1,1);
plot(N,MSE1,'b',N,eravg1,'k');
axis([0 d 0 maxEr]);
title('Khao sat su hoi tu cua bo can bang voi giai thuat LMS');
legend(['L = ', num2str(L1),',','step-size = ',num2str(delta1)]);
subplot(3,1,2);
plot(N,MSE2,'g',N,eravg2,'k');
axis([0 d 0 maxEr]);
ylabel('MSE');
legend(['L = ', num2str(L2),',','step-size = ',num2str(delta2)]);
subplot(3,1,3);
plot(N,MSE3,'y',N,eravg3,'k');
xlabel('So lan lap di lap lai');
legend(['L = ', num2str(L3),',','step-size = ',num2str(delta3)]);
axis([0 d 0 maxEr]);
end
没有合适的资源?快使用搜索试试~ 我知道了~
TEST_Convergence_LMS.zip_LMS CONVERGENCE_convergence_lms
共1个文件
m:1个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 154 浏览量
2022-09-14
21:55:14
上传
评论
收藏 1KB ZIP 举报
温馨提示
TEST Convergence LMS
资源推荐
资源详情
资源评论
收起资源包目录
TEST_Convergence_LMS.zip (1个子文件)
TEST_Convergence_LMS
TEST_Convergence_LMS.m 2KB
共 1 条
- 1
资源评论
weixin_42653672
- 粉丝: 94
- 资源: 1万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功