mainfun.asv

这是一个有关parzen窗估计的代码

clc;
clear;
x=randn(1,2000)
px=normpdf(x,0,1)
%Parzen窗h1=0.25时，N不同的估计
subplot(2,2,1);
plot(x,px,'.');
title('原始一维正态分布')

pNx=parzen(1,0.25,x);
subplot(2,2,2);
plot(x,pNx,'.');
title('Parzen窗法估计单一正态分布')
xlabel('h1=0.25,N=1')

pNx=parzen(512,0.25,x);
subplot(2,2,3);
plot(x,pNx,'.');
title('Parzen窗法估计单一正态分布')
xlabel('h1=0.25,N=512')

pNx=parzen(256,0.25,x);
subplot(2,2,4);
plot(x,pNx,'.');
title('Parzen窗法估计单一正态分布')
xlabel('h1=0.25,N=256')
%Parzen窗h1=1时，N不同的估计
figure(2);

subplot(2,2,1);
plot(x,px,'.');
title('原始一维正态分布')

pNx=parzen(1,1,x);
subplot(2,2,2);
plot(x,pNx,'.');
title('Parzen窗法估计单一正态分布')
xlabel('h1=1,N=1')

pNx=parzen(512,1,x);
subplot(2,2,3);
plot(x,pNx,'.');
title('Parzen窗法估计单一正态分布')
xlabel('h1=1,N=512')

pNx=parzen(256,1,x);
subplot(2,2,4);
plot(x,pNx,'.');
title('Parzen窗法估计单一正态分布')
xlabel('h1=1,N=256')
%Parzen窗h1=4时，N不同的估计
figure(3);

subplot(2,2,1);
plot(x,px,'.');
title('原始一维正态分布')

pNx=parzen(1,4,x);
subplot(2,2,2);
plot(x,pNx,'.');
title('Parzen窗法估计单一正态分布')
xlabel('h1=4,N=1')

pNx=parzen(512,4,x);
subplot(2,2,3);
plot(x,pNx,'.');
title('Parzen窗法估计单一正态分布')
xlabel('h1=4,N=16')

pNx=parzen(256,4,x);
subplot(2,2,4);
plot(x,pNx,'.');
title('Parzen窗法估计单一正态分布')
xlabel('h1=4,N=256')