MATLIB产生高斯白噪声并生成C语言代码

陈拓 2020/12/10-2020/12/10

 

我要在他处使用C语言产生高斯白噪声，先用MATLIB生成一个能产生高斯白噪声的C程序作为参考。

1. 高斯白噪声

百度百科，高斯白噪声(White Gaussian Noise，WGN)：如果一个噪声，它的瞬时值服从高斯分布（正态分布），而它的功率谱密度又是均匀分布的，则称它为高斯白噪声。

2. 用BATLIB产生高斯白噪声

• BATLIB帮助中心

https://ww2.mathworks.cn/help/matlab/ref/fft.html

• 新建一个函数文件

• 使用MATLIB的wgn函数产生高斯白噪声

查看语法说明：

wgn(m,n,power)产生一个m*n高斯白噪声样本的矩阵，单位为伏特。噪声的功率power以dbW为单位。

• 编写函数

保存

我的保存路径是E:\iMATLIB

• 调用iwgn函数产生高斯白噪声

测试命令：

fs=100;
T=10; %时间10s
n=round(T*fs);%采样点个数
y=iwgn(1,n,0);%信号功率0dbW（等于1W）
subplot(2,1,1),plot(y); %将图分成2行1列，在第一行绘图

• 高斯白噪声频谱

在下面命令中：

fft是快速傅里叶变换；fftshift通过将零频分量移动到数组中心，重新排列傅里叶变换。

linspace生成线性间距向量。MATLIB中的向量是只有一行元素的数组，向量中的单个项通常称为元素。Matlab中的向量索引值从1开始，而不是从0开始。

fft_y=fftshift(fft(y));
f=linspace(-fs/2,fs/2,n); % 生成-fs/2,fs/2之间的n个点
figure;%打开绘图窗口
plot(f,abs(fft_y));%绘图
title('高斯白噪声功率谱');%加标题
xlabel('频率');% x轴标题
ylabel('功率');% y轴标题

3. 生成C语言程序

• 在APP中找到Matlab Coder

• 打开MATLIB Coder

点击3个点，选择已经保存的函数文件iwgn.m，进入函数的入口

点击右下角Next。

• 定义输入参数类型

输入一个调用该函数的测试例子，例如：iwgn(1,1000,0)。

然后点击Autodefine Input Types按钮，matlab会根据测试代码自动确定形参的类型。

点击右下角Next。

• 检查能否转换为C函数

点击check for issues按钮。

没有问题，点击右下角Next。

• 生成代码

点击generate按钮，生成C函数代码。

代码生成成功。

可以看到在matlab中wgn函数调用randn函数产生高斯白噪声。

• 查看报告

点击右下角Next：

 

参考文档

1. MATLAB wgn函数产生高斯白噪声程序https://jingyan.baidu.com/article/e3c78d644dbdec3c4d85f56c.html