原创 示波器制作系列1：正弦插值算法

/*
pbuf    数据缓冲空间
buf_num 数据缓冲的大小     缓冲空间的大小 最好是WAVE_LEN的除数 按照正弦插值的理论算法，每个周期2.5采样点就可以重构一个正弦波形
color    线条的颜色
AD_CLK  AD采集时钟频率

///////////
数据点不能少于    SINP_LEN 个
*/
void GUI_Wave_SinP(uint8 *pbuf,uint16 buf_num,uint16 color)
{
    uint16  i = 0 ;     
    uint16  j = 0 ;
    float   sum      = 0;
    float   temp  = 0;
    uint8 data_buf[WAVE_LEN] = {0};           //存放插值以后的数据
    float sa_buf  [SINP_LEN] = {0};        // 正弦内插函数求得的系数

    uint16 interval = WAVE_LEN / buf_num ;
   
    for(i = 1 ; i < WAVE_LEN ; i++)      
    {
        temp = 0 ;      
        sum  = 0 ;        //不能少了啊
        for(j = 0 ; j < SINP_LEN ; j++) //计算内插系数
        {
         // sa_buf[j] =   sin(PI * ADC_CLK * (1/ ( ADC_CLK * interval) * i - j)       //原始计算式
            temp =     (float) ( PI *  ( (float) i / (float)interval - (float)j )) ;
            if (temp == 0)
                sa_buf[j] = 1.0000;
            else
                   sa_buf[j] =(float)((float) sin(temp) / (float)temp  );                                    //简化后的式子
        }

        //计算插值以后每个点的数据

    /*    if( i <= SINP_LEN )                    //如果计算的点小于SINP_LEN的一半，采用数据缓冲区的前60个点计算
        {
            for(j = 0 ; j < SINP_LEN ; j++)
                  sum +=  (float)sa_buf[j] * (float)(*(pbuf + j)) ;
        }
        else if( i >=  buf_num - SINP_LEN )   //如果计算的点大于buf_num - SINP_LEN / 2，采用数据缓冲区的后60个点计算
        {
            for(j = 0 ; j < SINP_LEN ; j++)
                  sum +=  (float)sa_buf[j] * (float)(*(pbuf + buf_num - SINP_LEN + j)) ;
        }
        else
        {
               for(j = 0 ; j < SINP_LEN ; j++)    //其余情况采用以这个点为中心的  SINP_LEN个点
                 sum +=  (float)sa_buf[j] * (float)(*(pbuf + i + j - (uint8)(SINP_LEN / 2) )) ;
        }   
          */
       for(j = 0 ; j < SINP_LEN ; j++)
               sum +=  (float)sa_buf[j] * (float)(*(pbuf + j)) ;
        data_buf  =   (uint8) sum ;
       }
   
    GUI_Wave(data_buf,color);
 
}