test.asv

一维序列概率密度的估计

clear all
x=linspace(-3,3,100);
x1=randn(1,16);%一维正态分布
x2=randn(1,1024);
p1=Parzen(x1,x,1,[]);
p2=Parzen(x2,x,1,[]);
figure(1)
subplot(1,2,1); plot(x,p1);ylim([0,1]);
title('一维正态分布概率密度估计(N=16)')
subplot(1,2,2); plot(x,p2);ylim([0,1]);
title('一维正态分布概率密度估计(N=1024)')


clear all
x1=rand(1,16);%一维均匀分布序列
x2=rand(1,1024);
x=linspace(-2,2,100);
p1=Parzen(x1,x,1,[]);
p2=Parzen(x2,x,1,[]);
figure(2)
subplot(1,2,1); plot(x,p1);ylim([0,1]);
title('一维均匀分布概率密度估计(N=16)')
subplot(1,2,2); plot(x,p2);ylim([0,1.5]);
title('一维均匀分布概率密度估计(N=1024)')


clear all
k=3;
n=1;m1=16;m2=1024;
B1=chi2rnd(k,n,m1);   %x^2分布
B2=chi2rnd(k,n,m2);
x=linspace(0,10,100);
p1=Parzen(B1,x,1,[]);
p2=Parzen(B2,x,1,[]);
figure(3)
subplot(1,2,1); plot(x,p1);
title('一维x^2分布分布概率密度估计(N=16)')
subplot(1,2,2); plot(x,p2);
title('一维x^2分布分布概率密度估计(N=1024)')

clear all
v=10;
n=1;m1=16;m2=1024;
C1=trnd(v,n,m1); %T分布
C2=trnd(v,n,m2);
x=linspace(-2,2,1024);
p1=Parzen(C1,x,1,[]);
p2=Parzen(C2,x,1,[]);
figure,plot(x,p); 
title('一维T分布分布概率密度估计');
clear all
v1=6;v2=6;
n=1;m=1024;
D=frnd(6,6,1,1024);  %F分布
x=linspace(0,8,100);
p=Parzen(D,x,1,[]);
figure,plot(x,p); 
title('一维F分布分布概率密度估计')