使用巴特沃斯带通滤波器对采样频率为10Hz的nirs信号，进行0.01-0.1的滤波的MATLAB代码怎么写

在MATLAB中，你可以使用`filter`函数配合`butter`函数来设计巴特沃思带通滤波器，并应用到你的NIRS信号上。假设你的NIRS信号保存在变量`signal`中，采样频率`Fs`为10 Hz，滤波器需要设计为0.01 - 0.1 Hz的带宽。以下是相应的代码：

% 设定滤波器参数
fc = [0.01 0.1]; % 频率范围 (低频、高频)
fs = 10; % 采样频率
order = 4; % 滤波器阶数，可以根据需要调整

% 计算截止频率，注意采样频率单位是Hz
[~, wp] = butter(order, fc/(fs/2), 'bandpass'); 

% 应用滤波器
filtered_signal = filter(wp, [], signal);

% 确保信号长度一致
final_length = length(signal);
if final_length > length(filtered_signal)
    filtered_signal = filtered_signal(1:final_length);
elseif final_length < length(filtered_signal)
    filtered_signal = filtered_signal(1:length(signal));
end

% 显示滤波后的信号
plot(signal, 'r', filtered_signal, 'b');
xlabel('样本索引');
ylabel('信号强度');
title(['0.01-0.1 Hz 巴特沃斯带通滤波后的NIRS信号']);

% 相关问题 --
% 1. 这段代码如何保证滤波器的稳定性？
% 2. 如果信号噪声大，我该如何提高滤波效果？
% 3. 过滤后的信号可能存在失真，怎样减小失真？

这段代码首先生成了一个巴特沃思带通滤波器，然后将该滤波器应用到NIRS信号上。请注意实际使用时可能需要根据你的具体需求调整滤波器的阶数（order）、截止频率以及后续处理步骤。