模式识别中Parzen窗法的matlab实现

% clear;clc opengl hardwarebasic % 高斯窗 X=normrnd(0,1,1,1000); x=-5:0.01:5; y=normpdf(x,0,1);%标准正态分布 % N=1,h=0.25,1,4 p1=Parzen(X,0.25,1,x); p2=Parzen(X,1,1,x); p3=Parzen(X,4,1,x); subplot(331) plot(x,p1) hold on plot(x,y,'r') legend('h=0.25,N=1') subplot(332) plot(x,p2) hold on plot(x,y,'r') legend('h=1,N=1') subplot(333) plot(x,p3) hold on plot(x,y,'r') legend('h=4,N=1') hold on % N=10,h=0.25,1,4 p4=Parzen(X,0.25,10,x); p5=Parzen(X,1,10,x); p6=Parzen(X,4,10,x); subplot(334) plot(x,p4) hold on plot(x,y,'r') legend('h=0.25,N=10') subplot(335) plot(x,p5) hold on plot(x,y,'r') legend('h=1,N=10') subplot(336) plot(x,p6) hold on plot(x,y,'r') legend('h=4,N=10') hold on % N=100,h=0.25,1,4 p7=Parzen(X,0.25,100,x); p8=Parzen(X,1,100,x); p9=Parzen(X,4,100,x); subplot(337) plot(x,p7) hold on plot(x,y,'r') legend('h=0.25,N=100') subplot(338) plot(x,p8) hold on plot(x,y,'r') legend('h=1,N=100') subplot(339) plot(x,p9) hold on plot(x,y,'r') legend('h=4,N=100')