function on separate interval

13 次查看（过去 30 天）
Oluwaseun Lijoka
Oluwaseun Lijoka 2016-3-25
评论： Oluwaseun Lijoka 2016-3-26
Hello there, I have a function f(x,t) defined on (0,1). The function is however discontinuous at x=1/2 such that we have two separate values f1(x,t) for x in (0,1/2) and f2(x,t) for x in (1/2, 1). How do I program this function with M function file that takes input x and t?. Arbitrary function can be used to explain. Thanks a lot
  4 个评论
显示 2更早的评论
Oluwaseun Lijoka
Oluwaseun Lijoka 2016-3-26
x is allowed to be a vector. Mostly computed via quadrature.@Robertson

请先登录，再进行评论。

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

采纳的回答

Ced
Ced 2016-3-25
编辑：Ced 2016-3-25
This means that you have defined a function with an output (in your case the variable f), but have not assigned it. This would happen if x does not fall into any of your cases.
Your if conditions are a bit strange.
(x<=1/2) will return 0 if false, 1 if true. So, if e.g. x = 0, then x(x<=1/2) = x(1) = 0, which is considered false instead of true.
This should work better:
function f = exact_soln(x,t,c_1,c_2)
fn=@(x)sin(2*pi*x);
if x <= 0.5 % domain section 1
f= 0.5*(fn(x+c_1*t) + fn(x-c_1*t));
else % domain section 2, 0.5 < x
f= 0.5*(fn(x+c_2*t) + fn(x-c_2*t));
end
end
Having an else statement ensures that no matter x, f will always get a value.
Note that the boundaries x == 0 and x == 1 are actually not checked! If you need your wave to stop there, you could either add more if/else statements, or simply bound x before hand, i.e.
x = max(min(x,1),0);
  1 个评论
Oluwaseun Lijoka
Oluwaseun Lijoka 2016-3-25
Thanks Ced, I actually extended the solution by f=0 elsewhere before but your idea is simple and better. I also like the idea of bounding x as well. Thanks a lot.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Matrix Indexing 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by