Passing Functions and m-files to functions

9 views (last 30 days)
Hi Mathworks,
I have an m-file which codes up a system of differential equations;
9D_system
This has various parameters in it, and I have a few parameter schemes for specifying the values to these parameters;
P1.m, P2.m, P3.m
I want to be able to specify which parameter scheme I endow the system with, and then solve this system of differential equations with ode15s.
How do I do this?
Moreover, the parameter schemes
P1.m, P2.m, P3.m
each call functions (some of the parameters are functions of other parameters, and so must be calculated by these called functions). Thus the ode will have a series of nested functions being called!!
Any help on this would be greatly appreciated. :)
Linford
  1 Comment
José-Luis
José-Luis on 12 Feb 2013
Pass them as a string and evaluate inside your function.

Sign in to comment.

Accepted Answer

Sean de Wolski
Sean de Wolski on 12 Feb 2013
@Linford, these are the ways we recommend passing extra parameters:
e.g:
ode45(@(t,y)myfun(t,y,@p1)
And
function x = myfun(t,y,pfun)
stuff = pfun();
x = do_things_with(stuff,t,y)
  1 Comment
Jan
Jan on 12 Feb 2013
Edited: Jan on 12 Feb 2013
Or add the parameters as additional arguments to the anoymous function directly. See also: http://www.mathworks.de/matlabcentral/answers/1971

Sign in to comment.

More Answers (2)

Azzi Abdelmalek
Azzi Abdelmalek on 12 Feb 2013

José-Luis
José-Luis on 12 Feb 2013
myStr = 'P1.m'
function [result] = myFun(someArg, myStr)
[str str] = fileparts(myStr);
run(str); %calling your script
%Do your thing
  1 Comment
Linford Briant
Linford Briant on 12 Feb 2013
Edited: Linford Briant on 12 Feb 2013
Hi Jose-Luis,
Thanks for this; it works a treat for my parameter functions, but it isn't working for my ODE files. For example, the following is the first 4 lines of a file making my ODE file.
function xprime = 9D_system(t,x,myStr,L)
[str str] = fileparts(myStr);
run(str); %calling your script
%Do your thing
When I call this in MATLAB I get the error:
[t,x] = ode15s(@(t,x) 9D_system(t,x,myStr,0.001),[0 2e6],x0);
??? Input argument "myStr" is undefined.
Even though I have previously defined
myStr='Parameter1.m'

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!