How to use substitution on function handle?

Question

I-Chun 2023-11-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2051312-how-to-use-substitution-on-function-handle

编辑： I-Chun 2023-11-23

在 MATLAB Online 中打开

I use an example to simplify my problem. Some steps are as follows:

(1) There is a system of linear equations such as

Use analytical method to solve

and

.

(2) And a differential equation is

Substitute

and

solved from step(1) into the differential equation, then use numerical method to solve

with

.

(3) Another differential equation is

Replace A with

solved from step(2), then use numerical method to solve

with

.

is numerically solved, so I use spline which will give a function handle form. Because my differential equation is actually very complicated, it’s not as simple as just multiplying by 6, so I must use substitution here. However, "subs" can only be used for symbolic or double form, are there other ways to use substitution on function handle? Thank you for any suggestions.

clear
clc
% Step 1
syms z x y;
A = [5 2; 1 1];
B = [12*z; 3*z];
Sol = A\B;
xz = matlabFunction(Sol(1));
yz = matlabFunction(Sol(2));
% Step 2
syms a;
ode_a = @(z, a) xz(z)*yz(z)*a;
initial_a = 1;
[z_a, values_a] = ode45(@(z, a) ode_a(z, a), [0, 1], initial_a);
a_sol = @(z) interp1(z_a, values_a, z, 'spline');  % Here a_sol is a function handle.
% Step 3
syms A
ode_b = @(z, b) 6*A;
ode = subs(ode_b, A, a_sol(z));  % This step mistakes, "subs" can not be used for function handle.
initial_b = 0;
[z_b, b_values] = ode45(@(z, b) ode(z, b), [0, 1], initial_b);

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

请先登录，再进行评论。

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

Answer 1

John D'Errico 2023-11-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2051312-how-to-use-substitution-on-function-handle#answer_1358602

编辑：John D'Errico 2023-11-23

ODE45 CANNOT use symbolic parameters. PERIOD. No matter how you try to trick things, you cannot pass a symbolic parameter (here A) into ODE45. For that matter, neither can you do that with INTERP1.

INTERP1 and ODE45 are purely numerical tools.

If you know the value of that parameter, then yes, you can do what you ask. But then it becomes a purely numerical problem.

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

I-Chun 2023-11-23

编辑：I-Chun 2023-11-23

在 MATLAB Online 中打开

Thank you. But if I don't use "subs", the code can run even though I use spline:

% Step 3
syms A
ode_b = @(z, b) 6*a_sol(z);
%ode = subs(ode_b, A, a_sol);
initial_b = 0;
[z_b, b_values] = ode45(@(z, b) ode_b(z, b), [0, 1], initial_b);

ode45 can indeed solve INTERP1 functions

请先登录，再进行评论。