Doubt in declaring a function

1 次查看（过去 30 天）
Pankaj
Pankaj 2016-1-28
评论： Walter Roberson 2016-1-28
I have a doubt regarding declaration of a function, kindly conside the following code
fun = @GVF
sqrt_estimt = fminsearch(@(theta)ODE_fit(fun,t,y,Q,theta(1),theta(2), iniVal), theta0);
..
function err = ODE_fit(fun, exp_t, exp_y, Q, theta(1), theta(2), iniVal)
% exp_y = Experimental observation at time exp_t
[t,y] = ode45(@(t,X)fun(t,X,n, S0, Q, theta(1),theta(2)),exp_t,iniVal);
err = sum(sum((y-exp_y).^2)); % compute error between experimental y and fitted y
end
..
function dy = GVF(x, y, Q, T, g)
% T = theta(1)
% g = theta(2);
A = y*T;
V = Q/A;
P = T+2*y;
R = A/P;
Sf = (.14*V)^2/(R^(4/3));
Fr = V/sqrt(g*y);
dy = (Sf)/(1-Fr^2);
end
..
The above does not work, but the following does: Is there a way to make the above way function? Thanks
sqrt_estimt = fminsearch(@(theta)ODE_fit(fun, t, y, Q, theta, iniVal), theta0);
..
function err = ODE_fit(fun, exp_t, exp_y, Q, theta, iniVal)
% exp_y = Experimental observation at time exp_t
[t,y] = ode45(@(t,X)fun(t,X,n, S0, Q, theta(1),theta(2)),exp_t,iniVal);
err = sum(sum((y-exp_y).^2)); % compute error between experimental y and fitted y
end
..
function dy = GVF(x, y, Q, T, g)
T = theta(1)
g = theta(2);
A = y*T;
V = Q/A;
P = T+2*y;
R = A/P;
Sf = (.14*V)^2/(R^(4/3));
Fr = V/sqrt(g*y);
dy = (Sf)/(1-Fr^2);
end

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

采纳的回答

Walter Roberson
Walter Roberson 2016-1-28
No, there is not. You cannot name an element of a matrix in a function header. You can use two different variables though.
function err = ODE_fit(fun, exp_t, exp_y, Q, theta1, theta2, iniVal)
% exp_y = Experimental observation at time exp_t
[t,y] = ode45(@(t,X)fun(t,X,n, S0, Q, theta1,theta2),exp_t,iniVal);
err = sum(sum((y-exp_y).^2)); % compute error between experimental y and fitted y
end
  2 个评论
Pankaj
Pankaj 2016-1-28
Thank you Walter, you cleared my doubt.
Walter Roberson
Walter Roberson 2016-1-28
Please Accept the answer to indicate you are finished with the Question.

请先登录，再进行评论。

更多回答（0 个）

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by