Doubt in declaring a function
1 回表示 (過去 30 日間)
古いコメントを表示
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
0 件のコメント
採用された回答
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 件のコメント
Walter Roberson
2016 年 1 月 28 日
Please Accept the answer to indicate you are finished with the Question.
その他の回答 (0 件)
参考
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!