在 matlab 中声明函数时使用“int”

Question

我正在尝试定义一个看起来有点像这样的函数：

function u = expectedutility(x)
u = int(utility(x,I)*density(I),I,-10,10)

function t = utility(x,I)
m = x+I
if m>=0 u = m^(1/2)
else u = -2*((-m)^(1/2))

end
end
function s = density(I)
s = 1/20
end

end

在上面的示例中，我使用了上限 10 和下限 -10 的均匀分布，但我希望能够对其进行概括，以便函数密度（I）可以是随机变量 I 的任何 pdf . 另外我知道我可以使用分段（。）来指定实用程序，但我使用的是 2016a，所以它不在我的 Matlab 版本上。

基本上我的问题是每当我点击"expectedutility(2)" 时都会收到错误提示

"Undefined function 'int' for input arguments of type 'double'."

我希望 'int' 集成关于 I 的函数，插入 -10 和 10 说，然后让我观察 u 的不同 x 值。

提前致谢

Answer 1

您的代码中有几处错误，但我认为您尝试做的是集成一个具有符号变量 I 和数字变量 x 的函数。

您的代码中的主要问题是您根本没有定义I。然而，即使你这样做了，你还有一个次要问题：符号函数中不允许使用条件。还有一种替代方法，即使用heaviside 函数。

如果您将utility 函数定义为：

function u=utility(x,I)
m = x+I;
u = m^(1/2)*heaviside(m); % for positives
u =u + -2*( (-m)^(1/2) )*heaviside(-m);    % for negatives  
end

然后你可以调用：

density=1/20;
syms I;
x=2;
u = int(utility(x,I)*density,I,-10,10)

它给出了：

u=
  (4*3^(1/2))/5 - (16*2^(1/2))/15

要获得数字结果，您只需调用double(u)

ans =

  -0.1229

---

如果您不想要单独的功能，只需这样做：

density=1/20;
syms I;
x=2;

m = x+I;
u = m^(1/2)*heaviside(m); % for positives
u =u + -2*( (-m)^(1/2) )*heaviside(-m);    % for negatives 

uint = int(u*density,I,-10,10)

Answer 2

整数的int 函数不采用任何双精度数据类型作为输入。您必须使用syms x（符号标量变量）使其兼容。您可以查看here。代码的另一个问题是变量I 未定义。