请生成一个matlab代码

以下是一个使用布莱克曼窗函数法设计FIR数字带通滤波器的MATLAB代码示例：

% 设计FIR数字带通滤波器
% 采用布莱克曼窗函数法，通带截止频率为0.3*pi，阻带截止频率为0.4*pi
% 通带最大衰减量为1dB，阻带最小衰减量为60dB

% 参数设置
fp = 0.3*pi; % 通带截止频率
fs = 0.4*pi; % 阻带截止频率
Ap = 1; % 通带最大衰减量
As = 60; % 阻带最小衰减量
M = 100; % 滤波器阶数

% 计算带宽和中心频率
BW = fs - fp;
fc = (fs + fp) / 2;

% 计算理想的频率响应
h_ideal = fir1(M, [fp, fs]/(pi/2), 'bandpass');

% 计算布莱克曼窗函数系数
w = blackman(M+1);
b = h_ideal .* w';

% 计算实际频率响应
[H, w] = freqz(b, 1, 1024);

% 归一化频率响应
Hn = abs(H) / max(abs(H));

% 绘制频率响应曲线
figure;
plot(w/pi, Hn);
xlabel('归一化频率');
ylabel('幅度');
title('FIR数字带通滤波器频率响应');

% 计算滤波器的差分方程
b = b / sum(b);
a = 1;

这段代码可以生成一个FIR数字带通滤波器，并绘制出其频率响应曲线。其中，使用了MATLAB内置的fir1函数来计算理想的频率响应，使用了blackman函数来计算布莱克曼窗函数系数，使用了freqz函数来计算频率响应。最后，通过归一化频率响应，得到了幅度响应曲线，并计算出滤波器的差分方程。