Question

我正在尝试使用'fwrite'并制作snd文件。我想制作IIR滤波器。我制作了一个FIR滤波器，我将代码用于IIR滤波器。（当然，改变系数）但我认为'fwrite'不起作用。因为IIR滤波器的结果仅为零。我想我在代码中犯了一些错误。

你能给我一个暗示吗？我现在完全处于恐慌状态。

当我检查输出[]时，它看起来很好。它有整数组，我尊重。

我不知道如何写它。

我知道代码太长了，看起来很难。

但我不知道如果你帮我一把，那就很愉快了。谢谢你的阅读。

#include <stdio.h>
#include <stdint.h>
#include <memory.h>
// Filter Code Definitions
#define MAX_INPUT_LEN 80
#define MAX_FLT_LEN 14
#define BUFFER_LEN (MAX_FLT_LEN - 1 + MAX_INPUT_LEN)

double insamp[BUFFER_LEN];
// IIR inititialization
void IirFloatInit(void)
{
    memset(insamp, 0, sizeof(insamp));
}

// the IIR filter function
void IirFloat(double *coeffs, double *input, double *output, int length, int 
filterLength)
{
    double acc;
    double *coeffp;
    double *inputp;
    int n;
    int k;

// put the new samples at the high end of the buffer
memcpy(&insamp[filterLength - 1], input, length * sizeof(double));
for (n = 0; n < length; n++) {
    coeffp = coeffs;
    inputp = &insamp[filterLength - 1 + n];
    acc = 0;
    for (k = 0; k < filterLength; k++) {
        acc += (*coeffp++) * (*inputp--);
    }
    output[n] = acc;
}
memmove(&insamp[0], &insamp[length], (filterLength - 1) * sizeof(double));
}

double coeffs[MAX_FLT_LEN];

void intToFloat(int16_t *input, double *output, int length)
{
    int i;
    for (i = 0; i < length; i++) {
        output[i] = (double)input[i];
    }
}
void floatToInt(double *input, int16_t *output, int length)
{
    int i;
    for (i = 0; i < length; i++)
    {

        input[i] += 0.5;
        if (input[i] > 32767.0)
        {
            input[i] = 32767.0;
        }
        else if (input[i] < -32768.0)
        {
            input[i] = -32768.0;
        } 
        output[i] = (int16_t)input[i];
    }
}
// number of samples to read per loop
int main(void)
{
    int size;
    int16_t input[MAX_INPUT_LEN];
    int16_t output[MAX_INPUT_LEN];
    double floatInput[MAX_INPUT_LEN];
    double floatOutput[MAX_INPUT_LEN];
    double coeffs[MAX_FLT_LEN];

    FILE *in_fid;
    FILE *out_fid;
    FILE *filter;

    // open the input waveform file
    in_fid = fopen("input.snd", "rb");
    if (in_fid == 0) {
        printf("couldn't open input.snd");
        return;
    }

    // open the output waveform file
    out_fid = fopen("outputFloatIIR.snd", "wb");
    if (out_fid == 0) {
        printf("couldn't open outputFloat.snd");
        return;
    }

    filter = fopen("coeffs_iir.txt", "r");
    if (filter == NULL) {
        puts("couldn't open coeffs_iir.txt");
        return;
    }    
    for (int i = 0; i < MAX_FLT_LEN; i++) {
        fscanf(filter, "%le", &coeffs[i]);
        printf("%le \n", coeffs[i]);
    }


    IirFloatInit();

    do {

        size = fread(input, sizeof(int16_t), MAX_INPUT_LEN, in_fid);

        intToFloat(input, floatInput, size);

        IirFloat(coeffs, floatInput, floatOutput, size, MAX_FLT_LEN);

        floatToInt(floatOutput, output, size);

        fwrite(output, sizeof(int16_t), size, out_fid);
    } while (size != 0);

    fclose(in_fid);
    fclose(out_fid);
    fclose(filter);
    getchar();
    return 0;
}

Answer 1

进入手册（man printf），我可以读到：

eE双参数被舍入并以样式[ - ] d.ddde + -dd转换，其中小数点前有一位数字字符和后面的位数等于精度; 如果缺少精度，它被视为6;如果精度为零，则不显示小数点字符。 E转换使用字母E' (rather than e'）来介绍指数。指数始终包含至少两位数字;如果值为零，则指数为00。

在存储系数的文件中，是否采用该格式？如果使用调试器，读取到fscanf的值会发生什么？

我的意思是，可能格式不是预期的，因此你得到0.也许你想使用fscanf(filter, "%lf", &coeffs[i]);？

fscanf返回的值是多少？进入手册（man fscanf）可以阅读：

返回值这些函数返回分配的输入项数。在a的情况下，这可能比提供的数量少，甚至为零匹配失败。零表示虽然有输入可用，不转换已分配;通常这是由于输入字符无效，例如'％d'的字母字符转换。如果发生输入故障，则返回值EOF 在任何转换之前发生了诸如文件结束之类的情况。如果转换开始后出现错误或文件结尾，则转换次数为何成功完成后返回。

如何在C中使用'fwrite'？

1 个答案: