发布时间 : 星期日 文章北邮scilab - 通信原理软件实验报告 - 图文更新完毕开始阅读041550835ef7ba0d4a733bf7
四、思考题
（1）观察HDB3码的频谱图，分析其与原码频谱的不同。
答：通过观察可得，HDB3码的频谱图频率集中在中频部分,高频和低频分量很小，而原码频谱集中在低频部分。
（2）HDB3码译码如何仿真实现？请给出SCICOS模块连接的工程图，以及相关波形图。 （3）线路码在数字通信系统中起什么作用？除了HDB3码，还有哪些线路码？试列出各自优缺点以及适用的场合。 答：在实际数字通信中，经常需要在数字通信设备之间通过同轴电缆或其他有线传输媒介来传输数字基带信号，由于在基带信道传输时，不同传输媒介具有不同的传输特性，所以需要使用不同的接口线路码型。除了HDB3码外还有AMI码，CMI码，数字双相码等线路码，其中AMI码其功率谱无离散的直流分量，低频和高频分量较小，能减小码间干扰，具有检错能力以及能提取出时钟，缺点是在连“0”时时钟提取困难。而HDB3码正好改善了AMI码的缺点，不会出现连“0”的情况，因此提取时钟简单。CMI码功率谱不仅含连续谱，还含有离散的时钟分量及其奇次谐波分量，无离散直流分量。数字分相码在收端利用简单的非线性变换后提取时钟方便。但是它所提取的时钟频率是符号速率的2倍，再由它分频得到的定时信号，必定存在相位的不确定问题。
综合类实验一 数字基带系统的仿真
一、实验目的
（1）了解在理想限带及加性白高斯噪声信道条件下数字基带系统的基本原理和设计方法，
第33页
完成在仿真平台上的系统搭建与仿真。
（2）进一步熟悉SCICOS下的复杂系统设计。
（3）运用工具库现有的通信工具模块搭建系统，利用其功能全面且封装性强的特点，针对数字基带系统进行横向功能分解，使系统设计更加精炼、准确。
（4）深入学习数字基带系统核心传输节点的性能，并掌握眼图示波器的使用方法，观察接收滤波器输出的眼图和功率谱密度，判断系统传输的正确性和精准性，调试以达到最佳传输效果。
二、实验原理
若使得在接收端抽样时刻码间干扰为零，则系统的合成传递函数必须满足以下条件：
则接收到的确定信号的频谱仅取决与发送滤波器的特性，所以接收滤波器的应与发送滤波器共轭匹配，即：
这样，在理想限带信道情况下，既要使接收端抽样时刻的抽样值无码间干扰，又要使得在抽样时刻抽样值的信噪比最大，则
综上所述，数字PAM 信号通过限带信道、并受到加性噪声干扰的情况下，在限带信道为理想低通条件下，最佳基带传输的发送及接收滤波的设计原则为：总的收发系统的传递函数要符合无码间干扰基带传输的升余弦特性；且又要考虑在抽样时刻信噪比最大的收、发滤波共轭匹配的条件。
第34页
可得无码间干扰的条件下，其系统框图如图5.150：
三、实验过程
(1)编程实现
模块连接图：
在Diagram--Context中进行如下内容设置： function X=t2f(x) H=fft（x）;
X=[H(mtlb_imp(mtlb_a(N/2,1),N)),H(mtlb_imp(1,N/2))]*dt; endfunction
function x=f2t(X)
S=[X(mtlb_imp(mtlb_a(N/2,1),N)),X(mtlb_inp(1,N/2))]; x=ifft(S)/dt; endfunction
dt=0.01; L=32; M=16; N=L*M Ts=L*dt;
第35页
Rb=1/Ts; df=1/(N*dt); T=N*dt; Bs=N*df/2; alpha=0.5;
t=linspace(-T/2,T/2,N);
f=linspace(-Bs,Bs,N)+%eps;
hr1=sin(%pi*t/Ts)./(%pi*t/Ts);
hr2=cos(((alpha*%pi)*t)/Ts)./(1-(((2*alpha)*t)/Ts).^2); hr=hr1.*hr2; HR=abs(t2f(hr)); GT=sqrt(HR); GR=GT;
Sending filter超级模块中的Scifunc模块的代码： function X=t2f(x) H=fft(x);
X=[H(mtlb_imp(mtlb_a(N/2,1),N)),H(mtlb_imp(1,N/2))]*dt; endfunction
function x=f2t(X)
S=[X(mtlb_imp(mtlb_a(N/2,1),N)),X(mtlb_imp(1,N/2))]; x=ifft(S)/dt; endfunction
u=u1’; S=t2f(u); S1=S.*GT; y=real(f2t(S1)); y1=y’
receiving filter超级模块中的Scifunc模块的代码： function X=t2f(x) H=fft(x);
X=[H(mtlb_imp(mtlb_a(N/2,1),N)),H(mtlb_imp(1,N/2))]*dt; endfunction
function x=f2t(X)
S=[X(mtlb_imp(mtlb_a(N/2,1),N)),X(mtlb_imp(1,N/2))]; x=ifft(S)/dt; endfunction
u=u1’; SR=t2f(u);
第36页