关于数值积分，向高手请教

coolxh

问题描述：
(1)  x：M个元素的向量,已知，平均值为0
(2) n：0、1或1.5
(3) 被积函数形式：（x-d）^n*f(x) ，式中，f(x)为x的分布（函数），但无法用已知分布函数拟合
(4) 积分限：[d Inf)
(5) d的取值范围：[0，2*sqrt(var(x))]

目标：
输出积分函数与d 的关系曲线

我的计算步骤如下，但结果与预想的差十万八千里
（1）定义 d
n=1.5;
M=length(x);
NN=100;
d=linspace(0，2*sqrt(var(x)),NN);
（2）逐个计算积分函数
for ii=1:NN
     tt=(x-d(ii))^n;
     ind=find(x>=d(ii));
     ff=tt*length(ind)/M;  %觉得这里计算的是概率而非概率密度，但不知该怎么算
  integrate_ff(ii)=trapz(x(ind),ff(ind)); %这里积分计算似乎有也问题
end
（3）画图
plot(d,integrate_ff)


目前毫无头绪，相当着急:funk: ，请高手指点:handshake

非常感谢！

[ 本帖最后由 sigma665 于 2009-4-22 08:58 编辑 ]

sigma665

（x-d）^n*f(x)
程序里面貌似没有乘f（x）