. .
对于仿真系统，若x是时域取样值矢量，X是对应的傅氏变换，那么x的功率谱便为
P??X.*conj?X??/T。针对随机过程X(t)，其平均功率谱密度定义为各样本功率谱密度
的数学期望
1Px(f)?K?P(f)
kk?1K3、作出仿真图
由于需要作出的图形较多，且图形间需要对比，故采用了两种视图进行绘图，一是各个占空比的RZ码波形图和其功率谱进行横向对比，二是分别作出各占空比下的单双极性归零码波形，以便于观察。
另外，各个占空比的RZ码波形和其频谱变换后的结果使用多行的矩阵进行存储，方便最后作图，因而代码显得有些冗余。可改用定义函数，输入参数的方式给出不同占空比下的计算与绘图。
四、程序框图
对于单极性归零码：
产生M个0、1等概随机码 产生各占空比单极性归零码波形 单极性归零码的功率谱密度 作图 对于双极性归零码：
产生M个-1、1等概随机码 产生各占空比双极性归零码波形 双极性归零码的功率谱密度 作图 .. ..
. .
五、仿真源代码
1、单极性
（1）占空比25%
..
..
. .
.. ..
. .
（2）占空比50%
.. ..
. .
.. ..
北邮通信原理软件实验报告
