matlab均衡仿真报告材料

1、均衡matlab实验仿真报告学院：电子信息工程学院班级：通信1106班姓名：胡梦春学号：11211132指导教师：杨维老师时间：2014年5月27日均衡MATLAB仿真报告要求：现给出迫零均衡(ZF )、最小均方误差均衡中的最小均方 算法(LMS )的matlab 程序，理解各程序，完成以下习题。将程 序运行结果及各题目的解答写入 word中：用 matlab 分别运行“ main_zf.m ”和“ mainms.m ”(a)在程序中标注“注释”处加上注释(英文或中文)。 main_zf.mM=1500; %码元数目P=0.5; %1码概率data=2*round(rand(1,M)+P-0.

2、5)-1;%产生一列 01 码h=0.02 0.05 0.1 -0.2 1 -0.2 0.1 0.06 0.01;%注释：归一化的多径信道系数 hr=conv (data,h);%迫零均衡N=5;C=force_zero(h,N);%获得均衡器系数 Cdataout=co nv(r,C); % 注释:将经过多径信道传输的信号通过均衡器figure(1)subplot(2,2,1)plot(1:le ngth(data),data,.)title(发送信号序列)subplot(2,2,2)plot(1:le ngth(r),r,.)title(多径信号序列)subplot(2,2,3)plot(

3、1:le ngth(dataout),dataout,.)title(均衡后的信号序列)figure(2)eyediagram(r,2);title(迫零均衡前的眼图);figure(3)eyediagram(dataout,2);title(迫零均衡后的眼图);%注释：画迫零均衡后的眼图%用不同阶数的迫零均衡器均衡后的误码率，并与理想误码率曲线比较。snr_in_dB=4:11;%注释：产生不同的信噪比条件，4dB至U 11dBN=1 2 3; % 注释：均衡器抽头数为2*N+1，分别表示3、5、7阶均衡器 err_rate=zeros(length(N),length(snr_in_dB)

4、;%注释：经过均衡误码率统计err_rate仁zeros(1,le ngth(s nr_in_dB);%注释：未经过均衡误码率统计for ii=1:le ngth(N)C=force_zero(h,N(ii);forjj=1:le ngth(s nr_in_dB)SNR=10A(snr_in_dB(jj)/10); %注释：将信噪比dB值转化为线性值err=0; %注释：经过均衡的误码个数err仁0; %注释：没有经过均衡的误码个数for kk=1:10A3%注释：kk从1到1000，每种信噪比下，试验1000次x=2*round(rand(1,M)+P-0.5)-1;%注释：%产生一列 01

5、 码，1 码概率为 Px2=awg n( x,SNR,measured,li near); %注释：x序列中加入高斯白噪声x1=conv(x2,h);%注释：序列x2通过多径信道y=co nv(x1,C);%注释：将经过多径信道传输的信号通过均衡器L=(le ngth(y)-M)/2;y=y(L+1:L+M); % 注释：取L+1到L+M个数的符号y=sign(y);%注释：根据y的正负，以零为判决门限err=er叶sum(abs(x-y)/2;%根据收发序列是否相同统计误码个数% 不经过均衡L1=(le ngth(x1)-M)/2; %注释：收到序列的起始点x1仁x1(L1+1:L1+M);

6、 %注释：取L+1到L+M个数的符号y1仁sign(x11); % 注释：根据y的正负，以零为判决门限err仁err1+sum(abs(x-y11)/2;%注释：根据收发序列是否相同统计误码个数enderr_rate(ii,jj)=err/(M*10A3)%注释：统计不同阶数均衡器下的误码率err_rate1(1,jj)=err1/(M*10A3) % 注释：统计不经过均衡器的误码率end end figure;semilogy(s nr_in_dB,0.5*erfc(sqrt(1*10.A(s nr_in_dB/10),g*-);hold on;semilogy(s nr_in_dB,err

7、_rate1(1,:),g-); hold on;semilogy(s nr_in_dB,err_rate(1,:),ko);hold on;semilogy(s nr_in_dB,err_rate(2,:),go); hold on;semilogy(s nr_in_dB,err_rate(3,:),ro);title(误码率);legend(高斯信道下理想误码率特性,未进过均衡的误码率,三阶迫零均衡误码 率,五阶迫零均衡误码率,七阶迫零均衡误码率);xlabel(SNR); main_LMS.m clear;clc;Q%*变量设置*Q%N=10000;%二进制信源的长度，测误码率时使用的N

8、h = 0.6 -0.3 0.1; %多径响应序列order = 5;%C的阶数(取大于1的奇数)mu = 0.01;%步长delta=2;%延迟SNRdB = 5:15;Loops = 100,100,100,100,100,100,100,100,100,100,100;实现信源长度10A6%循环数,Q%k*Q%C = zeros(order,1);%初始化自适应系数for i = 1:length(SNRdB)% 注释：在不同信噪比下，从5dB到15dBdisp(目前仿真到第,num2str(i),轮:,SNR二,num2str(SNRdB(i),dB);%注释：不同轮数，不同信噪比条件

9、下TotalError = 0;%注释：经过均衡的误码个数，初始化为0totalerror1=0;SNR = 10A(SNRdB(i)/10);%注释：将信噪比dB值转化为线性值for JJJ = 1:Loops(i)%注释：第一轮迭代循环x=randsrc(1,N,0,1;0.5,0.5);% 注释：产生有N个数的01随机序列，0、1概 率都为0.5x1=1-x*2;%注释：将x序列中的0转化为1,1转化为-1x2 = conv(x1,h);%注释：序列x1通过多径信道x3 = awgn(x2,SNR,measured,linear);%注释：x2 序列中加入高斯白噪声% iterate循环

10、，迭代for jjj=order+1:Nrr=x3(jjj:-1:jjj-order+1);% vector of received sig nale=x1(jjj-delta)-C*rr;% calculate errorC=C+mu*e*rr;% update equalizer coefficie ntsend x6=filter(C,1,x3); %注释：通过信道传输的序列x3经过均衡器Cx7 = (1-sign(x6)/2;%注释：判决x6序列，小于0的数判为1，大于零的数判为-1%注释：统计经过均衡器注释：统计不经过均衡x8=x3pnl=dlud图4迫零均衡10-1误码率高斯信道下

11、理想谋玛率特性未逬过均衡的误码率O三阶迫零均衡误码率五阶迫零均衡误玛率O七阶迫雲均衡误码率1010-1107SNR1011图5迫零均衡最小均方误差均衡中的最小均方算法 （LMS ）中的相关图如图6、7、8、9所示:多径信号序列发送信号序列0.50061D5C0010000%60001000015000均衡后信号序列家径加噪声信号序列图6最小均方误差（LMS算法）均衡均衝前的眼團0.5im4 心 pn.+=_duu图7最小均方误差（LMS算法）均衡均衡后的眼图-2J0 600.5Time图8最小均方误差（LMS算法）均衡图9最小均方误差（LMS算法）均衡结论：对迫零均衡，均衡后的序列+1和-1值差距很大，很容易判别，眼图 也很好看，眼睛张开得很大，体现了迫零均衡是以消除码间干扰