calc_fft

Posted 2022-03-29 firstparke

/******************************************************0x005E******************************************************/
/*
 *  sw  三相电流不平衡度
 */
void CurrentUnbanceShell(slot_enum _no)

    TSlotState SlotInfo = .Io1Set = 1, .Io4Set = 0, .Io4Out = 1,  .VddVol = 1, .AD0En = 1,  .AD1En = 1;
    unsigned char ADId[3][2] =   15, 13, 16, 4, 11, 3 ; 
    uint16_t Input_u16[256];
    float Input_f32[256], Kp = 19.96008;//1000/50.1;
    FFTRunData _fft; SUnbalance _sun;
    float Am[256]; float Phase[256]; float AmHarm[10]; float PhaseHarm[10]; float FreqHarm[10];
    float abc[3][2], unbalanceBasic[3], unbalZero, unbalNega;    
    SensorSlotCreate(_no, &SlotInfo);
    osDelay(500);
/* FFT算出基波和相位 */        
    GPIO4[_no]->BSRR = 1<<  (GPIOPos4[_no] + 16);//ADC_AII0 - ADC_AII2
    osDelay(500);
    GeneralGetADSamples(ADId[_no][0], 500, Input_u16, 256);
    for(int i = 0; i < 256; i++)  Input_f32[i] = (float)(Kp*(AD_CLAC(Input_u16[i]) - 1650.0));
    
    for(int i=0;i<256;i++)//生成信号序列
    
        Input_f32[i]=230*arm_sin_f32(2*PI*i*49.5/2000+5*PI/3); //生成实部
    
    
    FFT_NEW(&_fft, 2000, Input_f32, Am, Phase, AmHarm, PhaseHarm, FreqHarm);
    FFT_RUN(_fft);
    abc[0][0] = AmHarm[0]*0.70710678;//基波有效值
    abc[0][1] = PhaseHarm[0];//基波相位
    MSG_LOG("Ia_rms is %f, Ia_phase is %f\r\n", abc[0][0], abc[0][1]);
    
    GeneralGetADSamples(ADId[_no][1], 500, Input_u16, 256);
    for(int i = 0; i < 256; i++)  Input_f32[i] = (float)(Kp*(AD_CLAC(Input_u16[i]) - 1650.0));

    for(int i=0;i<256;i++)//生成信号序列
    
        Input_f32[i]=228*arm_sin_f32(2*PI*i*49.5/2000+PI/3); //生成实部
    
    
    FFT_NEW(&_fft, 2000, Input_f32, Am, Phase, AmHarm, PhaseHarm, FreqHarm);
    FFT_RUN(_fft);
    abc[2][0] = AmHarm[0]*0.70710678;
    abc[2][1] = PhaseHarm[0];
    MSG_LOG("Ic_rms is %f, Ic_phase is %f\r\n", abc[2][0], abc[2][1]);

    
    GPIO4[_no]->BSRR = 1<<  GPIOPos4[_no];//ADC_AII1 - ADC_AII3
    osDelay(500);
    GeneralGetADSamples(ADId[_no][0], 500, Input_u16, 256);
    for(int i = 0; i < 256; i++)  Input_f32[i] = (float)(Kp*(AD_CLAC(Input_u16[i]) - 1650.0));
    
    for(int i=0;i<256;i++)//生成信号序列
    
        Input_f32[i]=230*arm_sin_f32(2*PI*i*49.5/2000+PI); //生成实部
    
    
    FFT_NEW(&_fft, 2000, Input_f32, Am, Phase, AmHarm, PhaseHarm, FreqHarm);
    FFT_RUN(_fft);
    abc[1][0] = AmHarm[0]*0.70710678;
    abc[1][1] = PhaseHarm[0];
    MSG_LOG("Ib_rms is %f, Ib_phase is %f\r\n", abc[1][0], abc[1][1]);

