Re: [問題] 請問一題積分的寫法已回收
※ 引述《suntimes (阿肯)》之銘言：
: 想請教的是一題單變數(q)的積分
: g(k) = 積分(k~無窮) (q-k)*f(q) dq
: 其中f(q) = (q*s*根2pi)^(-1) * exp[-(ln(q)-m)^2/2s^2] <--lognormal(m,s)的pdf
: k,m,s都是會給定值的參數，想要計算出g(k)這個積分出來的結果
: 請問一下在MATLAB程式中應該怎麼打呢?
: 小弟是程式新手才剛剛開始接觸，麻煩解釋的詳細一些 > <
: 謝謝!!
m,v分別是輸入的參數
mean = lognstat(m,v);
k = 0.4*mean;%您的預估值
F = @(x) (k*lognpdf(x,m,v));
^^^^^^^^^^^^^^^^^^
Q = quad(F,k,10);
(我有試著把數學方程式寫在另一個function來呼叫，可是會有bug @@)
我的問題在於底線附註的地方，無法加入x
猜想應該是我用"quad"的問題所致!
因為出現的錯誤訊息：
??? Error using ==> mtimes
Inner matrix dimensions must agree.
Error in ==> @(x)(x*lognpdf(x,m,v))
Error in ==> quad at 77
y = f(x, varargin{:});
Error in ==> lognormal at 11
Q = quad(F,k,10);
因為原po要求期望值，即(x-k)*lognpdf(x,m,v)的積分，範圍是k~∞
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.37.49
→
11/21 20:31, , 1F
11/21 20:31, 1F
→
11/21 20:31, , 2F
11/21 20:31, 2F
討論串 (同標題文章)