zhengtaigaimiguji.m

parzen窗法

%利用Parzen窗法对总体概率密度函数进行估计
%窗函数选择正态窗函数
%课本140页，利用样本估计正态函数概密
clear all;clc;
%--------------------------1024点----------------------------------------%
xi=randn(1,1024);%产生正态随机数，1024点
x=linspace(-2,2,1024);%在-2到2之间产生1024点的线性空间矢量 
p=Parzen(xi,x,1);%调用parzen函数进行概密函数估计
figure(1)
subplot(2,1,1)
plot(x,p);
title('蓝线为估计的正态概密函数曲线，红线是标准的正态概密函数曲线(N=1024,h1=1)')
xlabel('x');ylabel('p');
hold on
pb=normpdf(x,0,1);%画标准正态函数曲线，对比用
plot(x,pb,'r');

xi=randn(1,1024);
x=linspace(-2,2,1024); 
p=Parzen(xi,x,4);
subplot(2,1,2)
plot(x,p);
title('蓝线为估计的正态概密函数曲线，红线是标准的正态概密函数曲线(N=1024,h1=4)')
xlabel('x');ylabel('p');
hold on
pb=normpdf(x,0,1);
plot(x,pb,'r');
%------------------------------------------------------------------------%


%-------------------------------N=2048-----------------------------------%
figure(2)
xi=randn(1,2048);
x=linspace(-2,2,2048); 
p=Parzen(xi,x,1);
subplot(2,1,1)
plot(x,p);
title('蓝线为估计的正态概密函数曲线，红线是标准的正态概密函数曲线(N=2048,h1=1)')
xlabel('x');ylabel('p');
hold on
pb=normpdf(x,0,1);
plot(x,pb,'r');


xi=randn(1,2048);
x=linspace(-2,2,2048); 
p=Parzen(xi,x,4);
subplot(2,1,2)
plot(x,p);
title('蓝线为估计的正态概密函数曲线，红线是标准的正态概密函数曲线(N=2048,h1=4)')
xlabel('x');ylabel('p');
hold on
pb=normpdf(x,0,1);
plot(x,pb,'r');
%------------------------------------------------------------------------%