【创龙TMS320C6748开发板试用】用CCS画出信号时域图&在对信号进行傅里叶变换时遇到的问题

                                                      图 6
           2、傅里叶变换算法，贴上a[k],b[k]，c[k]的代码，注意以下代码并不能实现需要的傅里叶变换，主要问题是：为什么得到的a[k]和b[k]均为0，贴出来是希望哪位能够指正以下。
  #define N 256
#defin #define PF(x)    ((x)*(x))
double *compute_ak()
{
        int n=256;
        float dt=0.02;
        double *a=(double *)malloc(128*sizeof(double));
        int k,j;

        int h;
                         for(h=0;h<128;h++)
                         {
                                 a[h]=0.0;
                         }
                         double *sin_ak=(sin_fun(n,dt));
        for(k=0;k<128;k++)
        {
                for(j=0;j<256;j++)
                {
                        a[k]=a[k]+(sin_ak[j])*(2/N)*cos(2*pi*k*j/N);
                }
        }

    return a;
}

double *compute_bk()
{
            int n=256;
                 float dt=0.02;
                 double *b=(double *)malloc(128*sizeof(double));
                 int k,j;
                 double *sin_bk=(sin_fun(n,dt));
                 int h;
                 for(h=0;h<128;h++)
                 {
                         b[h]=0.0;
                 }
                 for(k=0;k<128;k++)
                 {
                         for(j=0;j                          {
                                 b[k]=b[k]+(sin_bk[j])*(2/N)*sin(2*pi*k*j/N);
                         }
                 }

             return b;
}

double *compute_ck()
{
        double *c=(double *)malloc(128*sizeof(double));
        double *a=compute_ak();
        double *b=compute_bk();
        int k;
        for(k=0;k         {
                        c[k]=sqrt(PF(a[k])+PF(b[k]));
        }
        return c;
}

                        
                                                  
           3、上述希望得到的c即为经傅里叶变换后所得到的谱值(由于a和b均为0，得到的c也为0)。


            附MATLAB的傅里叶变换程序以及效果图
        

                           

                                 

1

7.png (29.97 KB, 下载次数: 7)

8.png (29.44 KB, 下载次数: 5)

9.png (21.9 KB, 下载次数: 10)