0

大家好。所以我正在研究隨機微積分,偶爾我們需要計算一些複雜分佈的積分(從無窮大到無窮大)。在這種情況下,它是Matlab無法計算無限積分?

Dist. being integrated

與正確的答案。這是我投入Matlab的代碼(和我有符號數學工具箱),這Matlab的根本無法處理:

>> syms x t 
>> f = exp(1+2*x)*(1/((2*pi*t)^0.5))*exp(-(x^2)/(2*t)) 
>> int(f,-inf,inf) 

ans = 

-((2^(1/2)*pi^(1/2)*exp(2*t + 1)*limit(erf((2^(1/2)*((x*1i)/t - 2i))/(2*(-1/t)^(1/2))), x, -Inf)*1i)/(2*(-1/t)^(1/2)) - (2^(1/2)*pi^(1/2)*exp(2*t + 1)*limit(erf((2^(1/2)*((x*1i)/t - 2i))/(2*(-1/t)^(1/2))), x, Inf)*1i)/(2*(-1/t)^(1/2)))/(2*pi*t)^(1/2) 

這個答案在最後看起來像nonesense,而鎢(通過其免費工具),即送我是上面圖片的答案。我錯過了在matlab中做這樣的集成的基礎知識,以至於基礎的數學工作頁面沒有涵蓋?請有人幫我理解我的錯誤處理方式嗎?

回答

0

int()第二個參數必須要集成在變量(它看起來像t在這種情況下):

syms x t 
f = exp(1+2*x)*(1/((2*pi*t)^0.5))*exp(-(x^2)/(2*t)) 
int(f,'t',-inf,inf) % <- integrate over t 
+0

這樣做,我得到另一個答案nonesense, ANS = 分段(x in Dom :: ImageSet(x * 1i,x,R_)&in(x^2,'real'),exp(2 * x + 1)*(Inf_limit(t^*(x^2/t)^(1/2)-t ^(1/2)* erfc((x^2 /(2 * t))^(1/2))*(x^2/t )^(1/2)+ .... (輸出太長,無法在此評論中發佈) 這裏發生了什麼問題? – Coolio2654