使用matlab：数据中包含3个不同频率正弦波，1个特定频率的方波，以及高斯白噪声，采样率为30000HZ， 使用“filterDesigner”toolbox设计FIR滤波器（4个滤波器），分别提取其中3个不同频率正弦 波，1个频率的方波。正弦波的频率分别为110Hz，2009Hz，5015Hz，分别画出滤波之后的3个不同频率正弦波，1个频率的方波的时域波形和fft频谱

好的，我会根据您的要求进行回答。以下是解决方案：

1. 首先使用Matlab中的“filterDesigner”工具箱设计4个FIR滤波器，分别用于提取110Hz、2009Hz、5015Hz的正弦波和1个特定频率的方波。在“filterDesigner”中，选择“FIR”滤波器类型，选择“Bandpass”滤波器响应类型，并设置相应的频率范围和通带、阻带参数。最后，导出每个滤波器的系数。

2. 在Matlab中，读入包含3个不同频率正弦波、1个特定频率方波和高斯白噪声的数据，并将数据进行归一化处理。

3. 使用Matlab中的“filter”函数，将每个滤波器的系数应用于原始数据，分别提取出3个不同频率正弦波和1个频率的方波。具体实现代码如下：

% 读入数据
data = load('data.mat');
x = data.x;

% 归一化处理
x = x / max(abs(x));

% 导入滤波器系数
h110 = data.h110;
h2009 = data.h2009;
h5015 = data.h5015;
hsq = data.hsq;

% 应用滤波器
y110 = filter(h110, 1, x);
y2009 = filter(h2009, 1, x);
y5015 = filter(h5015, 1, x);
ysq = filter(hsq, 1, x);

% 提取正弦波和方波
fs = 30000;
t = (0:length(x)-1) / fs;
sine110 = sin(2*pi*110*t);
sine2009 = sin(2*pi*2009*t);
sine5015 = sin(2*pi*5015*t);
square = square(2*pi*1000*t);

sine110_filtered = y110 .* sine110;
sine2009_filtered = y2009 .* sine2009;
sine5015_filtered = y5015 .* sine5015;
square_filtered = ysq .* square;

4. 绘制每个信号的时域波形和FFT频谱，代码如下：

% 绘制110Hz正弦波的时域波形和FFT频谱
figure;
subplot(2,1,1);
plot(t, sine110_filtered);
title('110Hz Sine Wave (Filtered)');
xlabel('Time (s)');
ylabel('Amplitude');
subplot(2,1,2);
N = length(sine110_filtered);
f = (-N/2:N/2-1) * fs / N;
Y = fftshift(fft(sine110_filtered)/N);
plot(f, abs(Y));
title('110Hz Sine Wave FFT (Filtered)');
xlabel('Frequency (Hz)');
ylabel('Magnitude');

% 绘制2009Hz正弦波的时域波形和FFT频谱
figure;
subplot(2,1,1);
plot(t, sine2009_filtered);
title('2009Hz Sine Wave (Filtered)');
xlabel('Time (s)');
ylabel('Amplitude');
subplot(2,1,2);
N = length(sine2009_filtered);
f = (-N/2:N/2-1) * fs / N;
Y = fftshift(fft(sine2009_filtered)/N);
plot(f, abs(Y));
title('2009Hz Sine Wave FFT (Filtered)');
xlabel('Frequency (Hz)');
ylabel('Magnitude');

% 绘制5015Hz正弦波的时域波形和FFT频谱
figure;
subplot(2,1,1);
plot(t, sine5015_filtered);
title('5015Hz Sine Wave (Filtered)');
xlabel('Time (s)');
ylabel('Amplitude');
subplot(2,1,2);
N = length(sine5015_filtered);
f = (-N/2:N/2-1) * fs / N;
Y = fftshift(fft(sine5015_filtered)/N);
plot(f, abs(Y));
title('5015Hz Sine Wave FFT (Filtered)');
xlabel('Frequency (Hz)');
ylabel('Magnitude');

% 绘制方波的时域波形和FFT频谱
figure;
subplot(2,1,1);
plot(t, square_filtered);
title('Square Wave (Filtered)');
xlabel('Time (s)');
ylabel('Amplitude');
subplot(2,1,2);
N = length(square_filtered);
f = (-N/2:N/2-1) * fs / N;
Y = fftshift(fft(square_filtered)/N);
plot(f, abs(Y));
title('Square Wave FFT (Filtered)');
xlabel('Frequency (Hz)');
ylabel('Magnitude');

以上就是实现您要求的所有功能的Matlab代码和结果展示。