test.m

一维序列概率密度的估计

clc
disp('Enter ...')
disp('0 (default) ——一维正态分布概率密度估计 ')
disp('1 ——一维其他分布概率密度估计')
chos=input('Enter your choice: ');
if isempty(chos), chos=0; end % default choice
if chos==0,

    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,0.7]);
    title('一维正态分布概率密度估计(N=16)')
    subplot(1,2,2); plot(x,p2);ylim([0,0.7]);
    title('一维正态分布概率密度估计(N=1024)')

else
    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(-3,3,1024);
    p1=Parzen(C1,x,1,[]);
    p2=Parzen(C2,x,1,[]);
    figure(4)
    subplot(1,2,1); plot(x,p1);
    title('一维T分布概率密度估计(N=16)')
    subplot(1,2,2); plot(x,p2);
    title('一维T分布概率密度估计(N=1024)')

    clear all
    v1=6;v2=6;
    n=1;m1=16;m2=1024;
    D1=frnd(6,6,n,m1);  %F分布
    D2=frnd(6,6,n,m2);
    x=linspace(0,8,100);
    p1=Parzen(D1,x,1,[]);
    p2=Parzen(D2,x,1,[]);
    figure(5)
    subplot(1,2,1); plot(x,p1);
    title('一维F分布概率密度估计(N=16)')
    subplot(1,2,2); plot(x,p2);
    title('一维F分布概率密度估计(N=1024)')
end