21ic电子技术开发论坛 › 单片机与嵌入式系统 › ST MCU › 用AD采集连续两个零点
12
返回列表 发新帖我要提问

[STM32F1] 用AD采集连续两个零点

[复制链接]
iyoum 发表于 2017-9-6 22:40 | 显示全部楼层
tmp这个数据是怎么产生的？
回复 赏

举报

wwppd 发表于 2017-9-6 22:40 | 显示全部楼层
楼主最好附上完成的代码
回复 赏

举报

jkl21 发表于 2017-9-6 22:40 | 显示全部楼层
采样复位设置一个变量。
回复 赏

举报

maqianqu 发表于 2017-9-6 22:40 | 显示全部楼层
如果不想采集，将变量置位为0.
回复 赏

举报

dspmana 发表于 2017-9-6 22:40 | 显示全部楼层
采样完成就不讲tmp换成你要数据就行。
回复 赏

举报

eefas 发表于 2017-9-6 22:40 | 显示全部楼层
把flag设置为0采样不会继续的。
回复 赏

举报

pl202 发表于 2017-9-6 22:40 | 显示全部楼层
需要采样的时候将变量置位true否则为false。
回复 赏

举报

 楼主| rainic 发表于 2017-9-7 11:51 | 显示全部楼层
P70904-101731_副本.jpg
复制
  1. #include "led.h"
    

  2. #include "sys.h"
    

  3. #include "usart.h"         
    

  4. #include "adc.h"
    

  5. #include "find_peak.h"
    

  6. #include "timer.h"
    

  7. #include "math.h"
    


  8. 

  9. int sig[MAX_SIGNAL_LEN];  //采集数据 
    

  10. int peaks[MAX_SIGNAL_LEN];   //峰值点位置数据
    

  11. int troughs[MAX_SIGNAL_LEN];   // 谷值点位置数据
    

  12. int i,n1,n2,n3;
    

  13. int tmp;
    

  14. double fre;
    


  15. 

  16. int main()
    

  17.  {
    

  18.          i=0;
    

  19.          NVIC_PriorityGroupConfig(NVIC_PriorityGroup_2);
    

  20.          LED_Init ();
    

  21.          uart_init(115200);
    

  22.          Adc_Init();
    

  23.          TIM_OPMGATE_Init();
    

  24.          while(1)
    

  25.                  {
    

  26.                          tmp=Get_Adc_Average(ADC_Channel_1,10);
    

  27.                          if(tmp==3103)
    

  28.                          {
    

  29.                                  for(i=0;i<MAX_SIGNAL_LEN;i++)
    

  30.                                  {
    

  31.                                          sig[i]=Get_Adc_Average(ADC_Channel_1,10);
    

  32.                                  //sig[i]=Get_Adc(ADC_Channel_1);
    

  33.                                          printf("%d \r\n",sig[i]);
    

  34.                                  }
    

  35.                          }
    

  36.                          
    

  37.         n1 = MAX_SIGNAL_LEN;  
    

  38.         
    

  39.         process(sig,n1);  //数据预处理
    

  40.         
    

  41.         n2=find_peaks(sig,n1,peaks);  //
    

  42.         n3=find_troughs(sig,n1,troughs);   //
    


  43. 

  44.         printf("\n峰值点位置\n");
    

  45.         for(i=0;i<n2;i++)
    

  46.                          {
    

  47.                                  printf("%u ",peaks[i]+1);  //下标从1开始 
    

  48.                                  TIM_Open();
    

  49.                          }
    

  50.         printf("\谷值点位置\n");
    

  51.         for(i=0;i<n3;i++)
    

  52.                          {
    

  53.                                  printf("%u ",troughs[i]+1);
    

  54.                                  TIM_Open();
    

  55.                          }
    

  56.                  }
    

  57.  }
