编写一段程序在matlab环境中计算一段语音波形的 90%带宽，98%带宽

在MATLAB环境中，计算语音波形的90%带宽和98%带宽通常涉及到频谱分析，可以使用快速傅立叶变换（FFT）对信号进行离散化处理，然后找出包含特定能量百分比的频率范围。以下是基本步骤：

1. **加载音频文件**：

   [y, Fs] = audioread('your_audio_file.wav'); % 替换为实际的音频文件名

2. **计算FFT**：

   Y = fft(y); % 对信号取离散傅立叶变换
   P = abs(Y).^2; % 获取功率谱
   freq = (0:length(P)-1) * Fs / length(P); % 生成频率向量

3. **计算90%带宽**：
使用累积分布函数（CDF）找到包含90%能量的两个频率点，即`freq(upperIndex90)`和`freq(lowerIndex90)`：

   upperIndex90 = find(cumsum(P) >= 0.9, 1);
   lowerIndex90 = max(1, upperIndex90 - 1); % 如果只有一个元素，防止下标越界
   bw90 = freq(upperIndex90) - freq(lowerIndex90);

4. **计算98%带宽**：
同样地，找到包含98%能量的频率点`freq(upperIndex98)`和`freq(lowerIndex98)`：

   upperIndex98 = find(cumsum(P) >= 0.98, 1);
   lowerIndex98 = max(1, upperIndex98 - 1);
   bw98 = freq(upperIndex98) - freq(lowerIndex98);

5. **显示结果**：

   fprintf('90%带宽: %.2f Hz\n', bw90);
   fprintf('98%带宽: %.2f Hz\n', bw98);

注意：这个过程假设信号是对称的，如果不对称，可能需要一些额外处理。此外，对于非均匀采样的信号，上述公式可能需要调整。