test10.c

来自「5402开发板附带的源程序」· C语言 代码 · 共 73 行

#include <math.h>

unsigned int N1,N2,n;
int m,i,k;
float sum_real,sum_imag,t_real,t_imag;
float x_real[500],x_imag[500];
float h_real[500],h_imag[500];
float y_real[500],y_imag[500];

main()
{	
   N1=8;                           /* x 长度*/
   N2=8;                           /* h 长度*/
   n=N1+N2-1;                      /* 输出 y 长度*/  

   for(i=0;i<500;i++)              /* 初始化数组 */
   {
     x_real[i]=0;
     x_imag[i]=0;
     h_real[i]=0;
     h_imag[i]=0;
     y_real[i]=0;
     y_imag[i]=0;
   }

   for(i=0;i<n;i++)              
   {
      if(i<N1)
      {
         x_real[i]=exp(-i);
         x_imag[i]=0;
      }
      else
      {
         x_real[i]=0;
         x_imag[i]=0;
      }
   }

   for(i=0;i<n;i++)
   {
      if(i<N2)
      {
         h_real[i]=1;
         h_imag[i]=0;
      }
      else
      {
         h_real[i]=0;
         h_imag[i]=0;
      }
   }

   for(i=0;i<n;i++)
   {
      sum_real=0;
      sum_imag=0;
      t_real=0;
      t_imag=0;

      for(k=0;k<=i;k++)
      {
         t_real=x_real[k]*h_real[i-k]-x_imag[k]*h_imag[i-k];
         t_imag=x_real[k]*h_imag[i-k]+x_imag[k]*h_real[i-k];
         sum_real=sum_real+t_real;
         sum_imag=sum_imag+t_imag;
      }

      y_real[i]=sum_real;
      y_imag[i]=sum_imag;             
   }
}