（1）首先计算已知信号序列（采样之后得到的信号）的平均功率。该序列在第n个点处的功率为：

如果已知的信号序列中的总共的点数为N个，则该序列的平均功率为：

在MATLAB中求平均功率的方法是：

Pav=sum(x.^2)/length(x);

（2）第二步是求单个符号的能量。能量的定义是功率乘以时间。对于单个符号来说，因为已经被采样了，每个符号可能对应多个采样点。因此，此处需要已知符号速率。

每个符号的能量为：

在MATLAB中求单位符号能量的方法是：

Eb=sum(x.^2)/(length(x)*fb);

（3）第三步是求噪声的单边功率谱密度。根据Eb/N0即可，注意首先要把dB转化为实际的倍数。若已知的Eb/N0是dB单位的，则转化过程如下：

ebn0=10^((EbN0)/10);

随后得到n0：

n0=Eb/ebn0;

（4）第四步是根据噪声的单边功率谱密度求噪声的方差。对于实信号来说，如果信号的采样频率是fs，则噪声带宽为采样频率的一半。噪声的方差即为噪声的平均功率，噪声的平均功率是单边功率谱密度和噪声带宽的乘积，因此噪声的方差为：

MATLAB中的计算过程：

pn=n0*fs/2;（从实际仿真的结果来看，乘以fs是符合理论值的，该内容需要进一步核实）

（5）第五步即生成噪声序列。

NOISE=sqrt(pn)*randn(1,length(x));

注：参考文献:

Bit-Error-Rate Simulation Using Matlab

作者：JE Gilley

http://xueshu.baidu.com/s?wd=paperuri:(18dc898216ce76b4c166f194fab1b349)&filter=sc_long_sign&sc_ks_para=q%3DBit-Error-Rate+Simulation+Using+Matlab&tn=SE_baiduxueshu_c1gjeupa&ie=utf-8&sc_us=12877983928575481934