matlab不同频率声音的声强如何归一

% 读取音频文件
[s1,fs1] = audioread('sound1.wav');
[s2,fs2] = audioread('sound2.wav');

% 对原始信号进行短时傅里叶变换，得到语音信号的功率谱和噪声功率谱
window_length = 1024; % 窗口长度
hop_size = 512; % 帧移
[~,F,T,Pxx_1] = spectrogram(s1,hamming(window_length),window_length-hop_size,window_length,fs1);
[~,~,~,Pxx_2] = spectrogram(s2,hamming(window_length),window_length-hop_size,window_length,fs2);

% 计算带噪语音的信噪比
SNR = zeros(size(Pxx_1));
for i = 1:size(Pxx_1,2)
    noise_Pxx = median(Pxx_1(:,max(1,i-10):min(size(Pxx_1,2),i+10)),2); % 计算噪声功率谱，使用邻近帧的平均值
    signal_Pxx = Pxx_1(:,i) - noise_Pxx; % 计算信号功率谱
    SNR(:,i) = signal_Pxx./noise_Pxx; % 计算信噪比
end

% 利用信噪比计算噪声功率谱的加权平均值
SNR_mean = mean(SNR(SNR>0)); % 去除负值
noise_Pxx_1 = median(Pxx_1(:,SNR<SNR_mean*0.9),2); % 信噪比小于平均值90%的帧作为噪声帧
noise_Pxx_2 = median(Pxx_2(:,SNR<SNR_mean*0.9),2);

% 计算纯净语音功率谱
Pxx_clean_1 = max(0, Pxx_1 - noise_Pxx_1); % 将小于0的幅度值置为0 
Pxx_clean_2 = max(0, Pxx_2 - noise_Pxx_2); 

% 合成增强后的语音信号
s1_enhanced = istft(sqrt(Pxx_clean_1).*exp(1j*angle(stft(s1,window_length,hop_size,window_length,fs1))),window_length,hop_size,window_length,fs1);
s2_enhanced = istft(sqrt(Pxx_clean_2).*exp(1j*angle(stft(s2,window_length,hop_size,window_length,fs2))),window_length,hop_size,window_length,fs2);

% 可视化增强前后的语音信号
time = (0:length(s1)-1)/fs1;
figure;
subplot(2,2,1), plot(time,s1), title('Original sound 1');
subplot(2,2,2), plot(time,s2), title('Original sound 2');
subplot(2,2,3), plot(time,s1_enhanced), title('Enhanced sound 1');
subplot(2,2,4), plot(time,s2_enhanced), title('Enhanced sound 2');