How can I compute numerical integration from a multivariate function ( f(t,C) ) where C'=[c_1,c_2,...,c_m] with respect to " t "

Question

Siavash Kardar Tehran 2013-10-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/90447-how-can-i-compute-numerical-integration-from-a-multivariate-function-f-t-c-where-c-c_1-c_2

评论： Yannick 2013-10-17

I want to find the numerical integration from a multivariate function ( f(t,C) ) where C'=[c_1,c_2,...,c_m] with respect to " t ",

what command is the most accurate and efficient Matlab command for doing this?

for example this is my function:

f(t,c_1,c_2,c_3,c_4) =
(t^2 + 1)*((5040315906112595*t^(3/5)*(156*c_1 - 156*3^(1/2)*c_2 + 156*5^(1/2)*c_3 - 156*7^(1/2)*c_4 + 195*3^(1/2)*c_2*t - 585*5^(1/2)*c_3*t + 1170*7^(1/2)*c_4*t + 450*5^(1/2)*c_3*t^2 - 2250*7^(1/2)*c_4*t^2 + 1250*7^(1/2)*c_4*t^3))/702561541869797376 + 1)^2 + (sin(t)*(c_1 - (exp(t)*((5040315906112595*t^(3/5)*(156*c_1 - 156*3^(1/2)*c_2 + 156*5^(1/2)*c_3 - 156*7^(1/2)*c_4 + 195*3^(1/2)*c_2*t - 585*5^(1/2)*c_3*t + 1170*7^(1/2)*c_4*t + 450*5^(1/2)*c_3*t^2 - 2250*7^(1/2)*c_4*t^2 + 1250*7^(1/2)*c_4*t^3))/702561541869797376 + 1))/2 + 3^(1/2)*c_2*(2*t - 1) + 5^(1/2)*c_3*(6*t^2 - 6*t + 1) + 7^(1/2)*c_4*(20*t^3 - 30*t^2 + 12*t - 1))^2)/cos(t)^2

note that the command " int ( f,t,0,1 ) " didn't work and it's said: Warning Explicit integral could not be found.

I'll be very appreciate if you help me with this.

Best regards,

Siavash.

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

请先登录，再进行评论。

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

Answer 1

Yannick 2013-10-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/90447-how-can-i-compute-numerical-integration-from-a-multivariate-function-f-t-c-where-c-c_1-c_2#answer_99927

编辑：Yannick 2013-10-16

Hi - a few things:

INT is a Symbolic Math Toolbox command, so it will try to do the integration symbolically. In this case, since the function is quite complicated, it's having a hard time doing so.
Instead, you probably want to try INTEGRAL if you have a recent enough MATLAB version, or QUAD otherwise.
To handle the extra parameters, check out this doc page: Parameterizing Functions

Hope this helps!

2 个评论
显示无隐藏无

Siavash Kardar Tehran 2013-10-17

Dear Yannick, I've tried the "integral" command but it didn't work out or maybe I couldn't use it correctly, my matlab version is R2013a, so please give me more specific answer, If you can just copy my function to your command window and then try to integrate the function in the way you think it might work then if it worked then tell me the exact commands and how I must use them, Thank you very much, regards, Siavash.

Yannick 2013-10-17

Hi Siavash,

Could you show your code here, as well as how you execute it and what result/warning/error you get? Then I can help figure out with you why it isn't working as expected.

Yannick

请先登录，再进行评论。