jysfgaimiguji.m

parzen窗法

%课本140页，用p窗法对均匀分布双峰概密函数进行估计
clear all;clc;
%----------------------1024点---------------------------------%
xi=randjs(1024);
x=linspace(-3,3,1024);%在-3到3之间产生1024点的线性空间矢量 
p=Parzen(xi,x,1);%调用parzen函数进行概密函数估计
figure(1)
subplot(2,1,1)
axis([-3,3,0,1.5]);
plot(x,p);
title('蓝线为估计的均匀分布双峰概密函数曲线，红线是标准的均匀分布双峰概密函数曲线(N=1024,h1=1)')
xlabel('x');ylabel('p');
hold on
pb=junyunshuangfengsd(x);%画标准均匀分布双峰概密函数曲线，对比用
plot(x,pb,'r');

p=Parzen(xi,x,4);%调用parzen函数进行概密函数估计
subplot(2,1,2)
axis([-3,3,0,1.5]);
plot(x,p);
title('蓝线为估计的均匀分布双峰概密函数曲线，红线是标准的均匀分布双峰概密函数曲线(N=1024,h1=4)')
xlabel('x');ylabel('p');
hold on
pb=junyunshuangfengsd(x);%画标准均匀分布双峰概密函数曲线，对比用
plot(x,pb,'r');
%-------------------------------------------------------------%

%-----------------------2048点---------------------------------%
xi=randjs(2048);
x=linspace(-3,3,2048);%在-3到3之间产生2048点的线性空间矢量 
p=Parzen(xi,x,1);%调用parzen函数进行概密函数估计
figure(2)
subplot(2,1,1)
axis([-3,3,0,1.5]);
plot(x,p);
title('蓝线为估计的均匀分布双峰概密函数曲线，红线是标准的均匀分布双峰概密函数曲线(N=2048,h1=1)')
xlabel('x');ylabel('p');
hold on
pb=junyunshuangfengsd(x);%画标准均匀分布双峰概密函数曲线，对比用
plot(x,pb,'r');

p=Parzen(xi,x,4);%调用parzen函数进行概密函数估计
subplot(2,1,2)
axis([-3,3,0,1.5]);
plot(x,p);
title('蓝线为估计的均匀分布双峰概密函数曲线，红线是标准的均匀分布双峰概密函数曲线(N=2048,h1=4)')
xlabel('x');ylabel('p');
hold on
pb=junyunshuangfengsd(x);%画标准均匀分布双峰概密函数曲线，对比用
plot(x,pb,'r');
%-------------------------------------------------------------%