Is any way in Matlab to perform definite integral with conditions in the integration domain?

Question

Oscar 2014-1-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/113892-is-any-way-in-matlab-to-perform-definite-integral-with-conditions-in-the-integration-domain

评论： Oscar 2014-1-28

采纳的回答： Mike Hosea

Basic example;

f(x)= x^2+3*x for x<100 f(x)=x^3 for x>100; a=10; b=200;

d=integral(f(x), a, b)

Note: this is a basic example, I know for something like the equation above with the int function in symbolic tool box can be resolve easily doing two independent indefinite integrals and then transforms its results in functions and then evaluate those functions (one function in the interval a the other function in the interval b) to obtain the result. However the equation I’m working is very complex, and when I try to find its indefinite integral with the int function, I'm receiving the following warning error " the explicit integral can't been found". For that reason I’m trying to figure out if the integral function could work having some limitation in the integration domain.

I will appreciate whatever help that can point me out to resolve this issue. Thanks

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

请先登录，再进行评论。

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

Answer 1

Mike Hosea 2014-1-27

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/113892-is-any-way-in-matlab-to-perform-definite-integral-with-conditions-in-the-integration-domain#answer_122366

编辑：Mike Hosea 2014-1-27

在 MATLAB Online 中打开

I'm not sure I understand your clarification correctly. Here is my best guess

>> f = @(x)(x.^2+3*x).*(x <= 100) + x.^3.*(x > 100);
>> Q1 = @(c)integral(f,10,c,'Waypoints',100);
>> Q2 = @(c)integral(f,c,200,'Waypoints',100);

Then you can calculate the whole integral using Q = Q1(200) = Q2(10), and you can also calculate the two pieces Q1(c) + Q2(c) = Q for any c. The use of Waypoints is for efficiency. If any of the pieces become nonfinite while the other interval is active and finite, you may need to write a MATLAB function instead of using the logical masking technique above, say fun.m:

function y = fun(x)
y = zeros(size(x),'like',x);
for k = 1:numel(x)
    if x(k) <= 100
        y(k) = x(k)^2 + 3*x(k);
    else
        y(k) = x(k)^3;
    end
end

Then Q1 and Q2 would be written

>> Q1 = @(c)integral(@fun,10,c,'Waypoints',100);
>> Q2 = @(c)integral(@fun,c,200,'Waypoints',100);

It occurs to me that you may mean to integrate the pieces as "basis functions" with finite support, in which case you would define two different functions to integrate in Q1 and Q2, respectively, each defined to be zero outside of their support, either using logical masking or using the same MATLAB function approach with if/else/end.

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

Oscar 2014-1-28

Thanks Mike for your answer. It was very helpful to clarify my issue.

请先登录，再进行评论。

Answer 2

Wayne King 2014-1-27

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/113892-is-any-way-in-matlab-to-perform-definite-integral-with-conditions-in-the-integration-domain#answer_122336

编辑：Wayne King 2014-1-27

在 MATLAB Online 中打开

Using integral()

integral(@(x) x.^2+3*x,10,100)
integral(@(x) x.^3,100,200)

You just integrate your first function from 10 to 100 since it is defined for x<100. x^2+3*x

Integrate the second function from 100 to 200 because it is defined for x>100, x^3

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

Oscar 2014-1-27

编辑：Oscar 2014-1-27

Thanks Wayne,

But I forgot explain something my interest is obtain the result only from the lower integration interval in the first f(x) and only the result from the upper integration integral in the second function. Something like

f(x)= x^2+3*x for x<100 f(x)=x^3 for x>100; a=10; b=200;

d= integral = (x^2+3*x, a, NaN or whatever can be used as dummy load) e=integral = (x^3, NaN,b)

For that reason I was trying to use the int functions instead the integral, because with the int function I can obtain the indefinite integral of each function independently and then use each results as functions. Example

d=int(x^2+3*x,x) e=int(x^3,x) g= matlabFuction(d) i= matlabFuction(e) j=g(10) k=i(200)

but I can't resolve my issue using the int function, because when I tried to obtain the indefinite integral of my equation I receive the following warning " the explicit integral can't been found". I’m assuming this error is due the equation is very complex.

Again I appreciate whatever help that can point me out to resolve this issue. Thanks

请先登录，再进行评论。