派克变换VC++源码（附文档）

派克变换的VC 实现

#include “pch.h“
#include “ParkTransform.h“
#include 
#include “fftpf.h“
using namespace std;
using namespace splab;
ParkTransform::ParkTransform（）
{
}


ParkTransform::~ParkTransform（）
{
}

void ParkTransform::park_transform（double * a double * b double * c int length double* id double* iq）
{
	for （int i = 0; i < length; i++） {
		id[i] = （sqrt（2） / sqrt（3））*a[i] - （1 / sqrt（6））*b[i] - （1 / sqrt（6））*c[i];
		iq[i] = （1 / sqrt（2））*b[i] - （1 / sqrt（2））*c[i];
	}
}
double getValue（complex num） {
	return sqrt（std::pow（num.real（） 2） + std::pow（num.imag（） 2））;
}
void ParkTransform::extended_parker_vector（double * id double * iq int length int Fs double * amplitudedouble* frequency
）
{
	double* park = new double[length];
	memset（park 0 sizeof（double）*length）;
	for （int i = 0; i < length; i++） {
		park[i] = sqrt（id[i] * id[i] + iq[i] + iq[i]）;
	}
	Vector> sn Sk;
	sn.resize（length）;
	Sk.resize（length）;

	for （int i = 0; i < length; ++i）
		sn[i] = park[i];
	/******************************** [ FFT ] ********************************/
	FFTPF Fourier;
	Fourier.fft（sn Sk）;
	for （int i = 0; i < length / 2 + 1; i++） {
		amplitude[i] = pow（getValue（Sk[i]） 2） * 2 / （Fs*length）;
		frequency[i] = （double）i * Fs / length;
	}

	Fourier.~FFTPF（）;
	sn.~Vector（）;
	Sk.~Vector（）;
	delete[] park;
}

void ParkTransform::extended_parker_vector（double * a double * b double * c int length int Fs double * amplitude double * frequency）
{
	double* id = new double[length];
	memset（id 0 sizeof（double）*length）;
	double* iq = new double[length];
	memset（iq 0 sizeof（double）*length）;
	for （int i = 0; i < length; i++） {
		id[i] = （sqrt（2） / sqrt（3））*a[i] - （1 / sqrt（6））*b[i] - （1 / sqrt（6））*c[i];
		iq[i] = （1 / sqrt（2））*b[i] - （1 / sqrt（2））*c[i];
	}

	double* park = new double[length];
	memset（park 0 sizeof（double）*length）;
	for （int i = 0; i < length; i++） {
		park[i] = sqrt（id[i] * id[i] + iq[i] + iq[i]）;
	}

	Vector> sn Sk;
	sn.resize（length）;
	Sk.resize（length）;

	for （int i = 0; i < length; ++i）
		sn[i] = park[i];
	/******************************** [ FFT ] ********************************/
	FFTPF Fourier;
	Fourier.fft（sn Sk）;
	for （int i = 0; i < length / 2 + 1; i++） {
		amplitude[i] = pow（getValue（Sk[i]） 2） * 2 / （Fs*length）;
		frequency[i] = （double）i * Fs / length;
	}

	Fourier.~FFTPF（）;
	sn.~Vector（）;
	Sk.~Vector（）;
	delete[] park;
	delete[] id;
	delete[] iq;
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2020-05-27 15:36  ParkTransformTest\
     目录           0  2020-05-27 15:36  ParkTransformTest\.vs\
     目录           0  2020-05-27 15:36  ParkTransformTest\.vs\ParkTransformTest\
     目录           0  2020-05-27 15:36  ParkTransformTest\.vs\ParkTransformTest\v15\
     文件       39424  2020-05-27 15:36  ParkTransformTest\.vs\ParkTransformTest\v15\.suo
     文件     6389760  2020-05-27 15:36  ParkTransformTest\.vs\ParkTransformTest\v15\Browse.VC.db
     目录           0  2020-05-20 11:48  ParkTransformTest\.vs\ParkTransformTest\v15\ipch\
     文件      327680  2020-05-27 15:25  ParkTransformTest\.vs\ParkTransformTest\v15\ipch\ec71f23617bdfdcb.ipch
     目录           0  2020-05-27 15:36  ParkTransformTest\Debug\
     文件      212480  2020-05-27 15:29  ParkTransformTest\Debug\ParkTransformTest.exe
     文件     2040668  2020-05-27 15:29  ParkTransformTest\Debug\ParkTransformTest.ilk
     文件     1683456  2020-05-27 15:29  ParkTransformTest\Debug\ParkTransformTest.pdb
     目录           0  2020-05-27 15:36  ParkTransformTest\ParkTransformTest\
     文件        1464  2020-05-20 11:48  ParkTransformTest\ParkTransformTest.sln
     文件        2206  2010-12-10 21:43  ParkTransformTest\ParkTransformTest\constants.h
     目录           0  2020-05-27 15:36  ParkTransformTest\ParkTransformTest\data\
     文件       24852  2018-11-23 17:12  ParkTransformTest\ParkTransformTest\data\1.txt
     文件       24858  2018-11-23 17:12  ParkTransformTest\ParkTransformTest\data\2.txt
     文件       24860  2018-11-23 17:12  ParkTransformTest\ParkTransformTest\data\3.txt
     文件       23501  2018-11-23 17:12  ParkTransformTest\ParkTransformTest\data\4.txt
     文件       23429  2018-11-23 17:12  ParkTransformTest\ParkTransformTest\data\5.txt
     文件       23850  2018-11-23 17:12  ParkTransformTest\ParkTransformTest\data\6.txt
     目录           0  2020-05-27 15:36  ParkTransformTest\ParkTransformTest\Debug\
     目录           0  2020-05-27 15:36  ParkTransformTest\ParkTransformTest\Debug\ParkTran.DD441EE7.tlog\
     文件        2816  2020-05-27 15:29  ParkTransformTest\ParkTransformTest\Debug\ParkTran.DD441EE7.tlog\CL.command.1.tlog
     文件       28394  2020-05-27 15:29  ParkTransformTest\ParkTransformTest\Debug\ParkTran.DD441EE7.tlog\CL.read.1.tlog
     文件        2420  2020-05-27 15:29  ParkTransformTest\ParkTransformTest\Debug\ParkTran.DD441EE7.tlog\CL.write.1.tlog
     文件        3262  2020-05-27 15:29  ParkTransformTest\ParkTransformTest\Debug\ParkTran.DD441EE7.tlog\link.command.1.tlog
     文件        8340  2020-05-27 15:29  ParkTransformTest\ParkTransformTest\Debug\ParkTran.DD441EE7.tlog\link.read.1.tlog
     文件         816  2020-05-27 15:29  ParkTransformTest\ParkTransformTest\Debug\ParkTran.DD441EE7.tlog\link.write.1.tlog
     文件         208  2020-05-27 15:29  ParkTransformTest\ParkTransformTest\Debug\ParkTran.DD441EE7.tlog\ParkTransformTest.lastbuildstate
............此处省略27个文件信息