matlab-基于nlms算法的判决反馈均衡器-源码

%基于nlms算法的判决反馈均衡器 %前溃阶数取5，后溃阶数取3 clc; clear; close all; warning off; addpath(genpath(pwd)); %**************************************************************************** %更多关于matlab和fpga的搜索“fpga和matlab”的CSDN博客: %matlab/FPGA项目开发合作 %https://blog.csdn.net/ccsss22?type=blog %**************************************************************************** %产生数字信号 p = .5; x = binornd(1, p, 1, 5000); %模拟经过信道传输，加入加噪声 for i=1:5000 z(2*i)=x(i); z(2*i-1)=x(i); end x1=z; x1(2:10000)=0.3*z(2:10000)+0.7*z(1:9999); v = 0.0316 * randn(1,10000); y = x1+ v; % y=x1; %%%取前1000个点，均衡器工作于训练阶段 %设置均衡器前溃，后溃参数 mu1 = 0.08; mu2 = 0.00008; %前溃，后溃抽头权重初值 c=1; param = 0; c0(5) = param; c1(5) = param; c2(5) = param; c3(5) = param; c4(5) = param; d1(5) = param; d2(5) = param; d3(5) = param; for n=5:997 e(n) = z(n) - [c4(n) c3(n) c2(n) c1(n) c0(n) -d1(n) -d2(n) -d3(n)] * [y(n+4); y(n+3); y(n+2); y(n+1); y(n); z(n-1); z(n-2); z(n-3)]; xes(n) = [c4(n) c3(n) c2(n) c1(n) c0(n) -d1(n) -d2(n) -d3(n)] * [y(n+4); y(n+3); y(n+2); y(n+1); y(n); z(n-1); z(n-2); z(n-3)]; if xes(n) > 0.5 xes2(n) = 1; elseif xes(n) < 0.5 xes2(n) = 0; end c4(n+1) = c4(n) + mu1 * e(n) * y(n+4)/(c+(y(n+4)'*y(n+4))); c3(n+1) = c3(n) + mu1 * e(n) * y(n+3)/(c+(y(n+3)'*y(n+3))); c2(n+1) = c2(n) + mu1 * e(n) * y(n+2)/(c+(y(n+2)'*y(n+2))); c1(n+1) = c1(n) + mu1 * e(n) * y(n+1)/(c+(y(n+1)'*y(n+1))); c0(n+1) = c0(n) + mu1 * e(n) * y(n)/(c+(y(n)'*y(n))); d1(n+1) = d1(n) + mu2 * e(n) * x(n-1)/(c+(z(n-1)'*z(n-1))); d2(n+1) = d2(n) + mu2 * e(n) * x(n-2)/(c+(z(n-2)'*z(n-2))); d3(n+1) = d3(n) + mu2 * e(n) * x(n-3)/(c+(z(n-3)'*z(n-3))); end mu1 = 0.000001; mu2 = 0.00002; for n=997:5000 e(n) = xes2(n) - [c4(n) c3(n) c2(n) c1(n) c0(n) -d1(n) -d2(n) -d3(n)] * [y(n+4); y(n+3); y(n+2); y(n+1); y(n); xes2(n-1); xes2(n-2); xes2(n-3)]; xes(n+1) = [c4(n) c3(n) c2(n) c1(n) c0(n) -d1(n) -d2(n) -d3(n)] * [y(n+4); y(n+3);y(n+2); y(n+1); y(n); xes2(n-1); xes2(n-2); xes2(n-3)]; if xes(n+1) > 0.5 xes2(n+1) = 1; elseif xes(n+1) < 0.5 xes2(n+1) = 0; end c4(n+1) = c4(n) + mu1 * e(n) * y(n+4)/(c+(y(n+4)'*y(n+4))); c3(n+1) = c3(n) + mu1 * e(n) * y(n+3)/(c+(y(n+3)'*y(n+3))); c2(n+1) = c2(n) + mu1 * e(n) * y(n+2)/(c+(y(n+2)'*y(n+2))); c1(n+1) = c1(n) + mu1 * e(n) * y(n+1)/(c+(y(n+1)'*y(n+1))); c0(n+1) = c0(n) + mu1 * e(n) * y(n)/(c+(y(n)'*y(n))); d1(n+1) = d1(n) + mu2 * e(n) * xes(n-1)/(c+(xes(n-1)'*xes(n-1))); d2(n+1) = d2(n) + mu2 * e(n) * xes(n-2)/(c+(xes(n-2)'*xes(n-2))); d3(n+1) = d3(n) + mu2 * e(n) * xes(n-3)/(c+(xes(n-3)'*xes(n-3))); end eyediagram(x(1000:length(x)),8); title('发射信号眼图') eyediagram(y(1000:length(y)),8); title('未均衡前眼图') eyediagram(xes(1000:length(xes2)),8); title('均衡后眼图');