clear,clc

%-----------參數設置---------------------------
P=0.3;%0,1兩點分布隨機數，1的概率
m=10000;%兩點分布隨機數的產生個數
n=2000;%目標值x的個數

%--------------計算目標值，畫直方圖--------------
x=zeros(1,n);
for i=1:n
    R=binornd(1,P,1,m);%產生m個，1的概率為P的0，1兩點分布
    for j=1:m
        x(i)=x(i)+R(j)/(2^j);%計算x，m個兩個分布構造一個x
    end
end
hist(x)%直方圖

%-----------parzen窗函數法----------------------
%-------------設置步長，h值------------------------
xx=0:0.01:1;%步長設置為0.01
h=1/sqrt(n);%設置h
%---------------估計概率-------------------------
len=length(xx);
p=zeros(len,1);%存儲xx的估計，最后估計得到點對[xx,p]
for j=1:len
    fy=0;
    for i=1:n
        %選擇正態窗函數
        fy=fy+1/h*1/sqrt(2*pi)*exp(-((xx(j)-x(i))/h)^2/2);
    end
    p(j)=1/n*fy;%xx(j)點的估計概率p(j)
end
figure,
plot(xx,p,'r')