/* 对称分量法算出正序负序和零序分量 */        
   _sun.a.real = abc[0][0]*cos(abc[0][1]*0.0174533);//PI/180
   _sun.a.image =  abc[0][0]*sin(abc[0][1]*0.0174533);
   
   _sun.b.real = abc[1][0]*cos(abc[1][1]*0.0174533);
   _sun.b.image = abc[1][0]*sin(abc[1][1]*0.0174533);
   
   _sun.c.real = abc[2][0]*cos(abc[2][1]*0.0174533);
   _sun.c.image = abc[2][0]*sin(abc[2][1]*0.0174533);   
   SymmetryWeight(&_sun);

/* 按国标算出零序不平衡和负序不平衡 */
    unbalanceBasic[1] = sqrt(_sun.z1.real*_sun.z1.real + _sun.z1.image*_sun.z1.image);
    unbalanceBasic[2] = sqrt(_sun.z2.real*_sun.z2.real + _sun.z2.image*_sun.z2.image);        
    unbalanceBasic[0] = sqrt(_sun.z0.real*_sun.z0.real + _sun.z0.image*_sun.z0.image);
    
    unbalNega = unbalanceBasic[2]/unbalanceBasic[1];
    unbalZero = unbalanceBasic[0]/unbalanceBasic[1];

   MSG_LOG("unbalance Zero is %f, unbalance Nega is %f\r\n", unbalZero, unbalNega);
   
   

    
    
    
    
    TSlotState SlotInfo = .Io1Set = 1, .Io4Set = 0, .Io4Out = 1,  .VddVol = 1, .AD0En = 1,  .AD1En = 1;
    unsigned char ADId[3][2] =   15, 13, 16, 4, 11, 3 ; 
    uint16_t Input_u16[256*2], Input_u16_A[256], Input_u16_C[256],fs = 2000;;
    float Input_f32[256], Kp = 0.1996008;//1000/50.1;
    FFTRunData _fft; SUnbalance _sun;
    float Am[256]; float Phase[256]; float AmHarm[10]; float PhaseHarm[10]; float FreqHarm[10];
    float abc[3][3], unbalanceBasic[3], unbalZero, unbalNega;    
    SensorSlotCreate(_no, &SlotInfo);
    osDelay(500);
/* FFT算出基波和相位 */        
    GPIO4[_no]->BSRR = 1<<  (GPIOPos4[_no] + 16);//ADC_AII0 - ADC_AII2
    osDelay(500);
    GeneralGetADSamples(&ADId[_no][0], 2, 500, Input_u16, 256);
#if 1
    for(int i= 0; i<256; i++)
    
        Input_u16_A[i] = Input_u16[2*i];
        Input_u16_C[i] = Input_u16[2*i+1];
    
#endif    
    for(int i = 0; i < 256; i++)  Input_f32[i] = (float)(Kp*(AD_CLAC(Input_u16_A[i]) - 1650.0));
    
    FFT_NEW(&_fft, fs, Input_f32, Am, Phase, AmHarm, PhaseHarm, FreqHarm);
    FFT_RUN(_fft);
    abc[0][0] = AmHarm[0]*0.70710678;//基波有效值
    abc[0][1] = PhaseHarm[0];//基波相位
    abc[0][2] = FreqHarm[0];//基波频率
    MSG_LOG("Ia_rms is %f, Ia_phase is %f, Ia_freq is %f\r\n", abc[0][0], abc[0][1], abc[0][2]);

    for(int i = 0; i < 256; i++)  Input_f32[i] = (float)(Kp*(AD_CLAC(Input_u16_C[i]) - 1650.0));
    
    FFT_NEW(&_fft, fs, Input_f32, Am, Phase, AmHarm, PhaseHarm, FreqHarm);
    FFT_RUN(_fft);
    abc[2][0] = AmHarm[0]*0.70710678;
    abc[2][1] = PhaseHarm[0];
    abc[2][2] = FreqHarm[0];//基波频率
    MSG_LOG("Ic_rms is %f, Ic_phase is %f, Ic_freq is %f\r\n", abc[2][0], abc[2][1], abc[2][2]);