Matlab实现正态分布

1）使用MatLab画出正态分布的概率密度函数图像。
x=[-10:0.01:10];
y=normpdf(x,0,1);%正态分布函数。
figure;
axes1=axes('Pos',[0.1 0.1 0.85 0.85]);
plot(x,y);
set(axes1,'YLim',[-0.01 0.43],'XLim',[-3 3]);
图1：

2）验证概率密度函数在区间(-∞,∞)上的积分为1。
这里取参数mu=3，sigma=5（注：下文全用这两个参数）。
y='exp(-1/2*((x-3)/5)^2)/(sqrt(2*pi)*5)';
s=int(y,-inf,inf) %int积分函数（inf代表无穷大）。
输出：s=1

3）验证x=mu时取最大值。
思路：求解函数一阶导数为零的点。

* 求一阶导数
y='exp(-1/2*((x-3)/5)^2)/(sqrt(2*pi)*5)';
d=diff(y);%微分函数。
sd=simplify(d)
输出结果：sd = -1/250*(x-3)*exp(-1/50*(x-3)^2)*2^(1/2)/pi^(1/2)

* 通过图像判断解的位置
x=[0:0.001:40];
sd=-1/250.*(x-3).*exp(-1/50.*(x-3).^2).*2.^(1/2)./pi.^(1/2);
axes1=axes('Pos',[0.1 0.1 0.85 0.85]);
plot(x,sd);
set(axes1,'YLim',[-0.01 0.01],'XLim',[0 40]);
图2：

从图中可以看出在3附近有解。

* 定义函数并求解
function y=f(x)
y=-1/250.*(x-3).*exp(-1/50.*(x-3).^2).*2.^(1/2)./pi.^(1/2);

r=fzero('f',3)
输出：r=3

从图上看在x > 20以后，几乎是一条直线，若用20：
r=fzero('f',20)
输出：r=3

这说明是无限趋近于0。

* 进一步说明该点为最大值点
该概率密度函数一阶导数为0的解为3，此值正好为mu，再取x=1，x=4与x=3时的函数值比较。
>> normpdf(1,3,5)
ans = 0.0737

>> normpdf(4,3,5)
ans = 0.0782

>> normpdf(3,3,5)
ans = 0.0798

显然在x=3的两边函数值都比x=3小，说明该点为极大值点。根据正态分布函数的图像特点可知该点是最大值点。

4）验证x=mu +- sigma （即8或-2）处曲线有拐点。
思路：求二阶导数为零的点。

* 先求二阶微分
y='exp(-1/2*((x-3)/5)^2)/(sqrt(2*pi)*5)';
d=diff(y,2);%微分函数。
sd=simplify(d)
输出：sd = 1/6250*exp(-1/50*(x-3)^2)*2^(1/2)*(-16+x^2-6*x)/pi^(1/2)

* 通过图像判断解的位置
x=[-20:0.001:20];
sd=1./6250.*exp(-1/50.*(x-3).^2).*2.^(1/2).*(-16+x.^2-6.*x)/pi.^(1/2);
axes1=axes('Pos',[0.1 0.1 0.85 0.85]);
plot(x,sd);
set(axes1,'YLim',[-0.005 0.005],'XLim',[-20 20]);
图3：
从上图可以看出，曲线在(-5,0)和(5,10)之间分别都与y=0有交点，因此有两个解。

* 定义函数并求解
function y=f(x)
y=1./6250.*exp(-1/50.*(x-3).^2).*2.^(1/2).*(-16+x.^2-6.*x)/pi.^(1/2);

r=fzero('f',-5)
r = -2

>> r=fzero('f',5)
r = 8.0000

从而得到了两个拐点x=8和x=-2，也即mu +- sigma。

5）验证曲线以x轴为渐近线渐近线求解：
A 垂直渐近线  x=a是y=f(x)的渐近线<==>lim f(x)=∞或lim f(x)=∞
                         x->a+0       x->a-0
  其中a在间断点中找——∞型第二类间断点B 水平渐近线
  x→+∞（-∞）时，y=b是y=f(x)的渐近线<==>lim f(x)=b （或lim f(x)=b)
                                          x->+∞         x->-∞

求其一阶倒数在x趋向于无穷大时的极限值b，若存在，即有水平渐近线y=b。
* 先定义函数：
function y=f(x)
syms x; %定义符号变量。
y=exp(-1/2*((x-3)/5)^2)/(sqrt(2*pi)*5);

* 求x趋向无穷大时一阶导数的极限
limit(f,inf)
ans = 0

6）验证 3 sigma 法则
思路：求解概率密度函数在[mu-3*sigma,mu+3*sigma]区间上的积分。
y='exp(-1/2*((x-3)/5)^2)/(sqrt(2*pi)*5)';
double(int(y,-12,18))
ans = 0.9973