用Kaiser窗设计一个FIR数字带阻滤波器，对模拟信号xa(t)滤波， xa(t) = cos (a t) + cos (b t) + cos (c t) 其中，a = 2*pi *6500， b = 2*pi *7000, c = 2*pi *9000 要求滤去7000Hz 的频率成分。系统采样率为 fs = 32000 Hz。这里采样点数比较大，可以用 N = 4096。 滤波器的 Rp = 0.25 dB， As = 50 dB， 过渡带宽可以用模拟频率（例如200Hz）也可以用数字频率指定。还可以改变As（比如30dB）观察滤波效果。

根据题目要求，可以按照以下步骤设计FIR数字带阻滤波器：

1. 确定滤波器的通带和阻带频率以及对应的通带和阻带衰减量，以及采样频率。本题中，通带范围为[6000, 8000]Hz，阻带范围为[0, 5000]Hz，通带最大衰减量Rp为0.25dB，阻带最小衰减量Rs为50dB，采样频率为32000Hz。

2. 根据通带和阻带频率，计算出通带和阻带中心频率，然后通过带宽公式计算出截止频率。通带中心频率为7000Hz，阻带中心频率为6000Hz和8000Hz，过渡带宽可以根据数字频率指定，例如200Hz。

3. 计算出滤波器的阶数，可以使用kaiserord函数实现，需要指定通带衰减量、阻带衰减量、过渡带宽和窗函数类型等参数。

4. 使用fir1函数设计出FIR滤波器系数。

5. 将模拟信号进行数字化，然后使用filter函数进行滤波。

下面是MATLAB中的代码实现：

% 设计参数
a = 2*pi*6500;
b = 2*pi*7000;
c = 2*pi*9000;
Wp = [a b];        % 通带频率范围
Ws = [a-100 b+100]; % 阻带频率范围
Rp = 0.25;         % 通带最大衰减量
Rs = 50;           % 阻带最小衰减量
dev = [(10^(Rp/20)-1)/(10^(Rp/20)+1)  10^(-Rs/20)]; % 计算通带和阻带衰减量
f = [0 Ws(1) Wp(1) Wp(2) Ws(2) 32000/2]/(32000/2);   % 计算频率向量
a = [0 0 1 1 0 0]; % 定义通带和阻带的增益
n = kaiserord(f,a,dev);  % 计算出滤波器的阶数
b = fir1(n, Ws/(32000/2), 'bandstop', kaiser(n+1, 5)); % 设计kaiser窗口FIR数字带阻滤波器

% 模拟信号xa(t)
t = 0:1/32000:4095/32000;
xa = cos(a*t) + cos(b*t) + cos(c*t);

% 对模拟信号进行数字化
x = xa/max(abs(xa));

% 应用数字带阻滤波器
y = filter(b, 1, x);

% 绘制信号和滤波后的信号
figure;
subplot(2,1,1);
plot(t, x);
title('原始信号');
xlabel('时间');
ylabel('幅值');
subplot(2,1,2);
plot(t, y);
title('滤波后信号');
xlabel('时间');
ylabel('幅值');

运行上述代码，可以得到滤波后的信号图像。可以尝试改变阻带最小衰减量Rs的值，观察滤波效果的变化。