Function handle is giving wrong results!

2 次查看（过去 30 天）
Ahmad
Ahmad 2016-5-28
评论： Walter Roberson 2016-5-28
This code is clearly giving wrong results:
clc;clear all
a0 = 2;
a1 = 2;
y = [1 1.4 2.3];
x = [1 2 3];
f = @(x,a0,a1) (a0.*x)/(a1+x);
partialf_a0 = @(x,a1) (x./(x+a1));
partialf_a1 = @(x,a0,a1) (-(a0.*x)/((a1+x).^2));
partialf_a0(a1,x)
partialf_a1(x,a0,a1)
f(x,a0,a1)
Substituting a0 and a1 with their values into the equations, gives correct results for partialf_a0:
clc;clear all
a0 = 2;
a1 = 2;
y = [1 1.4 2.3];
x = [1 2 3];
f = @(x,a0,a1) (a0.*x)/(a1+x);
partialf_a0 = @(x,a1) (x./(2+x));
partialf_a1 = @(x,a0,a1) (-(a0.*x)/((a1+x).^2));
partialf_a0(x)
partialf_a1(x,a0,a1)
  3 个评论
显示 1更早的评论
Ahmad
Ahmad 2016-5-28
编辑：Ahmad 2016-5-28
What I mean is substituting x and a1 and a2 values by hand, gives different results than matlab, for example, x = 1, a1 =2, x/(x+a1)=1/3, please run the code and see what u get!
Walter Roberson
Walter Roberson 2016-5-28
Any code that uses "clear all" is clearly wrong. "clear all" does a lot more than you expect.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（2 个）

the cyclist
the cyclist 2016-5-28
编辑：the cyclist 2016-5-28
Just a guess, but where you did
f = @(x,a0,a1) (a0.*x)/(a1+x);
maybe you really want
f = @(x,a0,a1) (a0.*x)./(a1+x);
[Notice the "./" instead of "/".]
Similarly for your other functions.
If you don't understand why, read this documentation page about array vs. matrix operations.
the cyclist
the cyclist 2016-5-28
Please, as I asked earlier, tell us what you get and what you expect to get. I ran this simplified version of your code
f = @(x,a0,a1) (a0.*x)./(a1+x);
f(1,1,2)
The displayed answer is 0.3333. This is the answer I would expect to get.
I am trying to understand what you think would be different. One possibility is that you are getting confused about variable values that are defined before you define the function, which may be irrelevant because the function arguments will be used instead.

类别

在 Help Center 和 File Exchange 中查找有关 Structures 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by