OpenEdv-开源电子网

 找回密码
 立即注册
正点原子全套STM32/Linux/FPGA开发资料,上千讲STM32视频教程免费下载...
查看: 236|回复: 3

FFT求取相位不稳定的问题如何解决

[复制链接]

3

主题

6

帖子

0

精华

初级会员

Rank: 2

积分
80
金钱
80
注册时间
2024-9-1
在线时间
12 小时
发表于 2024-10-22 17:55:35 | 显示全部楼层 |阅读模式
2金钱
主要流程:(ADC设置:RCC_PCLK2_Div6;ADC_RegularChannelConfig(ADC1, ADC_Channel_5 , 1, ADC_SampleTime_41Cycles5);     A和B信号的周期为500微秒)
1、定时器T2=1000ms时开启信号A,T2=1005ms开启ADC + DMA 连续采集信号A信号值 DMAbuf[256]个点,在DMA采集完成中断中 ,关闭ADC+DMA,并把DMAbuf[256]的数据复制给A[256]。

2、定时器T2=1020ms时关闭信号A,开启信号B,T2=1025ms开启ADC + DMA 连续采集信号B信号值 DMAbuf[256]个点,在DMA采集完成中断中 ,关闭ADC+DMA,并把DMAbuf[256]的数据复制给B[256]。


3、定时器T2=1040ms时关闭信号B,随后对A[256]和B[256]进行FFT变换,然后求取两个信号变换后幅值最大的相位。


附上4个图,


求取FFT变换后最大幅值分量的相位部分的程序:

float GetPowerMag(long * Inbuff)
{
      
        float Outbuuf;
        signed short lX,lY;
        float X,Y,Mag;
        unsigned long MAG[128] = {0};
        unsigned long MaxMag = 0,MaxMagBuf = 0;
        float phrase = 0;
        unsigned short i,j;
        for(i=0;i<128;i++)
        {
                lX = (Inbuff << 16) >> 16;        //实部
                lY = (Inbuff >> 16);                        //虚部
                X = 256 * ((float)lX) / 32768;
                Y = 256 * ((float)lY) / 32768;
                Mag = sqrt(X * X + Y * Y) / 256;                //求取幅值
                if(i == 0)
                {
                        MAG = (unsigned long)(Mag * 32768);   //直流分量
                }
                else
                {
                        MAG = (unsigned long)(Mag * 65536);                //除直流分量外的幅值
                }
        }
      
        for(j=1;j<128;j++)
        {
                if(MaxMag < MAG[j])                        //判断最大幅值的数据
                {
                        MaxMag = MAG[j];
                        MaxMagBuf = Inbuff[j];  //提取最大幅值对应的FFT变换后的数据
                }
        }
                lX = (MaxMagBuf << 16) >> 16;   
                lY = (MaxMagBuf >> 16);                       
                phrase =atan2(lY,lX);                                //相位弧度
                Outbuuf = phrase*180.0f/3.1415926f;  //相位角度换算
        return Outbuuf;
}




通过打印出的数据发现处理完后得到AB信号相位的数据不固定,在某一范围内浮动没有规律,这种情况怎么处理?

2.png
1.png
4.png
3.png
正点原子逻辑分析仪DL16劲爆上市
回复

使用道具 举报

10

主题

3280

帖子

1

精华

论坛元老

Rank: 8Rank: 8

积分
8182
金钱
8182
注册时间
2020-5-11
在线时间
3694 小时
发表于 2024-10-23 09:49:08 | 显示全部楼层
ADC采到的数据画成曲线,看是否有抖动。
专治疑难杂症
回复

使用道具 举报

3

主题

6

帖子

0

精华

初级会员

Rank: 2

积分
80
金钱
80
注册时间
2024-9-1
在线时间
12 小时
 楼主| 发表于 2024-10-23 10:57:00 | 显示全部楼层
本帖最后由 18633846727 于 2024-10-23 10:58 编辑
LcwSwust 发表于 2024-10-23 09:49
ADC采到的数据画成曲线,看是否有抖动。

A信号还好,B信号有波动。

02.png
01.png
03.png
回复

使用道具 举报

10

主题

3280

帖子

1

精华

论坛元老

Rank: 8Rank: 8

积分
8182
金钱
8182
注册时间
2020-5-11
在线时间
3694 小时
发表于 2024-10-23 17:06:10 | 显示全部楼层
18633846727 发表于 2024-10-23 10:57
A信号还好,B信号有波动。

我的意思时这次的波形与上次的波形在时间轴上看(比如看波峰位置),是否有左右抖动,或者看每次的波形能否重叠。
专治疑难杂症
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则



关闭

原子哥极力推荐上一条 /2 下一条

正点原子公众号

QQ|手机版|OpenEdv-开源电子网 ( 粤ICP备12000418号-1 )

GMT+8, 2024-11-22 06:28

Powered by OpenEdv-开源电子网

© 2001-2030 OpenEdv-开源电子网

快速回复 返回顶部 返回列表