用simpson法求对数正态函数的概率密度函数与坐标轴之间的面积的matlab程序对数正态函数的概率密度函数与坐标轴围成的总面积应该为1,但是我用matlab运行出来的一段区间面积就比1大,我找不
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/04 03:14:27
![用simpson法求对数正态函数的概率密度函数与坐标轴之间的面积的matlab程序对数正态函数的概率密度函数与坐标轴围成的总面积应该为1,但是我用matlab运行出来的一段区间面积就比1大,我找不](/uploads/image/z/3502562-50-2.jpg?t=%E7%94%A8simpson%E6%B3%95%E6%B1%82%E5%AF%B9%E6%95%B0%E6%AD%A3%E6%80%81%E5%87%BD%E6%95%B0%E7%9A%84%E6%A6%82%E7%8E%87%E5%AF%86%E5%BA%A6%E5%87%BD%E6%95%B0%E4%B8%8E%E5%9D%90%E6%A0%87%E8%BD%B4%E4%B9%8B%E9%97%B4%E7%9A%84%E9%9D%A2%E7%A7%AF%E7%9A%84matlab%E7%A8%8B%E5%BA%8F%E5%AF%B9%E6%95%B0%E6%AD%A3%E6%80%81%E5%87%BD%E6%95%B0%E7%9A%84%E6%A6%82%E7%8E%87%E5%AF%86%E5%BA%A6%E5%87%BD%E6%95%B0%E4%B8%8E%E5%9D%90%E6%A0%87%E8%BD%B4%E5%9B%B4%E6%88%90%E7%9A%84%E6%80%BB%E9%9D%A2%E7%A7%AF%E5%BA%94%E8%AF%A5%E4%B8%BA1%2C%E4%BD%86%E6%98%AF%E6%88%91%E7%94%A8matlab%E8%BF%90%E8%A1%8C%E5%87%BA%E6%9D%A5%E7%9A%84%E4%B8%80%E6%AE%B5%E5%8C%BA%E9%97%B4%E9%9D%A2%E7%A7%AF%E5%B0%B1%E6%AF%941%E5%A4%A7%2C%E6%88%91%E6%89%BE%E4%B8%8D)
用simpson法求对数正态函数的概率密度函数与坐标轴之间的面积的matlab程序对数正态函数的概率密度函数与坐标轴围成的总面积应该为1,但是我用matlab运行出来的一段区间面积就比1大,我找不
用simpson法求对数正态函数的概率密度函数与坐标轴之间的面积的matlab程序
对数正态函数的概率密度函数与坐标轴围成的总面积应该为1,但是我用matlab运行出来的一段区间面积就比1大,我找不出原因,望大家帮个忙,快要上交了.
这是我写的matlab程序:函数一是simpson法的定义,函数2是概率密度函数.
function s=MySimp(a,b,k)
h=(b-a)/k;
x=[a:h:b];
n=length(x);
ys=f(x(2));
yf=f(x(n));
ymid=0;
yin=0;
for i=1:n-1
ymid=ymid+4*f((x(i)+x(i+1))/2);
end
for i=2:n-1
yin=yin+2*f(x(i));
end
s=(h/6)*(ys+ymid+yin+yf);
function y=f(x)
y=(1/x*0.891*sqrt(2*pi))*exp((-(log(x)-7.59)^2)/2*0.7938);
a=0;
n=8;
b=2000;
ssimp=MySimp(a,b,n)
用simpson法求对数正态函数的概率密度函数与坐标轴之间的面积的matlab程序对数正态函数的概率密度函数与坐标轴围成的总面积应该为1,但是我用matlab运行出来的一段区间面积就比1大,我找不
simpson法的精度取决于你划分的区间
你把n改成1000,看看是不是更靠近1了
看了一下,你的积分程序并没有看出什么问题来(不过我也很久没用过matlab了,近来都用C或者maple)
但是你给出的函数y=(1/x*0.891*sqrt(2*pi))*exp((-(log(x)-7.59)^2)/2*0.7938)
并不是正态分布的概率密度函数。正态分布密度函数里,不会出现log(x)。
另外,不同的补充语言对log(x)的规定不同,有的不是自然对数,...
全部展开
看了一下,你的积分程序并没有看出什么问题来(不过我也很久没用过matlab了,近来都用C或者maple)
但是你给出的函数y=(1/x*0.891*sqrt(2*pi))*exp((-(log(x)-7.59)^2)/2*0.7938)
并不是正态分布的概率密度函数。正态分布密度函数里,不会出现log(x)。
另外,不同的补充语言对log(x)的规定不同,有的不是自然对数,有的表示常用对数,一定别弄混淆了。
收起