MATLAB设计IIR滤波器设计小星星钢琴曲保留1音，滤除其他音

首先，需要明确小星星钢琴曲的音调和频率。小星星的音符为C4、C4、G4、G4、A4、A4、G4，频率分别为261.63 Hz、261.63 Hz、392 Hz、392 Hz、440 Hz、440 Hz、392 Hz。

在MATLAB中，可以使用fdatool工具箱来设计IIR滤波器。具体步骤如下：

1. 打开MATLAB，进入fdatool工具箱。
2. 在“Filter Type”中选择“IIR”滤波器。
3. 在“Design Method”中选择“Elliptic”设计方法，因为它可以实现比其他方法更陡峭的滤波器，同时可以达到更高的阻带衰减。
4. 在“Filter Order”中选择适当的阶数，一般建议选择8-10阶。阶数越高，滤波器的性能越好，但计算量也会增加。
5. 在“Passband Specification”中设置通带频率范围，因为只需要保留C4音，所以选择260-263 Hz作为通带范围。
6. 在“Stopband Specification”中设置阻带频率范围，因为需要滤除其他音，所以选择200-600 Hz作为阻带范围。
7. 点击“Design Filter”按钮，生成滤波器的系数。
8. 点击“View Filter Response”按钮，查看滤波器的频率响应，确保滤波器的性能符合要求。
9. 将生成的滤波器系数导入MATLAB代码中，对小星星钢琴曲进行滤波。

下面是MATLAB代码示例：

% 小星星钢琴曲频率
freq = [261.63 261.63 392 392 440 440 392];

% 设计IIR滤波器
n = 8; % 阶数
fs = 8000; % 采样率
wp = [260 263]/(fs/2); % 通带频率范围
ws = [200 600]/(fs/2); % 阻带频率范围
Rp = 1; % 通带最大衰减
Rs = 60; % 阻带最小衰减
[b, a] = ellip(n, Rp, Rs, wp, 'bandpass'); % 生成滤波器系数

% 读取小星星钢琴曲音频文件
[x, fs] = audioread('little_star.wav');

% 对小星星钢琴曲进行滤波
y = filter(b, a, x);

% 保存保留C4音的小星星钢琴曲
audiowrite('little_star_filtered.wav', y, fs);

运行以上代码后，会生成一个保留C4音的小星星钢琴曲文件“little_star_filtered.wav”。