fft计算P.Q.PF_fft计算电压有效值

void FFT_UIC()
{
	u8 i;
	fCosU=0;
	fSinU=0;
	fCosI=0;
	fSinI=0;
	fCos=0;
	fSin=0;
//	T_ON();
	memset(fft_inputbuf,0,FFT_LENGTH);
	memset(fft_outputbuf,0,FFT_LENGTH);
	for(i=0;i<FFT_LENGTH;i++)//生成信号序列
	{
		 fft_inputbuf[2*i]=ADC_3[i];	//生成输入信号实部
		 fft_inputbuf[2*i+1]=0;//虚部全部为0
	}
	arm_cfft_radix2_f32(&scfft,fft_inputbuf);	//FFT计算（基4）
	arm_cmplx_mag_f32(fft_inputbuf,fft_outputbuf,FFT_LENGTH);	//把运算结果复数求模得幅值 45us
	
	fftu1=fft_outputbuf[1]/103220;/*基波有效值 vpp*/
	fCosU=fft_inputbuf[2]/fft_outputbuf[1];//实部
	fSinU=fft_inputbuf[3]/fft_outputbuf[1];//虚部
		
	memset(fft_inputbuf,0,FFT_LENGTH);
	memset(fft_outputbuf,0,FFT_LENGTH);
	
	for(i=0;i<FFT_LENGTH;i++)//生成信号序列
	{
		 fft_inputbuf[2*i]=ADC_7[i];	//生成输入信号实部
		 fft_inputbuf[2*i+1]=0;//虚部全部为0
	}
	arm_cfft_radix2_f32(&scfft,fft_inputbuf);	//FFT计算（基4）
	arm_cmplx_mag_f32(fft_inputbuf,fft_outputbuf,FFT_LENGTH);	//把运算结果复数求模得幅值 
	
	ffti1=fft_outputbuf[1]/103210;/*基波有效值 vpp*/
	fCosI=fft_inputbuf[2]/fft_outputbuf[1];//实部
	fSinI=fft_inputbuf[3]/fft_outputbuf[1];//虚部
	
	fCos=fCosU*fCosI + fSinU*fSinI;
	fSin=fSinU*fCosI - fCosU*fSinI;
	SC=rms3*rms6;
	PC=SC*fCos;
	QC=SC*fSin;
	PFC=PC/SC;
//	T_OFF();
}

void FFT_HRI()
{
	u16 i;
	fCosU=0;
	fSinU=0;
	fCosI=0;
	fSinI=0;
	fCos=0;
	fSin=0;
//	T_ON();
	memset(fft_inputbuf,0,FFT_LENGTH);
	memset(fft_outputbuf,0,FFT_LENGTH);
	for(i=0;i<FFT_LENGTH;i++)//生成信号序列
	{
		 fft_inputbuf[2*i]=ADC_9[i];	//生成输入信号实部
		 fft_inputbuf[2*i+1]=0;//虚部全部为0
	}
	arm_cfft_radix2_f32(&scfft,fft_inputbuf);	//FFT计算（基2）
	arm_cmplx_mag_f32(fft_inputbuf,fft_outputbuf,FFT_LENGTH);	//把运算结果复数求模得幅值 45us
	
	fft_hria=fft_outputbuf[1]/524.626;/*基波有效值 vpp*/

		
	memset(fft_inputbuf,0,FFT_LENGTH);
	memset(fft_outputbuf,0,FFT_LENGTH);
	
	for(i=0;i<FFT_LENGTH;i++)//生成信号序列
	{
		 fft_inputbuf[2*i]=ADC_10[i];	//生成输入信号实部
		 fft_inputbuf[2*i+1]=0;//虚部全部为0
	}
	arm_cfft_radix2_f32(&scfft,fft_inputbuf);	//FFT计算（基2）
	arm_cmplx_mag_f32(fft_inputbuf,fft_outputbuf,FFT_LENGTH);	//把运算结果复数求模得幅值 
	
	fft_hrib=fft_outputbuf[1]/526.628;/*基波有效值 vpp*/
	
}
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86

可由公式 atan2(imag, real) 得到相位值。
其中：imag为虚数值，real为实数值，atan2为求反正切值。