matlab 伪随机干扰信号

以下是Matlab生成伪随机干扰信号的示例代码：

% 设置参数
fs = 1000; % 采样率
f0 = 50; % 基频
N = 1000; % 采样点数
M = 10; % 干扰音数量
f_min = 100; % 干扰音最小频率
f_max = 200; % 干扰音最大频率
P_min = 0.1; % 干扰音最小功率
P_max = 0.5; % 干扰音最大功率

% 生成基频信号
t = (0:N-1)/fs;
u = sin(2*pi*f0*t);

% 生成干扰信号
f = linspace(f_min, f_max, M);
P = linspace(P_min, P_max, M);
phi = 2*pi*rand(1, M);
v = zeros(1, N);
for i = 1:M
    v = v + sqrt(P(i))*sin(2*pi*f(i)*t + phi(i));
end

% 生成干扰信号和基频信号的和
x = u + v;

% 绘制信号波形和功率谱
figure;
subplot(2,1,1);
plot(t, x);
title('伪随机干扰信号波形');
xlabel('时间(s)');
ylabel('幅度(V)');
xlim([0 N/fs]);
subplot(2,1,2);
periodogram(x, [], length(x), fs, 'onesided');
title('伪随机干扰信号功率谱');

该代码生成了一个基频为50Hz，采样率为1000Hz，采样点数为1000的正弦信号，并在其中加入了10个随机干扰音。干扰音的频率在100Hz到200Hz之间均匀分布，功率在0.1到0.5之间均匀分布，相位是在0到2π之间均匀分布的随机值。最后，将基频信号和干扰信号相加得到伪随机干扰信号，并绘制了信号波形和功率谱。