arithmetic.h

介绍一些常用的数值计算方法

//用FFT计算离散傅立叶(Fourier)变换
void math_fft(double *pr,double *pi,double *fr,double *fi,int,int,int,int);
//pr，当l=0时，存放n个采样值的实部，返回时存放离散傅立叶变换的模
//当l=1时，存放傅立叶变换的n个实部，返回时存放逆傅立叶变换的模
//pi，当l=0时，存放n个采样值的虚部，返回时存放离散傅立叶变换的幅角
//当l=1时，存放傅立叶变换的n个虚部，返回时存放逆傅立叶变换的幅角
//n，输入的点数
//k，满足n=k个2相乘
//fr，当l=0时，返回傅立叶变换的实部，当l=1时，返回逆傅立叶变换的实部
//fi，当l=1时，返回傅立叶变换的虚部，当l=1时，返回逆傅立叶变换的虚部
//l,当l=0时，表示本函数计算傅立叶变换，当l=1时，表示计算逆傅立叶变换
//il，当il=0，表示不要求计算傅立叶变换或逆傅立叶变换的模和幅角
//当il=1，表示要求计算傅立叶变换或逆傅立叶变换的模和幅角


int math_kalman(double *f,double *q,double *r,double *h,double *y,double *x,
				double *p,double *g,int n,int m,int k);
//f，系统状态转移矩阵；q，模型噪声Wk的协方差阵
//r，观测噪声的协方差矩阵；h,观测矩阵
//y，观测向量序列；x，初值和状态向量估计序列
//p,存放初值P0，返回时存放最后时刻的估计误差协方差阵
//g，返回最后时刻的稳定增益矩阵
//n，动态系统的维数；m，观测矩阵的维数；k，观测序列的长度


//alpha-beida-ganma滤波
void albeigan(double *x,double *y,double t,double a,double b,double c,int n);
//x,n个等间隔点上的量测值；y，返回n个等间隔点上的滤波估值
//t，采样周期；a，b，c，滤波器的参数alpha,beida,ganma