探索者DSP_FFT测试

hi我歌月徘徊

使用了探索者实验47——2DSP FFT测试，
  测试完例程ok，然后我又AD采样一路直流电压，然后DMA发送1024个数据存起来。
然后根据原来的例程，把原来要进行fft的入口参数进行修改
[mw_shl_code=applescript,true]for(i=0;i<FFT_LENGTH;i++)
                        {
                                 fft_inputbuf[2*i]=ADC1ConvertedValue;
                                 fft_inputbuf[2*i+1]=0;
                        }[/mw_shl_code]

ADC1ConvertedValue[1024]是我AD采样结果
其他均不变 然后下载进开发板，按键key0 串口返回fft结果
1024 point FFT runtime:0.585ms
FFT Result:
fft_outputbuf[0]:2045429.000000
fft_outputbuf[1]:1230.973755
fft_outputbuf[2]:135.815155
fft_outputbuf[3]:1285.064697
fft_outputbuf[4]:1249.406982
fft_outputbuf[5]:172.963455
fft_outputbuf[6]:1081.007080
fft_outputbuf[7]:2533.878418
fft_outputbuf[8]:1983.940918
fft_outputbuf[9]:1012.018494
fft_outputbuf[10]:944.124329

因为我AD采样的是直流1.8v电压，根据分析：第 0 点，即直流分量，其 FFT 变换后的模值应该是原始信号幅值的 N 倍， N=1024，然后其他点，模值应该是原始信号幅值的 N/2 倍，
计算一下串口返回的数据： 第0点  2045429/1024=1997.489   直流电压的直流分量不就应该反映电压的真实值吗？？？  还有直流电压频率为0 其他点不都应该为0吗 为什么还有数值存在呢？？？  刚接触fft  求解

采样的AD值转化为电压之后
[mw_shl_code=applescript,true]fft_inputbuf[2*i]=(float)ADC1ConvertedValue*3.3/4096;
                 fft_inputbuf[2*i+1]=0;[/mw_shl_code]fft之后的数据
fft_outputbuf[0]:2108.405029  2108.405029  0.000000   -131.702850
fft_outputbuf[1]:0.297549  0.297549  0.000286   -131.702850
fft_outputbuf[2]:0.515750  0.045758  0.513716   -131.702850
fft_outputbuf[3]:0.770250  0.241686  -0.731350   -131.702850
fft_outputbuf[4]:0.442877  0.067056  -0.437771   -131.702850
fft_outputbuf[5]:0.669347  0.479659  -0.466854   -131.702850
fft_outputbuf[6]:0.261266  -0.155807  -0.209725   -131.702850
fft_outputbuf[7]:0.859252  -0.002719  -0.859248   -131.702850
fft_outputbuf[8]:0.636206  0.437347  -0.462045   -131.702850
fft_outputbuf[9]:0.639701  0.504262  0.393620   -131.702850
fft_outputbuf[10]:0.548882  0.489041  -0.249219   -131.702850

直流分量为2108/1024=2.058v    实测电压为2.06v      这样fft才对嘛~~

hi我歌月徘徊

直流的问题解决了 ，采样交流信号还是遇到了问题       使用电压传感器 将220v电压转换为2ma电流 之后加入1.62v（图上标的1.5v）直流电压将采样的交流信号幅度抬高   经过运放送给ADC
因为市电为50hz，所以根据采样定理设为256点的fft（因为只有16/64/256/1024/4096），其他程序不用调整。
串口返回的数据如下
                        模值           real              imag         之前算的一个相位（之前测试fft_inputbuf[2*i]=2+3*arm_cos_f32(2*PI*50*i/FFT_LENGTH-PI*30/180)；频率50hz，相位为-30，fft后数据都正确）
fft_outputbuf[0]:487.360718  487.360718  0.000000   -12.062287
fft_outputbuf[1]:7.286572  5.948340  4.208488   -12.062287
fft_outputbuf[2]:4.363874  -4.285675  0.822431   -12.062287
fft_outputbuf[3]:2.730595  1.053548  -2.519164   -12.062287
fft_outputbuf[4]:1.508674  0.127408  1.503284   -12.062287
fft_outputbuf[5]:1.111992  -0.665747  -0.890678   -12.062287
fft_outputbuf[6]:1.298024  0.632711  -1.133377   -12.062287
fft_outputbuf[7]:0.381704  0.047216  0.378773   -12.062287
fft_outputbuf[8]:0.573807  -0.557957  -0.133936   -12.062287
fft_outputbuf[9]:0.864974  0.793240  -0.344893   -12.062287
fft_outputbuf[10]:0.728437  -0.364859  0.630475   -12.062287

按道理这里面的直流分量应该为1.62v吧~~~    但是487/256=1.90v   还有第一个数据应该为基波电压吧 7.2肯定不对嘛。。然后是一次谐波 二次谐波。。。又不懂了   。。。

正点原子

帮顶