作者a5170040 (Piggy)
看板NCTU-STAT97G
标题[心得] 函数指标
时间Thu Mar 26 12:43:55 2009
哈~统算写到这个部分来分享一下心得
因为pdf有很多
总不可能每个pdf都写一个积分吧?
(也就是底下程式码Integral函数的内容复制好几次)
这样程式码会变得非常长
所以我找了一下书并且做了实验发现函数指标可以解决这个问题
所以来和大家分享一下
double p1,p2//外部变数(也就是以下的程式都可以用)
double (*pdf)(double); //宣告函数指标 我现在只有要放x而已，所以引数只有一个
//定义gamma pdf
double pdf_gamma(double x){
return (pow(x,p1-1.0)/g1)*pow(p2,-p1)*exp(-x/p2);//我的g1已经先积好了
}
//卡方pdf
double pdf_chi_square(double x){
return exp(-0.5*x)*pow(x,0.5*p1-1.0)/(pow(2,0.5*p1)*g1);
}
int main(void){
...
pdf = pdf_gamma;//把要的pdf给函数指标
//或者可以放
//pdf = pdf_chi_square;
Integral(); 呼叫积分函数
...
}
void Integral(void){
for(i=0;i<分割数;i++){
....
pdf(x1)+pdf(x2);//这边就可以用函数指标作运算了
....
}
}
大致上这样
哈
--
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.113.191.166
1F：推 hiloin:推一个，受益良多。 03/26 15:21
2F：推 liangxun:感谢感谢~~0rz 03/26 20:17
※ 编辑: a5170040 来自: 140.113.191.166 (03/26 20:21)
3F：→ teddd:有看有推... 03/27 19:08
4F：推 samlee04066:推 03/27 22:05
5F：推 yeya:噢我看不懂 03/27 22:27