c - 这个程序有什么错误？

Question

这是我用于查找两个信号的卷积的代码，但是每次有人可以解释我的代码中的错误时，我的输出都变为零了吗？我尝试编译它，但我无法找出我的整个代码是正确的问题是什么提前谢谢

#include<stdio.h>
#include<stdlib.h>
#define array_len(x) (sizeof(x)/sizeof(double))
void convolution(double *signal, int nt, double *wind, int r, double *rm)
{
    int i,j;
    printf("%u\n", sizeof(wind));
    int l = (nt+r-1);
    for(i=r;i<=l;i++)
    {
        wind[i]=0;
    }
        for(i=nt;i<=l;i++)
        {
            signal[i]=0;
        }
        printf("signal\n");
        for(i=0; i<=l; i++)
        {
            printf("%lf\n",signal[i]);
        }
        for(i=0;i<=l;i++)
        {
            rm[i]=0;
            for(j=0;j<=i;j++)
            {
                rm[i] = (rm[i]+(signal[j]*wind[i-j]) );
            }   
        }
}



void main()
{
    double a[] = {1,2,3,4};
    int i;
    int la = array_len(a);
    printf("\nc1\t%ld",la);
    double b[] = {1,1,1,1} ;
    int lb = array_len(b);
    printf("\nc2\t%ld\n",lb);
    double r[la+lb-1];
    int lr = array_len(r);
    printf("\nc3\t%ld\n",lr);

    printf("entering convolution\n");
    convolution(a,la,b,lb,r);

    for(i=0;i<(lr);i++)
    {
        printf("rm[%d]=%lf\n",i,r[i]);
    }
}

Answer 1

你试试这个！

void convolution(double *signal, int nt, double *wind, int r, double *rm)
{
    int i,j;
    printf("%u\n", sizeof(wind));  // Why you do this? this just returns the size of the pointer only
    int l = (nt+r-1);
    double one[l];
    double two[l];
    for(i=0;i<l;i++)
    {
        if (i < nt)
            one[i] = signal[i];
        else
            one[i] = 0;
        if (i < r)
            two[i] = wind[i];
        else
            two[i] = 0;
    }
    printf("signal\n");
    for(i=0; i<l; i++)
        printf("%lf\n",one[i]);
    for(i=0;i<l;i++)
    {
        rm[i]=0;
        for(j=0;j<=i;j++)
        {
            rm[i] = (rm[i]+(one[j]*two[i-j]) );
        }   
    }
}

Answer 2

可能不是您为什么将 0 视为答案的答案，但是，

for(i=r;i<=l;i++)
{
    wind[i]=0;
}

for(i=nt;i<=l;i++)
{
    signal[i]=0;
}

长度变量“l”为 (nb+r-1) 即 7，而数组传递给卷积函数时，“b”和“a”只有 4 个元素，并且仅分配内存以仅容纳 4 个元素。因此，使用您的代码，您正在跨越风和信号的阵列边界。我想知道为什么它没有崩溃。