我的目的并不是找到峰谷值点。输入信号是类正弦信号（带噪声，周期会变化但不会突变），我想在检测到峰谷点的时候立即输出一个脉冲波。但是对峰值检测理解不深，虽然有阈值法检测程序。但是有很大的问题：首先必须检测一段数据才可以进行处理，假如数据包含了峰谷点（或许不值一个），那么这个特征点已经过去了，再加上要处理，实时性根本做不到。
方法一：所以想控制ADC采集的数据长度，以半个周期采集一次（连续两个过零点），其中必有一个峰值或者谷值点。用定时器控制ADC采样，然后findpeak（）找到峰值点位置，计算出点数乘于定时器时间就是图中的t0，由于周期不会突变，就把t0的时间给谷值点，认为这两段时间近似相等，依次类推。
方法二：最简单误差有点大，过零点把类正弦转换为方波，然后输入捕获，计算半个周期时间除于2得时间t，开启定时器定时t，时间到发送脉冲波。也刚好在峰值点或者谷值点出现脉冲。只是误差大。
很想用第一种方法，但是感觉很复杂，各位老铁给分析下。
回复 赏

举报

 楼主| rainic 发表于 2017-9-7 11:53 | 显示全部楼层
wwppd 发表于 2017-9-6 22:40
楼主最好附上完成的代码

贴了代码和目的，请老铁查看。
回复 赏

举报

 楼主| rainic 发表于 2017-9-7 11:57 | 显示全部楼层
还是这个论坛好，好多人回复。谢谢大家。
回复 赏

举报

山东电子小菜鸟 发表于 2017-9-7 12:15 | 显示全部楼层
回复 赏

举报

山东电子小菜鸟 发表于 2017-9-7 12:17 | 显示全部楼层
本帖最后由 山东电子小菜鸟 于 2017-9-7 12:18 编辑
rainic 发表于 2017-9-7 11:57
还是这个论坛好，好多人回复。谢谢大家。

int ad_buf,ad_buf2;//AD缓存
int err;//零点误差范围
int zero_flag1,zero_flag2;
int i=0,sig[]={0};
ad_buf=Get_Ad();//第一次过零检测
if(abs(ad_buf-3103)<err)
{
    zero_flag1=1;
}
if(zero_flag==1)
{
ad_buf2=Get_Ad();
if(abs(ad_buf-3103)<err)//第二次过零检测
{
   i=0;
}
else if(abs(ad_buf-3103)>=err)//中间数据
{
  seg【i】=ad_buf2;
   i++;
}


}如果觉得可以，请给分
回复 赏

举报

评分

参与人数 1威望 +2 收起 理由
rainic + 2 很给力!

查看全部评分

 楼主| rainic 发表于 2017-9-7 14:13 | 显示全部楼层
峰值检测.rar (3.45 MB, 下载次数: 1)
对一段数据处理，如果有峰谷值会输出相应的位置。程序虽然可用，但是不能满足峰谷值检测的速度。我的目的：检测到峰谷立即输出脉冲波。实时性好像很难达到

P70703-161525.jpg

测试效果，比较差劲，变频率就飘了，而且还会出现多个脉冲，老铁门有什么好的解决方案？？
回复 赏

举报

 楼主| rainic 发表于 2017-9-7 19:12 | 显示全部楼层
本帖最后由 rainic 于 2017-9-7 19:26 编辑

方法二：过零检测（用的比较器），然后TIM5输入捕获（2个上升沿），TIM2设置为单脉冲输出，TIM3作为定时器。TIM5计算连续两个过零点出时间t0，在第二个过零点开启TIM3时间为t0/2，时间到则产生中断在中断里开启TIM2输出脉冲波。 这是我的逻辑方法，老铁门感觉对不对？？
复制
  1. #include "timer.h"
    

  2. #include "led.h"
    

  3. #include "usart.h"
    

  4.   
    

  5. u16 zhouqi=300;//*0.1ms
    

  6. u16 gaodianping=100;//*0.1ms
    


  7. 

  8. void TIM2_OPMGATE_Init()
    

  9. {        
    

  10.         TIM_TimeBaseInitTypeDef TIM_TimeBaseInitStruct;
    

  11.         TIM_OCInitTypeDef TIM_OCInitStruct;
    

  12.         GPIO_InitTypeDef  GPIO_InitStruct;
    

  13.         
    

  14.         RCC_APB2PeriphClockCmd(RCC_APB2Periph_GPIOA, ENABLE);//
    

  15.         RCC_APB1PeriphClockCmd(RCC_APB1Periph_TIM2,ENABLE);//
    

  16.         RCC_APB2PeriphClockCmd(RCC_APB2Periph_AFIO,ENABLE);              
    

  17.         GPIO_InitStruct.GPIO_Mode = GPIO_Mode_AF_PP; 
    

  18.         GPIO_InitStruct.GPIO_Pin = GPIO_Pin_1;                        
    

  19.         GPIO_InitStruct.GPIO_Speed = GPIO_Speed_50MHz;
    

  20.         GPIO_Init(GPIOA,&GPIO_InitStruct);        
    

  21. //        GPIO_InitStruct.GPIO_Pin = GPIO_Pin_6;
    

  22. //        GPIO_Init(GPIOA,&GPIO_InitStruct); 
    

  23.         
    

  24.        //TIM2工作在单脉冲下
    

  25.         TIM_TimeBaseInitStruct.TIM_CounterMode = TIM_CounterMode_Up;//向上计数
    

  26.         TIM_TimeBaseInitStruct.TIM_Prescaler = 7199;
    

  27.         TIM_TimeBaseInitStruct.TIM_Period = zhouqi-1; //
    

  28.         TIM_TimeBaseInit(TIM2,&TIM_TimeBaseInitStruct);        
    

  29.         
    

  30.         TIM_SelectOnePulseMode(TIM2,TIM_OPMode_Single);       
    

  31.         TIM_SelectOutputTrigger(TIM2,TIM_TRGOSource_OC1Ref); 
    

  32.         
    

  33.         TIM_OC1PreloadConfig(TIM2,TIM_OCPreload_Enable);
    

  34.         TIM_OCInitStruct.TIM_OCMode = TIM_OCMode_PWM2;
    

  35.         TIM_OCInitStruct.TIM_OutputState = TIM_OutputState_Enable;
    

  36.         TIM_OCInitStruct.TIM_OCPolarity = TIM_OCPolarity_High;
    

  37.         TIM_OCInitStruct.TIM_Pulse = zhouqi-gaodianping;
    

  38.         TIM_OC1Init(TIM2,&TIM_OCInitStruct);                
    

  39.         TIM_Cmd(TIM2,DISABLE); //Enable the TIM Counter
    

  40.         
    

  41. }
    


  42. 

  43. void TIM2_Open()
    

  44. {
    

  45.         TIM_Cmd(TIM2,ENABLE);//使能TIM2
    

  46. }
    

  47. /*************************************************************/
    


  48. 

  49. void TIM3_Int_Init(u16 arr3,u16 psc3)
    

  50. {
    

  51.     TIM_TimeBaseInitTypeDef  TIM_TimeBaseStructure;
    

  52.         NVIC_InitTypeDef NVIC_InitStructure;
    


  53. 

  54.         RCC_APB1PeriphClockCmd(RCC_APB1Periph_TIM3, ENABLE);
    

  55.         
    

  56.         //¶¨Ê±Æ÷TIM3³õʼ»¯
    

  57.         TIM_TimeBaseStructure.TIM_Period = arr3;       
    

  58.         TIM_TimeBaseStructure.TIM_Prescaler =psc3; 
    

  59.         TIM_TimeBaseStructure.TIM_ClockDivision = TIM_CKD_DIV1; 
    

  60.         TIM_TimeBaseStructure.TIM_CounterMode = TIM_CounterMode_Up;  
    

  61.         TIM_TimeBaseInit(TIM3, &TIM_TimeBaseStructure); 
    

  62.  
    

  63.         TIM_ITConfig(TIM3,TIM_IT_Update,ENABLE ); 
    


  64. 

  65.         //ÖжÏÓÅÏȼ¶NVICÉèÖÃ
    

  66.         NVIC_InitStructure.NVIC_IRQChannel = TIM3_IRQn;  
    

  67.         NVIC_InitStructure.NVIC_IRQChannelPreemptionPriority = 0; 
    

  68.         NVIC_InitStructure.NVIC_IRQChannelSubPriority = 3; 
    

  69.         NVIC_InitStructure.NVIC_IRQChannelCmd = ENABLE; 
    

  70.         NVIC_Init(&NVIC_InitStructure);  
    


  71. 


  72. 

  73.         TIM_Cmd(TIM3, DISABLE);                                        
    

  74. }
    

  75. void TIM3_Open()
    

  76. {
    

  77.         TIM_Cmd(TIM3,ENABLE);
    

  78. }
    

  79. //中断服务函数
    

  80. void TIM3_IRQHandler(void)   
    

  81. {
    

  82.         if (TIM_GetITStatus(TIM3, TIM_IT_Update) != RESET)  //检测到中断标志
    

  83.                 {
    

  84.                         TIM2_Open();
    

  85.                 }
    

  86.                 TIM_ClearITPendingBit(TIM3, TIM_IT_Update  );  //清除标志
    

  87. }
    


  88. 


  89. 


  90. 


  91. 

  92. /***********************************************************/
    

  93. //TIM5输入捕获
    


  94. 

  95. TIM_ICInitTypeDef  TIM5_ICInitStructure;
    


  96. 

  97. void TIM5_Cap_Init(u16 arr5,u16 psc5)
    

  98. {         
    

  99.         GPIO_InitTypeDef GPIO_InitStructure;
    

  100.         TIM_TimeBaseInitTypeDef  TIM_TimeBaseStructure;
    

  101.            NVIC_InitTypeDef NVIC_InitStructure;
    


  102. 

  103.         RCC_APB1PeriphClockCmd(RCC_APB1Periph_TIM5, ENABLE);      
    

  104.          RCC_APB2PeriphClockCmd(RCC_APB2Periph_GPIOA, ENABLE);  
    

  105.         
    

  106.         GPIO_InitStructure.GPIO_Pin  = GPIO_Pin_0;  
    

  107.         GPIO_InitStructure.GPIO_Mode = GPIO_Mode_IPD; 
    

  108.         GPIO_Init(GPIOA, &GPIO_InitStructure);
    

  109.         GPIO_ResetBits(GPIOA,GPIO_Pin_0);                                               
    

  110.         
    

  111.         //³õʼ»¯¶¨Ê±Æ÷5 TIM5         
    

  112.         TIM_TimeBaseStructure.TIM_Period = arr5; 
    

  113.         TIM_TimeBaseStructure.TIM_Prescaler =psc5;       
    

  114.         TIM_TimeBaseStructure.TIM_ClockDivision = TIM_CKD_DIV1; 
    

  115.         TIM_TimeBaseStructure.TIM_CounterMode = TIM_CounterMode_Up; 
    

  116.         TIM_TimeBaseInit(TIM5, &TIM_TimeBaseStructure); 
    

  117.   
    

  118.         //初始化输入捕获参数
    

  119.         TIM5_ICInitStructure.TIM_Channel = TIM_Channel_1; //CC1S=01      
    

  120.           TIM5_ICInitStructure.TIM_ICSelection = TIM_ICSelection_DirectTI;
    

  121.           TIM5_ICInitStructure.TIM_ICPrescaler = TIM_ICPSC_DIV1;       
    

  122.           TIM5_ICInitStructure.TIM_ICFilter = 0x00;
    

  123.           TIM_ICInit(TIM5, &TIM5_ICInitStructure);
    

  124.         
    

  125.         //中断分组初始化
    

  126.         NVIC_InitStructure.NVIC_IRQChannel = TIM5_IRQn;  
    

  127.         NVIC_InitStructure.NVIC_IRQChannelPreemptionPriority = 2;  
    

  128.         NVIC_InitStructure.NVIC_IRQChannelSubPriority = 0;
    

  129.         NVIC_InitStructure.NVIC_IRQChannelCmd = ENABLE;
    

  130.         NVIC_Init(&NVIC_InitStructure); 
    

  131.         
    

  132.         TIM_ITConfig(TIM5,TIM_IT_Update|TIM_IT_CC1,ENABLE);     
    

  133.         
    

  134.            TIM_Cmd(TIM5,ENABLE );       
    

  135.    
    

  136. }
    


  137. 


  138. 

  139. u8  TIM5CH1_CAPTURE_STA=0;     //输入捕获状态                                              
    

  140. u16        TIM5CH1_CAPTURE_VAL;   //输入捕获值    
    

  141. u32 temp=0;
    

  142. //中断程序     
    

  143. void TIM5_IRQHandler(void)
    

  144. {        
    

  145.          if((TIM5CH1_CAPTURE_STA&0X80)==0)//还未成功捕获    
    

  146.         {          
    

  147.                 if (TIM_GetITStatus(TIM5, TIM_IT_Update) != RESET)
    

  148.                  
    

  149.                 {            
    

  150.                         if(TIM5CH1_CAPTURE_STA&0X40)//已经捕获到高电平
    

  151.                         {
    

  152.                                 if((TIM5CH1_CAPTURE_STA&0X3F)==0X3F)//高电平时间过长
    

  153.                                 {
    

  154.                                         TIM5CH1_CAPTURE_STA|=0X80;//标记捕获成功
    

  155.                                         TIM5CH1_CAPTURE_VAL=0XFFFF;
    

  156.                                 }else TIM5CH1_CAPTURE_STA++;
    

  157.                         }         
    

  158.                 }
    

  159.         if (TIM_GetITStatus(TIM5, TIM_IT_CC1) != RESET)//²¶»ñ1·¢Éú²¶»ñʼþ
    

  160.                 {        
    

  161.                         if(TIM5CH1_CAPTURE_STA&0X40)                               
    

  162.                         {
    

  163.                 
    

  164.                                 TIM5CH1_CAPTURE_STA|=0X80;               
    

  165.                                 TIM5CH1_CAPTURE_VAL=TIM_GetCapture1(TIM5);
    

  166.                                    TIM_OC1PolarityConfig(TIM5,TIM_ICPolarity_Rising); //设置为上升沿捕获
    

  167.                         }else                                                               
    

  168.                         {
    

  169.                                 TIM5CH1_CAPTURE_STA=0;                     
    

  170.                                 TIM5CH1_CAPTURE_VAL=0;
    

  171.                                  TIM_SetCounter(TIM5,0);
    

  172.                                 TIM5CH1_CAPTURE_STA|=0X40;             
    

  173.                                    TIM_OC1PolarityConfig(TIM5,TIM_ICPolarity_Rising);                //设置为上升沿捕获
    

  174.                         }                    
    

  175.                 }                                                                                    
    

  176.          }
    

  177.         
    

  178.         if(TIM5CH1_CAPTURE_STA&0X80)
    

  179.                 {
    

  180.                         //TIM3_Open();
    

  181.                         temp=TIM5CH1_CAPTURE_STA&0X3F;
    

  182.                         temp*=65536;//溢出时间总和
    

  183.                         temp+=TIM5CH1_CAPTURE_VAL;//两个零点时间
    

  184.                         printf("HIGH:%d us\r\n",temp);//打印输出
    

  185.                         TIM5CH1_CAPTURE_STA=0;//清零 下次捕获用
    

  186.                 }
    

  187.  
    

  188.     TIM_ClearITPendingBit(TIM5, TIM_IT_CC1|TIM_IT_Update); //清除中断标志
    

  189.  
    

  190. }
回复 赏

举报

 楼主| rainic 发表于 2017-9-7 19:29 | 显示全部楼层
本帖最后由 rainic 于 2017-9-7 19:31 编辑
复制
  1. #include "led.h"
    

  2. #include "key.h"
    

  3. #include "sys.h"
    

  4. #include "usart.h"
    

  5. #include "timer.h"
    


  6. 

  7.  extern u32 temp;
    


  8. 

  9.  int main(void)
    

  10.  {         
    

  11.         NVIC_PriorityGroupConfig(NVIC_PriorityGroup_2);        
    

  12.         uart_init(115200);      
    

  13.          LED_Init();        
    

  14.         TIM2_OPMGATE_Init();
    

  15.         TIM3_Int_Init(temp/2,72-1);
    

  16.          TIM5_Cap_Init(0XFFFF,72-1);       
    

  17.          //TIM2_Open();
    

  18.            while(1)
    

  19.         {}}
回复 赏

举报

zhb3373 发表于 2017-9-7 21:30 | 显示全部楼层
学习了学习了，好好
回复 赏

举报

12
返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 在线客服 返回列表 返回顶部
快速回复 在线客服 返回列表 返回顶部