How do I use a step function in Boundary value problem ?

2022 11 10
1 个回答
更新时间：2022 11 11
1 次查看（30 天）

0 个投票

Hi,
. I am using sol = bvp4c(odefun,bcfun,solinit). I don't know how to include a step function in my differential equation defined in odefun. Is it possible to do so or odefun can only handle predefined functions? Thankyou for reading. Any suggestion is highly appreciated.

9 个评论
显示 7更早的评论 隐藏 7更早的评论

Torsten
Torsten 2022-11-10
Please include the problem you are trying to solve with equations and boundary conditions.
Gaurav Singh
Gaurav Singh 2022-11-10
Thanks for replying on this thread. Here is the differential equation:
y''''-2*y''+(T+beta)*y=0; where y,T are functions of x; -b<x<b ; beta is a parameter.
T is a step function.
T= -Tc -a<x<a
T= Tc+gamma elsewhere
Tc and gamma are constants.
Bcs: at x=+b, x=-b;
y'''-1.7*y'=0 at x=b, x=-b;
y''-0.3*y=0 at x=b, x=-b
If I arbitrarily chose T as a constant, I get some solution. I have guessed an initial solution of the form sin(p*x). As the equation involves a parameter beta, I have to guess one additional bc to solve for beta. I have just chosen y'=1, not sure if this is right bc though.
Torsten
Torsten 2022-11-10
Your ODE coefficients have a jump at x = +/- a. So without transmission conditions at x = +/-a , you won't be able to fix a solution.
Gaurav Singh
Gaurav Singh 2022-11-10
编辑：Gaurav Singh 2022-11-10
Suppose I make a smoth transition around point 'a', is it possoble to put this T(x) in the differential equation through odefun?
Torsten
Torsten 2022-11-10
If around -a and a, then yes.
Gaurav Singh
Gaurav Singh 2022-11-10
Thanks for confirming. Can you share any example how to do that using odefun, in bvp4c?
Torsten
Torsten 2022-11-10
编辑：Torsten 2022-11-10
What is the problem ? Making the discontinuous function T(x) smooth around x=-a and x=a ?
https://www.johndcook.com/blog/2022/03/29/smoothed-step-function/
But I think the solution to your ODE will greatly depend on how you choose the smoothing function. Thus the solution is kind of arbitrary.
Bjorn Gustavsson
Bjorn Gustavsson 2022-11-11
Fourth derivatives brings back memories (very very vague ones) of the one course in solid mechanics and bending beams and bridges. There we had tables of solutions for different kinds of loads and end-conditions. If this is your problem it might be possible to piece together a solution from such characteristic solutions that fits these conditions.
Gaurav Singh
Gaurav Singh 2022-11-11
Thanks Bjorn. The step function makes life a bit difficult.

请先登录，再进行评论。

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

回答（1 个）

Gaurav Singh
Gaurav Singh 2022-11-11

0 个投票

Thankyou all. I guess, I have figured it out. For now, I am able to use heaviside function in my problem. Here is the code;
function dydx = mat4ode(x,y,lambda) % equation being solved
global k
dydx = [y(2)
y(3)
y(4)
2*k^2*y(3)-(1*(-heaviside(x+1)+heaviside(x-1)+0.5)+lambda)*y(1)];
end

2 个评论
显示 无 隐藏 无

Torsten
Torsten 2022-11-11
编辑：Torsten 2022-11-11
This will introduce jumps at x+1 and x-1 for every point x of your grid vector. I doubt that this is what you want.
Gaurav Singh
Gaurav Singh 2022-11-11
编辑：Gaurav Singh 2022-11-11
Thanks for your comment Torsten. Can you please elaborate more on this. May be I am missing something important here. From my understanding, as its a forth order DE, the solver gives me continuous function up to order three. The third derivative is having sharp corner at x=+1 and -1.
I belive fourth derivative will be discontinuous. Is there any fundamental err in suppyling heaviside function directly in DE in Matlab?

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

产品

版本

R2022a

标签

Community Treasure Hunt

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

Start Hunting!

Translated by