如果是均值为1的白噪声通过一个带通系统，请再写一遍程序

好的，下面的程序演示了如何生成均值为1的白噪声，并通过一个带通系统进行滤波：

Fs = 44100; % 采样率
T = 5; % 信号持续时间
t = 0:1/Fs:T;
P = 1; % 平均功率
x = sqrt(P)*randn(size(t)) + 1; % 生成均值为1的白噪声信号

% 带通滤波器设计
f1 = 1000; % 通带下限频率
f2 = 5000; % 通带上限频率
Wn = [f1, f2]*2/Fs;
[b,a] = butter(4, Wn, 'bandpass');

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

% 绘制时域和频域图像
subplot(2,1,1), plot(t, y);
title('时域图像');
subplot(2,1,2), pwelch(y, [], [], [], Fs);
title('频域图像');

这个程序首先使用randn函数生成一个长度为T秒的随机信号，然后加上一个常量1，使其成为均值为1的白噪声信号。接着，使用butter函数设计一个4阶带通滤波器，通带下限频率为1000Hz，通带上限频率为5000Hz。最后，使用filter函数对白噪声信号进行滤波处理，得到的输出信号y是经过带通滤波器滤波后的有色噪声信号。最后绘制出该信号的时域和频域图像。