[求助]matlab求定积分的问题，高手帮忙啊-经管之家官网！

免费学术公开课,扫码加入

各位大侠，

我已知一个Lognormal分布的概率分布函数，小弟想利用matlab计算出x在0，-0.05。。。-0.45下的概率值是多少？

但是我跑出的结果却是复数形式，这是怎么回事啊？？郁闷，请高手指点一二。

我的程序是：

function y=lognormalprob1(x)
% aim to evaluate the loss prob for lognormal(2P or 3P)分布
%pdf=@(x)(1./(sqrt(2.*pi).*b.*(x-r)).*exp(-0.5.*((log(x-r)-u)./b).^2)); % lognormal分布pdf的形式
%F(x)=1./(sqrt(2.*pi).*quad(e.^(-t.^2)./2),0,(log(x-r)-u)./b) %lognormal分布的CDF形式
％u=input('please input your mean parameter:mu=');
％b=input('please input your variable parameter:sigma=');
％r=input('please input your location parameter:r=');%r=0时，即为2参数分布

%为简单起见，这里哥们直接给出了三个参数的值。

u=-0.6051;

b=0.7017;

r=-0.6396;
%=========按照PDF来的
pdf=@(x)(1./(sqrt(2.*pi).*b.*(x-r)).*exp(-0.5.*((log(x-r)-u)./b).^2));
rt=[0
-0.05
-0.1
-0.15
-0.2
-0.25
-0.3
-0.35
-0.4
-0.45
];
n=length(rt);
for i=1:n
y(i)=quad(pdf,-1,rt(i));
end
y

－－－－－－－－－－运行结果

y =

1.0e+003 *

Columns 1 through 6

1.0529 + 1.4131i 1.0528 + 1.4131i 1.0528 + 1.4131i 1.0527 + 1.4131i 1.0527 + 1.4131i 1.0526 + 1.4131i

Columns 7 through 10

1.0525 + 1.4131i 1.0525 + 1.4131i 1.0524 + 1.4131i 1.0524 + 1.4131i

[em16]

「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP