MATLABWEEK10-11LABSOLUTIONPROBABILITIESANDDISTRIBUTION
1.1functionPMF=TwoDiceSumPMF(X)
ifX>=2&X<=7
PMF=(X-1)/36;
elseifX>7&X<=12
PMF=(13-X)/36;
else
PMF=NaN;
end
end
1.2functionCDF=TwoDiceSumCDF(X)
CDF=0;
forX_n=2:min(X,12)
CDF=CDF+TwoDiceSumPMF(X_n);
end
end
%%Q2
A=rand(100,1)9+1;
B=round(rand(100)99+1);?
C=randn(1,10000);
D=randn(10,1000)sqrt(20)+5;
%%Q3
%(A)
lamda=0.1;
mu=1/lamda;
P_2yr=expcdf(2,mu);
%(B)
P_1yr=expcdf(1,mu);
%%Q4
%(i)
P_2boys=0.520.52;
%(ii)
P_1boy1girl=0.52(1-0.52);用binopdf做
%Q5
x=0:20;
lamda=5;
pdf=poisspdf(x,lamda);
cdf=poisscdf(x,lamda);
subplot(1,2,1),bar(x,pdf,1,''w''),title(''PoissonPMF:lamda=5'')
xlabel(''X''),ylabel(''f(X)'')
subplot(1,2,2),bar(x,cdf,1,''w''),title(''PoissonCDF:lamda=5'')
xlabel(''X''),ylabel(''F(X)'')
%%Q6
x=-100:0.5:100;
mu=[00101010];
sigma=[11011050];
N_pdf=length(mu);
pdf=zeros(N_pdf,length(x));
cdf=zeros(N_pdf,length(x));
fori=1:N_pdf
pdf(i,:)=normpdf(x,mu(i),sigma(i));
cdf(i,:)=normcdf(x,mu(i),sigma(i));
subplot(1,2,1),plot(x,pdf(i,:)),title(''NormalPDF'')
holdon
subplot(1,2,2),plot(x,cdf(i,:)),title(''NormalCDF'')
holdon将所有的pdf/cdf放在一张图上了
end
legend(''mu=0,sigma=1'',''mu=0,sigma=10'',''mu=10,sigma=1'',''mu=10,sigma=10'',''mu=10,sigma=50'')
holdoff
functionPDF=Laplacianpdf(X,LAMDA)
PDF=LAMDA/2exp(-LAMDAabs(X));
end
functionCDF=Laplaciancdf(X,LAMDA)
[r,c]=size(X);
CDF=zeros(r,c);
fori=1:r;
forj=1:c;
ifX(i,j)<0;
CDF(i,j)=exp(LAMDAX(i,j))/2
else
CDF(i,j)=1-exp(-LAMDAX(i,j))/2
end
end
end
end（将X作为矩阵数据进行输入后的处理）
x=-10:0.2:10;
lamda=1;
pdf=Laplacianpdf(x,lamda);
plot(x,pdf),title(''LaplacianPDFforLamda=1'')
|
|
