/*-------------------------------------------*
*                   離散卷積                *
*-------------------------------------------*/

#include<stdio.h>

float convolute(float* f, /*離散單邊信號f的序列表示*/ 
        float* g, /*離散單邊信號g的序列表示*/
        unsigned int k, /*卷積的自變量*/
        const int length /*信號的序列長度,本例子假設為10*/
        )
{
  unsigned int n = 0;
  float sum = 0; /*累加的結果，最終作為返回值返回*/
  if( k > length )
  {
  k = length;
  }

  for( n = 0 ; n != k ; n++ ) /*卷積公式*/
  {
    sum+=( f[n]*g[k-n] );
	
  }
  return sum;
}
int main()
{
float f[10] = {3,5,4,6,9,8,7,4,5,12};  /*離散單邊信號f的序列表示*/
float g[10] = {2,5,1,3,6,8,7,8,2,22};  /*離散單邊信號g的序列表示*/

float convolution = 0; /*該變量用于存儲卷積的結果*/

unsigned int k = 0;   /*卷積的自變量*/

printf( "lease input the value of k to convolute\n" );/*輸入卷積自變量*/

scanf( "%d", &k );

convolution = convolute( f, g, k ,10);  /*調用函數計算卷積公式*/

printf( "\nThe result: %f", convolution ); /*輸出結果*/


return 0;
}