MATLAB生成带宽小于40MHz的平稳随机信号 - 代码示例 - 常规

下面是一个示例代码，生成一个带宽为20MHz的平稳随机信号，并使用一个带通滤波器进行滤波，最后将信号与噪声混合以得到一个信噪比为15的信号：

% 生成平稳随机信号
Fs = 200e6;  % 采样率
T = 1/Fs;  % 采样间隔
N = 1e6;  % 信号长度
f0 = 10e6;  % 信号中心频率
B = 20e6;  % 信号带宽
t = (0:N-1)*T;  % 时间序列
x = randn(1,N).*exp(2*pi*1i*f0*t);  % 平稳随机信号

% 设计带通滤波器
f1 = f0 - B/2;  % 通带下限频率
f2 = f0 + B/2;  % 通带上限频率
Wp = [f1 f2]/(Fs/2);  % 归一化通带频率
Rp = 1;  % 通带最大衰减
Rs = 60;  % 阻带最小衰减
[n, Ws] = cheb2ord(Wp, [f1-10e6 f2+10e6]/(Fs/2), Rp, Rs);  % 计算滤波器阶数和归一化阻带频率
[b,a] = cheby2(n, Rs, Ws);  % 设计Chebyshev Type II滤波器

% 滤波
y = filter(b,a,x);

% 添加噪声
SNR = 15;  % 信噪比
Psignal = mean(abs(y).^2);  % 信号功率
Pnoise = Psignal/10^(SNR/10);  % 噪声功率
noise = sqrt(Pnoise)*randn(1,N);  % 生成噪声
z = y + noise;  % 混合信号

代码解析

生成平稳随机信号: 该代码段使用 randn 函数生成一个长度为 N 的随机信号，并将其乘以一个复指数项以使其具有中心频率 f0。
设计带通滤波器: 使用 cheb2ord 函数计算 Chebyshev Type II 滤波器的阶数 n 和归一化阻带频率 Ws。然后使用 cheby2 函数设计滤波器系数 b 和 a。
滤波: 使用 filter 函数对生成的随机信号进行滤波。
添加噪声: 该代码段通过控制信噪比 SNR 来计算噪声功率 Pnoise，然后生成一个均值为零、方差为 Pnoise 的高斯噪声，最后将噪声叠加到滤波后的信号上。

修改后的代码:

该代码可以轻松修改以满足带宽小于 40MHz 的要求，只需更改 B 的值即可。此外，可以调整其他参数以满足特定的需求，例如：

采样率 Fs: 更高的采样率可以提高信号的精度，但也会导致更高的计算量。
信号长度 N: 较长的信号长度可以提供更多的信息，但也会增加计算时间。
信噪比 SNR: 较高的信噪比可以获得更清晰的信号，但也会降低信号的真实性。

总结

该代码演示了如何使用 MATLAB 生成带宽小于 40MHz 的平稳随机信号，并对其进行滤波和噪声添加。用户可以根据自己的需要调整代码参数以生成符合要求的信号。