数字通信原理程序v1.doc

数字通信原理实验 参考书： 《通信原理MATLAB仿真教程》，人民邮电出版社 第六章 数字基带系统 一、眼图 % 升余弦滚降系统眼图% 数字基带信号波形及其眼图% 文件名：eye.mTs=1;N=15;eye_num=6;a=1;N_data=1000;dt=Ts/N;t=-3*Ts:dt:3*Ts;% 产生双极性数字信号d=sign(randn(1,N_data));dd=sigexpand(d,N);? %% 基带系统冲击响应（升余弦）ht=sinc(t/Ts).*(cos(a*pi*t/Ts))./(1-4*a^2*t.^2/Ts^2+eps);st=conv(dd,ht);tt=-3*Ts:dt:(N_data+3)*N*dt-dt;subplot(211)plot(tt,st);axis([0 20 -1.2 1.2]);xlabel(t/Ts);ylabel(基带信号);subplot(212)% 画眼图ss=zeros(1,eye_num*N);ttt=0:dt:eye_num*N*dt-dt;for k=3:50??? ss=st(k*N+1:(k+eye_num)*N);??? drawnow;??? plot(ttt,ss);??? hold on;end;xlabel(t/Ts);ylabel(基带信号眼图); --------------------------------------- % 将输入的序列扩成间隔为N-1个0的序列function[out]=sigexpand(d,M) N=length(d);out=zeros(M,N);out(1,:)=d;out=reshape(out,1,M*N); 双极性信号的功率谱。f=0:0.01:5;Ts=1;x=f*Ts;y=sin(pi*x);y=y./(pi*x);y(1)=1;dnrz=y.*y;dnrz=Ts*dnrz;? % 计算双极性非归零码的功率谱y=sin(pi*x/2);y=y./(pi*x/2);y(1)=1;drz=y.*y;drz=Ts*drz/4;? % 计算双极性归零码的功率谱plot(x,dnrz,:,x,drz,-);title(双极性非归零码DNRZ和归零码DRZ的功率谱);xlabel(f);ylabel(双极性(P=1/2));legend(DNRZ,DRZ); 比较单极性非归零码NRZ和双极性非极性归零码NRZ的输入信噪比与误比特率关系特性， % 并计算当要求基带传输系统的误比特率为10^(-6)时所需要的信噪比sn=0.1:0.01:100; % 定义信噪比序列snlg=20*log10(sn); % 将信噪比转化为dB表示sdouble=sqrt(sn);ssingle=sqrt(sn/2);bdouble=erfc(sdouble)/sqrt(2); % 求双极性的误比特率序列bsingle=erfc(ssingle)/sqrt(2); % 求单极性的误比特率序列semilogy(snlg,bdouble);hold;? % 保持住上一曲线semilogy(snlg,bsingle,--);hold;legend(bdouble,bsingle);axis([-20 30 0.0000001 1]);title(单、双极性非归零码NRZ);xlabel(输入信噪比);ylabel(误比特率); % 通过观察曲线，大致判断出双极性小于给定误比特率的信噪比位置i=1100;while(ilength(sn))??? if(bdouble(i)10^(-6))??????? bsn=snlg(i);? % 记下符合条件的信噪比??????? i=length(sn)+1;??? end;??? i=i+1;end;% 通过观察曲线，大致判断出单极性小于给定误比特率的信噪比位置i=2300;while(ilength(sn))??? if(bsingle(i)10^(-6))??????? ssn=snlg(i);? % 符合所需的条件时，记录下此时对应的信噪比值??????? i=length(sn)+1;??? end;??? i=i+1;end; disp(双极性NRZ码所需的信噪比为：(dB));bsn?? % 显示双极性的信噪比disp(单极性NRZ码所需的信噪比为：(dB));ssn?? % 显示单极性的信噪比 升余弦滚降特性（数字基带信号?码间串扰） % 数字基带信号传输? 码间串扰 升余弦滚降系统的频谱及其时域波形