Parzen窗matlab编程报告.doc

Parzen窗matlab仿真实验一、Parzen窗基本原理Parzen窗估计法是一种具有坚实理论基础和优秀性能的非参数函数估计方法，它能够较好地描述多维数据的分布状态。其基本思想就是利用一定范围内各点密度的平均值对总体密度函数进行估计。一般而言，设x为d维空间中任意一点，N是所选择的样本总数，为了对x处的分布概率密度进行估计，以x为中心作一个边长为的超立方体,则其体积为，为计算落入中的样本数，构造一个函数使得(1)并使满足条件，且，则落入体积V中的样本数为，则此处概率密度的估计值是： (2)式(2)是Parzen窗估计法的基本公式，称为窗函数，或核函数。在Parzen窗估计法的基本公式中，窗宽是一个非常重要的参数。当样本数N有限时，对估计的效果有着较大的影响。二、Matlab实现Parzen窗仿真实验为了编写Parzen窗，首先要确定一个窗函数。一般可以选择的窗函数有方窗、正态窗等。接下来将分别用方窗、正态窗实现Parzen窗法。2.1、方窗 （3）其中，为超立方体的棱长。因为编程中所用的信号是一维的正态信号，所以此处的取值为1。2.1.1 方窗Parzen仿真实验流程图 生成均值为0，方差为1，长度为N的一维正态随机信号用判断是否在以为中心为长度的对称区间内则落入为中心的区间内的样本个数加一YN用（2）式求出绘制估计概率密度函数曲线2.1.2 方窗Parzen仿真实验mtalab代码r=normrnd(0,1,1,10000);f=r(1:1000);f=sort(f);lth=1;h=0.25,1,4;V1=h(1);V2=h(2);V3=h(3);V=V1,V2,V3;a=zeros(1000,3);p=zeros(1000,3);b=0;for m=1:3for x=1:1000 for n=1:lth if abs(r(:,n)-f(x)/h(:,m)=1/2 q=1; else q=0; end b=q+b; enda(x,m)=b;b=0;endendfor m=1:3for x=1:1000p(x,m)=1/(lth*V(m)* a(x,m);endsubplot(4,3,m)plot(f(1:1000),p(:,m)ylabel(N=1)endhold on lth=16;b=0;for m=1:3for x=1:1000 for n=1:lth if abs(r(:,n)-f(x)/h(:,m)=1/2 q=1; else q=0; end b=q+b; end a(x,m)=b;b=0;endendfor m=1:3for x=1:1000p(x,m)=1/(lth*V(m)* a(x,m);endsubplot(4,3,m+3)plot(f(1:1000),p(:,m)ylabel(N=16)endhold on lth=256;b=0;for m=1:3for x=1:1000 for n=1:lth if abs(r(:,n)-f(x)/h(:,m)=1/2 q=1; else q=0; end b=q+b; end a(x,m)=b;b=0;endendfor m=1:3for x=1:1000p(x,m)=1/(lth*V(m)* a(x,m);endsubplot(4,3,m+6)plot(f(1:1000),p(:,m)axis(-5,5,0,0.6);ylabel(N=256)endlth=10000;b=0;for m=1:3for x=1:1000 for n=1:lth if abs(r(:,n)-f(x)/h(:,m)=1/2 q=1; else q=0; end b=q+b; end a(x,m)=b;b=0;endendfor m=1:3for x=1:1000p(x,m)=1/(lth*V(m)* a(x,m);endsubplot(4,3,m+9)plot(f(1:1000),p(:,m)axis(-5,5,0,0.5);ylabel(N=10000)end2.1.3 方窗Parzen仿真实验运行结果 h=0.25 h=1 h=42.2正态窗正态窗函数为 （4）概率密度的估计式为（5） 式中 。2.2.1 正态Parzen窗仿真实验流程图生成均值为0，方差为1，长度为N的一维正态随机信号根据h和求得h1并计算出通过核函数的值用（5）式求出绘制估计概率密度函数曲线2.2.2 正态Parzen窗仿真实验mtalab代码r=normrnd(0,1,1,10000);f=r(1:1000);f=sort(f);lth=1;h=0.25,1,4;a=zeros(1000,3);p=zeros(1000,3);b=0;for m=1:3 h1=h(m)/sqrt(lth);for x=1:1000 for n=1:lth b=b+exp(r(:,n)-f(x)/h1).2/(-2)/sqrt(2*pi)/h1; end a(x,m)=b/lth; b=0;endendfor m=1:3 for x=1:1000p(x,m)=a(x,m); endsubplot(4,3,m)plot(f(1:1000),p(:,m)ylabel(N=1)endhold on lth=16;h=0.25,1,4;a=zeros(1000,3);p=zeros(1000,3);b=0;for m=1:3 h1=h(m)/sqrt(lth);for x=1:1000 for n=1:lth b=b+exp(r(:,n)-f(x)/h1).2/-2)/sqrt(2*pi)/h1; end a(x,m)=b/lth;b=0;endendfor m=1:3 for x=1:1000p(x,m)=a(x,m); endsubplot(4,3,m+3)plot(f(1:1000),p(:,m)ylabel(N=16)endhold on lth=256;h=0.25,1,4;a=zeros(1000,3);p=zeros(1000,3);b=0;for m=1:3 h1=h(m)/sqrt(lth);for x=1:1000 for n=1:lth b=b+exp(r(:,n)-f(x)/h1).2/-2)/sqrt(2*pi)/h1; end a(x,m)=b/lth; b=0;endendfor m=1:3 for x=1:1000p(x,m)=a(x,m); endsubplot(4,3,m+6)plot(f(1:1000),p(:,m)axis(-5,5,0,0.6);ylabel(N=256)end lth=10000;h=0.25,1,4;a=zeros(1000,3);p=zeros(1000,3);b=0;for m=1:3 h1=h(m)/sqrt(lth);for x=1:1000 for n=1:lth b=b+exp(r(:,n)-f(x)/h1).2/-2)/sqrt(2*pi)/h1; end a(x,m)=b/lth; b=0;endendfor m=1:3 for x=1:1000p(x,m)=a(x,m); endsubplot(4,3,m+9)plot(f(1:1000),p(:,m)axis(-5,5,0,0.5);ylabel(N=10000)end2.2.3 正态Parzen窗仿真实验运行结果 h=0.25 h=1 h=42.3 结论从方窗和正态窗两个Parzen窗仿真实验可以看出，估计的概率密度