涉及3個變量

Question

功能二重積分假設3個變量的函數爲F(T,x,y)=ln(T-x)exp(T+y^3)，我想計算像

H(T)=\int_0^T \int_0^x ln(T-x)exp(T+y^3) dy dx

事情我想的是，對於T的不同值，比如0：0.25： 2，我可以立即獲得H（0：0.25：2）。我使用的代碼如下：

fun=@(T,x,y) ln(T-x)exp(T+y^3) 
ymax=@(x)x 
H=@(T)arrayfun(@(s)integral2(fun,0,s,0,ymax),T) 

但它不起作用。 Matlab給出錯誤。如何解決這個問題呢？

Answer 1

在Matlab中，自然對數不是ln()，而是log()。

首先，改變你的功能如下圖。

fun=@(x,y) log(T-x).*exp(T+y.^3) 

在這裏，我用於產品的操作者*，並且爲了描述矩陣乘積在*和^前加入.。

一個例子是下面，根據T.獲得積分

T=0:0.25:5; 

N=length(T); 

all_H=zeros(1,N); 

s=2; 

for k=1:N; 

    fun=@(x,y) log(T(k)-x).*exp(T(k)+y.^3) 
    ymax=@(x)x 
    H=integral2(fun,0,s,0,ymax) 

    all_H(k)=H; 

end 

plot(abs(all_H))