MATLAB产生特定功率谱密度的高斯白噪声的两种方法

最新推荐文章于 2025-02-27 00:43:09 发布

Retrospector

最新推荐文章于 2025-02-27 00:43:09 发布

阅读量2.1w

点赞数 23

分类专栏： MATLAB 文章标签： MATLAB 高斯白噪声 WGN 功率谱密度

本文链接：https://blog.csdn.net/Ocean_waver/article/details/102759840

版权

产生方法

假设我们想要的高斯白噪声信号的功率谱密度为k（这里只考虑k为常数的情形）。

1、wgn函数：

noise = wgn(1,N,k*B,'linear');             % 功率P = 功率谱密度*B = k*B；

其中前两个参数表示行列数，第三个参数表示信号功率P，要注意区分功率和功率谱密度的概念。此外，linear表示线性强度（linear power），单位为Watt。如果输入其他：‘dBw’或缺省则表示用dBw作为功率单位，则输入的第三个参数应为10log10(P)=10log10(k*B)‘dBm’表示以dBm作为功率单位，则输入的第三个参数应为10log10(P*1000)=10log10(k*B*1000)。下面会再次对此进行总结。

P(信号功率)=k(信号功率谱密度)* B(带宽);			% 这里k为常数；

具体原因是：我们所能产生的噪声总是离散且有限长的，不妨记长度为N，有限长的时域信号对应无限宽的频谱，但是我们只能通过DFT获得有限长的频谱，频谱长度由时域采样间隔td决定，即带宽为B = 1/td，td无法做到无穷小，那么B就不可能无限宽，这表明从频域来看，我们只能获得噪声的长度为 B 的频谱，那么功率谱的频率范围也为B，这样就可以得到P(功率)=k(功率谱密度)* B(带宽)。
我们可以通过下面的代码块来验证：

clear all; close all;
N = 2048;           % 噪声点数2048点；
B = 20000;          % 带宽20000Hz；
k = 1;              % 功率谱密度1；
f = 0-B/2:B/N:B-B/N-B/2;    % 频谱列表
noise = wgn(1,N,k*B,'linear');      % 功率P = 功率谱密度*B = k*B；
Nf =

最低0.47元/天解锁文章