How to pass an array into the initial condition function while using pdepe

Question

Sachin Hegde 2023-4-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1950553-how-to-pass-an-array-into-the-initial-condition-function-while-using-pdepe

编辑： Torsten 2023-4-20

I have a PDE which is solved using pdepe. However after the solution (1xn array) is achieved, i perform some other calculation and then i want to feed the end result back into the PDE (in the next iteration in a for loop) as initial condition. How can this be achieved.

Additinally how can additional paramters be passed into the pde function?

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

请先登录，再进行评论。

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

Answer 1

Torsten 2023-4-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1950553-how-to-pass-an-array-into-the-initial-condition-function-while-using-pdepe#answer_1219928

编辑：Torsten 2023-4-20

在 MATLAB Online 中打开

However after the solution (1xn array) is achieved, i perform some other calculation and then i want to feed the end result back into the PDE (in the next iteration in a for loop) as initial condition.

Define

icfun = @(x)interp1(X,U,x)

where X is the mesh vector and U is the vector of initial values corresponding to X you want to feed back into the solver for the next iteration.

Additinally how can additional paramters be passed into the pde function?

https://de.mathworks.com/help/optim/ug/passing-extra-parameters.html

4 个评论
显示 2更早的评论隐藏 2更早的评论

Sachin Hegde 2023-4-20

That means naming the mesh vector with ‚x‘ is not a good idea. So It’s better to rename the meshvector with X and then use x for interpolation. Right?

Torsten 2023-4-20

编辑：Torsten 2023-4-20

x is just a formal parameter - MATLAB does not expect it to have a value when you define the function "icfun". It will not conflict with variables explicitly set in your code. Simply make sure that X and U in your definition of "icfun" are your solution mesh vector and your solution from the first call to "pdepe". The name "x" can be replaced by any other name in the definition of "icfun".

Maybe before proceeding you should pass the free MATLAB online tutorial:

https://matlabacademy.mathworks.com/details/matlab-onramp/gettingstarted

请先登录，再进行评论。