Question

这是我的问题。我正在连续时间和离散时域中实现高通和低通滤波器，而不使用MATLAB内置函数或信号处理工具箱。我需要过滤上传的wav文件。我有以下代码实现低通连续时间过滤器。主要问题是，当我去播放文件时，我收到以下错误：

警告：在写入文件期间剪切数据：flute_new   在wavwrite＆gt; PCM_Quantize at 286

在wavwrite＆gt; write_wavedat at 302

在139的wavwrite

在23的test2中

我一直在读这个问题的常见解决方案是：signal = signal / max（abs（signal）），但仍然无法解决剪切问题。我不太明白是什么产生了这个错误，所以任何帮助都将不胜感激。谢谢你的阅读。代码如下：

    [y,Fs,nbits]=wavread('flute.wav');
    m = length(y); %96375
    Fc = 500;               %% frequency cut
    nbits
    n = -(m-1)/2 : (m-1)/2;            %% kernel length / samples
    %Fs = Fs;    %22050      %% sample frequency / (nyquist - %2*m)
    Tc = 2*pi*(Fc/Fs);      %% theta c
    Hm = sin(Tc*n) ./ (pi*n);
        if (rem(m, 2) == 1)
          index = (m+1)/2;
          Hm(index) = Tc/pi;
        end

    [H w] = freqz(Hm, 1, 1024); % 1024 - Amount of Fourier Transform points


    %do playback
    sound(w,Fs);
    %save and playback the file
    wavwrite(w,'flute_new')

更新：根据我的理解，因为我的wav文件是16位，我只能有-1＆lt; = w＆lt; 1。我按照评论中的建议更新了我的代码，但是当我检查w值时，我最后仍然有一个w = 1。建议的代码确实消除了剪辑警告，但通过扬声器播放音频文件仍然有剪辑。也许我理解错误的概念。对不起，如果是这样的话。

播放代码（已更新）：

%Write and play file
w = w./max(abs(w(:)))*(1-(2^-(nbits-1)));
w
%do playback
    sound(w,Fs);
    wavwrite(w,'flute_new');

终端输出：

Answer 1

这不是错误，只是警告。

如果您编写具有32位分辨率的.wav文件，则允许的数据范围为-1 <= w <= 1.但是对于8,16或24位，限制为-1 <= w ＆LT; 1，使得w == 1产生剪切警告（1位短'1'，即16位模式下1 - 1/32768 = 0.99996948242188。），并且它不能写入文件。你可以忽略它，或者：

format long
w = w./max(abs(w(:)))*(1-(2^-(nbits-1)));
wavwrite(w,'flute_new');

在MATLAB中播放音频文件时的Wav文件剪辑

1 个答案: