matlab代码导入原始信号音频和噪声两个音频文件，将之混合叠加后，使用rls方法从混合音频中提取出原始音频信号，结果画出原始音频，混合音频和提取出的原始音频的时域图和频域图共6个图

下面是一份matlab代码示例，用于演示如何使用rls方法提取原始音频信号。请注意，这只是一个示例代码，您需要适当修改以适应您的具体情况。

%% 导入原始音频和噪声
[y, Fs_y] = audioread('original_audio.wav');
[n, Fs_n] = audioread('noise_audio.wav');

%% 混合两个音频
mix = y + n;

%% RLS提取原始音频信号
% 设置参数
L = 100; % RLS滤波器长度
delta = 0.01; % RLS滤波器的遗忘因子
lambda = 0.99; % RLS滤波器的收敛速度
P0 = 1000 * eye(L); % 初始协方差矩阵
w = zeros(L, 1); % 初始权重向量

% RLS滤波器
for i = L:length(mix)
    x = flip(mix(i-L+1:i)); % 输入信号
    d = y(i); % 目标信号
    P = lambda * P0 + x * x'; % 更新协方差矩阵
    K = P * x / (delta + x' * P * x); % 计算增益
    e = d - w' * x; % 计算误差
    w = w + K * e; % 更新权重向量
    P0 = (eye(L) - K * x') * P; % 更新协方差矩阵
end

%% 画出时域图和频域图
t = (0:length(y)-1) / Fs_y;
figure;
subplot(3, 2, 1);
plot(t, y);
title('原始音频');
subplot(3, 2, 3);
plot(t, mix);
title('混合音频');
subplot(3, 2, 5);
plot(t, w' * fliplr(mix(L:end)));
title('提取出的原始音频');
xlabel('时间/s');

f = (0:length(y)-1) / length(y) * Fs_y;
Y = fft(y);
MIX = fft(mix);
W = fft(w' * fliplr(mix(L:end)));
subplot(3, 2, 2);
plot(f, abs(Y));
title('原始音频');
subplot(3, 2, 4);
plot(f, abs(MIX));
title('混合音频');
subplot(3, 2, 6);
plot(f, abs(W));
title('提取出的原始音频');
xlabel('频率/Hz');

注意：在运行此代码之前，请将原始音频和噪声音频文件（格式为.wav）放在matlab当前工作目录下，并将代码中的文件名替换为您的文件名。