How to define variables in the main program instead of in the function file

Question

dcydhb dcydhb 2019-8-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/474719-how-to-define-variables-in-the-main-program-instead-of-in-the-function-file

评论： dcydhb dcydhb 2019-8-4

采纳的回答： madhan ravi

在 MATLAB Online 中打开

I want to use a loop in the main program to make

a

the different value, and then get different graphics, but without defining

a

in the function file it can't run, how to solve, thank you!

codes are as this

function file

function f = rigid(t,y)
syms a
f = zeros(2,1);    % a column vector
f(1) = a*y(2) * y(1);
f(2) = -y(1) * y(2);

main program

a=1;
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);
[T,Y] = ode45(@rigid,[0 12],[1 1],options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.')

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

madhan ravi 2019-8-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/474719-how-to-define-variables-in-the-main-program-instead-of-in-the-function-file#answer_386029

在 MATLAB Online 中打开

function f = rigid(t,y,a)
% syms a %% not needed
f = zeros(2,1);    % a column vector
f(1) = a*y(2) * y(1);
f(2) = -y(1) * y(2);
%-------
a=1;
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);
[T,Y] = ode45(@(t,y)rigid(t,y,a),[0 12],[1 1],options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.')

3 个评论
显示 1更早的评论隐藏 1更早的评论

madhan ravi 2019-8-3

编辑：madhan ravi 2019-8-3

在 MATLAB Online 中打开

That was also an example you can change other parameters too , it's called parameterization see the link below:

https://in.mathworks.com/help/matlab/math/parameterizing-functions.html

for a = 1:10
    options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);
    [T,Y] = ode45(@(t,y)rigid(t,y,a),[0 12],[1 1],options);
    plot(T,Y(:,1),'-',T,Y(:,2),'-.')
end