C中的声波锯齿

时间:2018-11-27 19:32:15

标签: c audio wav

我正在尝试在c中生成锯齿波。 我已经快要结束工作了,但是我有一些必须解决的问题。 我在下面附上了衣服。

#include <stdbool.h>
#include <stdio.h>
#include <stdlib.h>
#include <tgmath.h>

/******************************
*  Magic file format strings. *
******************************/
const char fChunkID[]     = {'R', 'I', 'F', 'F'};
const char fFormat[]      = {'W', 'A', 'V', 'E'};
const char fSubchunk1ID[] = {'f', 'm', 't', ' '};
const char fSubchunk2ID[] = {'d', 'a', 't', 'a'};

/********************************
* WriteWavePCM() configuration: *
* - 2 channels,                 *
* - frequency 44100 Hz.         *
********************************/
const unsigned short N_CHANNELS = 2;
const unsigned int SAMPLE_RATE = 48000;
const unsigned short BITS_PER_BYTE = 8;

bool WriteWavePCM(short* sound, size_t pairAmount, char* fileName){
    const static unsigned int fSubchunk1Size = 16;
    const static unsigned short fAudioFormat = 1;
    const static unsigned short fBitsPerSample = 16;

    unsigned int fByteRate = SAMPLE_RATE * N_CHANNELS * fBitsPerSample / BITS_PER_BYTE;

    unsigned short fBlockAlign = N_CHANNELS * fBitsPerSample / BITS_PER_BYTE;
    unsigned int fSubchunk2Size;
    unsigned int fChunkSize;

    FILE* fout;
    size_t ws;

    if (!sound || !fileName || !(fout = fopen( fileName, "w" ))) return false;

    fSubchunk2Size = pairAmount * N_CHANNELS * fBitsPerSample / BITS_PER_BYTE;
    fChunkSize = 36 + fSubchunk2Size;

    // Writing the RIFF header:
    fwrite(&fChunkID, 1, sizeof(fChunkID),      fout);
    fwrite(&fChunkSize,  sizeof(fChunkSize), 1, fout);
    fwrite(&fFormat, 1,  sizeof(fFormat),       fout);

    // "fmt" chunk:
    fwrite(&fSubchunk1ID, 1, sizeof(fSubchunk1ID),      fout);
    fwrite(&fSubchunk1Size,  sizeof(fSubchunk1Size), 1, fout);
    fwrite(&fAudioFormat,    sizeof(fAudioFormat),   1, fout);
    fwrite(&N_CHANNELS,      sizeof(N_CHANNELS),     1, fout);
    fwrite(&SAMPLE_RATE,     sizeof(SAMPLE_RATE),    1, fout);
    fwrite(&fByteRate,       sizeof(fByteRate),      1, fout);
    fwrite(&fBlockAlign,     sizeof(fBlockAlign),    1, fout);
    fwrite(&fBitsPerSample,  sizeof(fBitsPerSample), 1, fout);

    /* "data" chunk: */
    fwrite(&fSubchunk2ID, 1, sizeof(fSubchunk2ID),      fout);
    fwrite(&fSubchunk2Size,  sizeof(fSubchunk2Size), 1, fout);

    /* sound data: */
    ws = fwrite(sound, sizeof(short), pairAmount * N_CHANNELS, fout);
    fclose(fout);
    return true;
}

////////////////////////////////////////////////

const unsigned int N_SAMPLE_PAIRS = 50000;

