涉及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))