概率密度的窗函数估计－Parzen窗法试验报告-Read.PDF

概率密度的窗函数估计－Parzen 窗法实验报告 一、 算法原理介绍（略） 二、 算法实现步骤： （1） 由计算机的随机函数rand （）产生的0 －RAND_MAX 上的均匀分布随机数产 生（0 ，1）上的均匀分布，再经过变换产生实验中所需的（0 ，1）正态分布和 双峰分布。 （2 ） 建立四种窗函数：方窗、正态窗、指数窗、三角窗，方便后续计算中的窗函数 调用。 （3 ） 由（1）中的函数产生实验所需服从一定分布的随机数据，利用教材（5－4 －12） 式计算在一定x 值下的概率密度估计。 （4 ） 变换x 值，产生一概率密度估计的曲线，并进行绘制。 （5 ） 不断变换参数h1 窗宽度和N 模式的数量以及窗函数，绘制出概率密度估计图。 三、 算法实验结果： （1） 参数设置： 实验中设计一参数设计对话框对参数进行设计，如图： 注：后续实验采用正态窗进行。 （2 ）（0，1）正态分布的P－窗估计结果 (3)双峰分布的P－窗估计结果，分布函数见教材140 页 四、 源程序介绍： （1）介绍算法实现部分的数据定义与函数定义，实现部分略： // Algorithm.h: interface for the Algorithm class. // ////////////////////////////////////////////////////////////////////// //Parzen 窗法进行概率密度估计的算法实现类 /////////////////////////////////////////////////////////////////////// class Algorithm { public: Algorithm(); virtual ~Algorithm(); public: int N;//模式的数量 double *pPat;//模式指针 double Step;//概率密度估计时的x 变量步长 //double x;//概率密度估计时的变量x double X_Start,X_End;//计算概率密度时x 的起点和终点 double Y_Start,Y_End;//绘图时y 轴的范围 double *pP_Est;//概率密度的估计 double h1;//窗宽度 public: static double Abs(double);//绝对值函数 double RandNorm(double mean,double variance);//产生正态随机数 double RandMean(double start,double end);//产生从 start 到 end 的均匀分布随机数 [start,end) static double SquareWinFun(double u);//方窗函数 static double NormalWinFun(double u);//正态窗函数 static double ExpWinFun(double u);//指数窗函数 static double TriangleWinFun(double u);//三角窗函数 void SetPatNum(int n){N=n;}//设置模式的数量 void SetWinWidth(double width){h1=width;}//设置窗的宽度 void GeneratePat(int num,int kind=0);//随机产生num 个模式,kind 为0 正态,1 双峰均 值 double Rand2Mean();//产生双峰均匀分布 void CalculateP(double(*win