私はcでxcodeを使用していますが、次のコードで奇妙な問題が発生しています。
アルゴリズムがconv関数のy[k]= 0.0行に到達し、x配列が0で埋められると、理由を理解する必要があります。私はすでにこのエラーを独自の方法で修正しています。main関数で2行目をコメントし、3行目のコメントを外しても、問題は発生しません(以下を参照)。generateSquareSignalで配列を再割り当てすることは知っていますが、これがこのバグの理由ではないと思います。
int length = 100;
double *output = (double*) malloc(10 * length * sizeof(double));
//double *output;
output = generateSquareSignal(length);
double *input1 = (double*) malloc(length * sizeof(double));
double *input2 = (double*) malloc(length * sizeof(double));
for (int i = 0; i < length; i++) {
input2[i] = output[i];
input1[i] = output[i];
//printf("-%d=%lf\n",i ,input1[i]);
}
conv(input1, length, input2, length, output, 2 * length );
double* generateSquareSignal(int length) {
printf("double* generateSquareSignal(int length)\n");
double *signal = (double*) malloc(length * sizeof(double));
int length_period = length / kPeriodSignal;
for(int i=0; i < length; i++) {
if (i % (length_period) < (length_period / 2)) {
signal[i] = 1.0;
} else {
signal[i] = -1.0;
}
//printf("%d. - %lf\n", i, signal[i]);
}
return signal;
}
void conv( double *x, int N_signal,
double *h, int N_kernel,
double *y, int N_output) {
int k;
for(k = 0; k < N_signal + N_kernel; k++) {
y[k] = 0.0;
}
for (int i = 0; i < N_signal; i++) {
printf("%lf-%lf\n", x[i], y[i]);
}
for(int i = 0; i < N_signal; i++) {
for(int j = 0; j < N_kernel; j++) {
double xx = x[i];
double hh = h[j];
double yy = y[i + j];
y[i + j] += x[i] * h[j];
//printf("=%lf\n", y[i + j]);
}
}
}