您现在的位置:首页 >> 技术文章 >> 图像语音处理 >> 内容

FPGA代做-光纤通信系统中的克服PMD效应的电域均衡技术

时间:2019-2-26 23:07:26 点击:

  核心提示:FPGA代做-光纤通信系统中的克服PMD效应的电域均衡技术...

主要是研究光纤通信系统中的克服PMD效应的电域均衡技术,构建电域均衡器的模型并进行仿真实验。找到补偿PMD效应的较理想的电域均衡器结构。

%%%%***********************************************************************

%%%%%%%%基于lms算法的判决反馈均衡器

%%%%%%%%前溃阶数取5,后溃阶数取3

%%%%%%%

clear; close all;

%产生数字信号

p = .5;

x = binornd(1 ,p, 1, 5000);

%模拟经过信道传输,加入加噪声

x1=x;

x1(2:5000)=0.3*x(2:5000)+0.7*x(1:4999);


v = 0.15 * rand(1,5000);



y = x1+v;

%%%取前1000个点,均衡器工作于训练阶段

%设置均衡器前溃,后溃参数

mu1 = 0.08;

mu2 = 0.000008;

%前溃,后溃抽头权重初值

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) = x(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); x(n-1); x(n-2); x(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); x(n-1); x(n-2); x(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);

    c3(n+1) = c3(n) + mu1 * e(n) * y(n+3);

    c2(n+1) = c2(n) + mu1 * e(n) * y(n+2);

    c1(n+1) = c1(n) + mu1 * e(n)* y(n+1);

    c0(n+1) = c0(n) + mu1 * e(n) * y(n);

    d1(n+1) = d1(n) + mu2 * e(n) * x(n-1);

    d2(n+1) = d2(n) + mu2 * e(n)* x(n-2);

    d3(n+1) = d3(n) + mu2 * e(n) * x(n-3);

end

% figure 

% subplot(411), plot(x(1:1000)); title('发射信号')

% subplot(412), plot(y(1:1000)); title('接受到信号')

% subplot(413), plot(xes); title('均衡后量化前信号')

% subplot(414), plot(xes2); title('均衡后量化后信号')

% %计算均衡及量化误差

% for i = 1:997

%     error1(i) = (x(i) - xes2(i))^2;

%     error2(i) = (x(i) - xes(i))^2;

% end

% figure,

% subplot(211), plot(error2); title('均衡后量化前误差信号');

% subplot(212), plot(error1); title('均衡后量化后误差信号');

%%%%%%%均衡器工作在判决阶段

%设置均衡器前溃,后溃参数

mu1 = 0.0000002;

mu2 = 0.000002;

%%%均衡算法

for n=997:4996

    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);

    c3(n+1) = c3(n) + mu1 * e(n)* y(n+3);

    c2(n+1) = c2(n) + mu1 * e(n)* y(n+2);

    c1(n+1) = c1(n) + mu1 * e(n)* y(n+1);

    c0(n+1) = c0(n) + mu1 * e(n)* y(n);

    d1(n+1) = d1(n) + mu2 * e(n)* xes(n-1);

    d2(n+1) = d2(n) + mu2 * e(n)* xes(n-2);

    d3(n+1) = d3(n) + mu2 * e(n)* xes(n-3);

    

end


% figure 

% subplot(411), plot(x(1001: 4990)); title('发射信号')

% subplot(412), plot(y(1001:4990)); title('接收到信号')

% subplot(413), plot(xes(1002:4991)); title('均衡后量化前信号')

% subplot(414), plot(xes2(1002:4991)); title('量化后信号')

%     

% %计算均衡及量化误差

% for i = 998:4997

%     error1(i) = (x(i-1) - xes2(i))^2;

%     error2(i) = (x(i-1) - xes(i))^2;

% end

% figure,

% subplot(211), plot(error1); title('量化后误差信号');

% subplot(212), plot(error2); title('量化前误差信号');

%画眼图

eyediagram(x(997:length(x)),4,0.5);

title('发射信号眼图')

eyediagram(y(997:length(y)),4,0.5);

title('未均衡前眼图')

eyediagram(xes(998:length(xes2)),4,0.5);

title('5阶FFE+3阶DFE均衡后眼图');xlabel('时间');ylabel('幅度');

% %计算误码率

for n=1:5000

 if y(n) > 0.5

       y(n) = 1;

   elseif y(n) < 0.5

       y(n) = 0;

    end

end

xes2(4998:5000)=y(4997:4999);

[h1,h2]=biterr(x(1000:5000),y(1000:5000))

[h3,h4]=biterr(x(999:4999),xes2(1000:5000))


 

联系:highspeedlogic

QQ :1224848052

微信:HuangL1121

邮箱:1224848052@qq.com

网站:http://www.mat7lab.com/

网站:http://www.hslogic.com/

微信扫一扫:

作者:电域均衡 来源:电域均衡
本站最新成功开发工程项目案例
相关文章
  • 没有相关文章
相关评论
发表我的评论
  • 大名:
  • 内容:
  • matlab代做|matlab专业代做|matlab淘宝代做(www.hslogic.com) © 2019 版权所有 All Rights Reserved.
  • Email:highspeed_logic@163.com 站长QQ: 1224848052

    专业代做/代写/承接、MATLAB、SIMULINK、FPGA项目、博士/硕士/本科毕业设计、课题设计、论文,毕业论文,Coursework、Eassy、Assignment