关于数值积分，向高手请教
问题描述：(1)x：M个元素的向量,已知，平均值为0
(2) n：0、1或1.5
(3) 被积函数形式：（x-d）^n*f(x) ，式中，f(x)为x的分布（函数），但无法用已知分布函数拟合
(4) 积分限：[d Inf)
(5) d的取值范围：
目标：
输出积分函数与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 编辑 ] （x-d）^n*f(x)
程序里面貌似没有乘f（x）
页:
[1]