int main(int argc, char* argv[]){
    short* sound;
    int i;
    int j;
    bool status;
    char* file_name;
    int l; // dodane

    sound = (int*) malloc(sizeof(int) * N_SAMPLE_PAIRS * N_CHANNELS );

    if (!sound)
    {
        puts("Could not allocate space for the sound data.");
        return (EXIT_FAILURE);
    }

    int amplitude = 10000;
    int frequency = 80;
    short record = 0;
    short waveNumber = 1;
    int samplesPerWavelength = SAMPLE_RATE / (frequency/N_CHANNELS);
    int soundLen = 10 * samplesPerWavelength;
    int ampStep = (int)((amplitude*2)/(int)samplesPerWavelength);
    short step = 5*samplesPerWavelength;
    short muteRate = amplitude/(soundLen/samplesPerWavelength);
    int totalSamplesWritten = 0;
    int tempSample =0;

    for (i=0, j=0; i<N_SAMPLE_PAIRS*N_CHANNELS; i+=2, j++) {
        ampStep = (int)((amplitude*2)/(int)samplesPerWavelength);
        tempSample = (int)((totalSamplesWritten%samplesPerWavelength)*ampStep);
        sound[i] = tempSample;
        sound[i+1] = tempSample;
        totalSamplesWritten++;
    }


    file_name = argc > 1 ? argv[1] : "Default2.wav";
    status = WriteWavePCM(sound, N_SAMPLE_PAIRS, file_name);


    free(sound);

    if (status)
    {
        printf("Discotheque is ready in \"%s\"\n", file_name);
    }
    else
    {
        puts( "Something seems to have gone wrong." );
        return (EXIT_FAILURE);
    }

    return 0;
}

这是我的结果。有锯齿,但行为奇怪。我不知道为什么会这样。我检查了样品样本,一切看起来都很好。

enter image description here

有人可以帮助我吗?我被卡住了,我不知道该怎么办。我在2个不同的程序中进行了检查:audiocity和wavepadaudio,然后在这两个程序中都存在这种行为。所以我肯定对我的代码有疑问。请帮助我。

enter image description here

2 个答案:

答案 0 :(得分:9)

当我查看文件转储(MSVC编译)时,我注意到一些FileNotFoundError数据对。然后我改变了这个

python3 -c "print(isinstance(FileNotFoundError(), IOError))"

为此指定一个二进制文件

True

现在小故障消失了。

请注意

  

如果try: content = open("somefile.txt").read() except IOError: # Works in both Python 2 & 3 print("Oops, we can not read this file") 中未提供0D 0Afout = fopen( fileName, "w" ) ,则默认翻译模式为   由全局变量fout = fopen( fileName, "wb" ) 定义。

因此,默认的翻译模式必须为“文本”。

答案 1 :(得分:1)

这是我得到的输出的一些屏幕截图。对我来说似乎很好。

使用内置规格。 COLLECT_GCC = gcc COLLECT_LTO_WRAPPER = / usr / lib / gcc / i686-pc-cygwin / 5.4.0 / lto-wrapper.exe 目标:i686-pc-cygwin 配置为:/cygdrive/i/szsz/tmpp/gcc/gcc-5.4.0-1.i686/src/gcc-5.4.0/configure --srcdir = / cygdrive / i / szsz / tmpp / gcc / gcc- 5.4.0-1.i686 / src / gcc-5.4.0 --prefix = / usr --exec-prefix = / usr --localstatedir = / var --sysconfdir = / etc --docdir = / usr / share / doc / gcc --htmldir = / usr /共享/文档/ gcc / h tml -C --build = i686-pc-cygwin --host = i686-pc-cygwin --target = i686-pc-cygwin --without-libiconv-prefix --without-libintl-prefix --libexecdir = / usr / lib --enable-shared --enable-shared-libgcc --enable-static --enable-version-specific-runtime-libs --enable-bootstrap --enable -__ cxa_ate xit --with-dwarf2 --with-arch = i686 --with-tune = generic --disable-sjlj-exceptions --enable-languages = ada,c,c ++,fortran,java,lto,objc,obj-c ++ --enable-graphite --enable-threads = posix --enable-libatomic --enable-libcilkrts --enable-libgomp --enable-libitm --enable-libquadmath --enab le-libquadmath-support --enable-libssp --enable-libada --enable-libjava --enable-libgcj-sublibs --disable-java-awt --disable-symvers --with- ecj-jar = / usr / share / java / ecj.jar --with-gnu-ld --with-gnu-as --with-cloog-include = / usr / include / cloog-isl --without-libiconv-前缀-不- libintl前缀--with-system-zlib --enable-linker-build-id --with-default-libstdcxx-abi = gcc4兼容 螺纹型号:posix gcc版本5.4.0(GCC)

uname -a CYGWIN_NT-6.1-WOW DE-E03895 2.8.0(0.309 / 5/3)2017-04-01 20:42 i686 Cygwin

Here it is

Zoomed in

这是输出中损坏的样本的屏幕截图;

enter image description here