[問題] 極值問題已回收
我是要算
0.5
16∫ x^a (1-x)^(a-1) dx
B(2a-1,2a-1) 0
——————— - ———————————— 在 (0.5,3) 上的最小值
[B(a,a)]^2 B(a,a)
我的程式是：
p=@(a) beta(2*a-1,2*a-1)/beta(a,a)^2;
s=@(x,a) x^a*(1-x)^(a-1);
q=@(a) quad(s,0,0.5)/beta(a,a);
fminbnd(p-16*q,0.5,3)
結果這樣算行不通
我猜是s和q出了問題
quad(s,0,0.5)是把 s=@(x,a) x^a*(1-x)^(a-1) 對x在(0.0.5)上積分
我的意思是 quad(s,0,0.5) 是把 s=@(x,a) x^a*(1-x)^(a-1) 對x在(0.0.5)上積分
因此quad(s,0,0.5)只是a的函數而已
但是這樣表示法又似乎不對
請高手們解說一下
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 58.115.20.204
※ 編輯: buttermilk 來自: 58.115.20.204 (10/23 21:40)
討論串 (同標題文章)
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章