PDA

查看完整版本 : [讨论]我的函数为何得不到想要的结果?


linzuhan2009
2009-03-20, 18:12
简单说明一下:
这里,y是已知,x是所求,gy是x的函数,e是y*gy在0到正无穷的积分,故e也是x的函数。
function x=fire(y)
gy='2/(sqrt(2*pi)*x))*exp(-y^2/(2*x^2))';
e='int(y*gy,y,0,+inf)';
x=solve('e=2*x/sqrt(2*pi)','x');
结果如下:
ans =

1/2*e*pi^(1/2)*2^(1/2)

而我需要的答案是没有e的,(这里e是x的函数)。
请大家帮我看看,应该怎样写程序。

slgu
2009-03-21, 07:55
正态分布的积分没有原函数,不能以函数的方式表示曲线下的面积!int对该被积函数(正态概率密度函数)无效。正态分布的积分应用诸如辛普森积分法才行。当然,matlab本身就带有正态分布的累积概率函数(概率)命令:normcdf

linzuhan2009
2009-03-21, 10:28
正态分布的积分没有原函数,不能以函数的方式表示曲线下的面积!int对该被积函数(正态概率密度函数)无效。正态分布的积分应用诸如辛普森积分法才行。当然,matlab本身就带有正态分布的累积概率函数(概率)命令:normcdf

你好,问题还是没有得到解决。