這是我的代碼,用於查找兩個信號的卷積,但是我的輸出每次都變成零,任何人都可以解釋我的代碼中的錯誤? 我試圖compling,但我無法findout是什麼問題我的整個代碼是正確的 由於事先這個程序中的錯誤是什麼?
#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]);
}
}
請嘗試縮小問題 – rethab
是的,我能夠編譯它 –
你的輸出不是零,無論如何,我得到了一些異常的結果,這是由於訪問數組越界! –