How can I solve this integral equation?

3 views (last 30 days)
Roberto
Roberto on 11 Apr 2017
Answered: Zeeshan Salam on 1 Dec 2019
Hello everyone, my name is Jose from Sevilla.
I have a function f(T) given by:
f(T)=a+bT+cT^2, where a,b and c are known numbers and T denotes Temperature.
Now I have this equation, where d is yet another constant:
I need to solve for T2, since T1 is also known. In fact, the only unknown here is T2.
I thought about using the trapz function, but I don't know how to include the T2 unknown. Any help will be greatly appreciated! (Important: I don't have the Symbolic Math Toolbox, so I can't do it symbolically. I don't have access to the Optimization Toolbox either, so fsolve and solve are ruled out, too).
Thank you!

Answers (3)

Roger Stafford
Roger Stafford on 12 Apr 2017
I would suggest utilizing a little calculus here:
integral of a+b*T+c*T^2 w.r. T from T = T1 to T = T2
is equal to:
I = a*(T2-T1)+b/2*(T2^2-T1^2)+c/3*(T2^3-T1^3)
or
c/3*T2^3+b/2*T2^2+a*T2-c/3*T1^3-b/2*T1^2-a*T1-I = 0
You have said everything is known except T2, so you can express T2 as the real solution (or solutions) to:
T2 = roots([c/3,b/2,a,-c/3*T1^3-b/2*T1^2-a*T1-I]);
  2 Comments
Jose Lopez
Jose Lopez on 12 Apr 2017
Thank you, since f(T) is always a polynomial your solution will work in my case withouth integrating f.
Never heard of root. Thank you again!
Roger Stafford
Roger Stafford on 12 Apr 2017
You need to be prepared for multiple roots from 'roots'. Presumably just one of them will be the one you want. An n-th order polynomial always yields n roots even though some of them may be complex-valued. That is not a difficulty introduced by 'roots' or matlab. It is inherent in the statement of your problem.

Sign in to comment.


Roberto
Roberto on 29 Apr 2017
One more question. I think the answer wouldn't be as straightforward in case that f(T) was of the form:
f(T) = a + b*T + c*T^2 + d/T
What would be the best approach in this case? Now the roots function wouldn't work since f(T) isn't a polynomial..
  1 Comment
Roger Stafford
Roger Stafford on 30 Apr 2017
I disagree. ‘roots’ will still work here. The equation
f(T) = a + b*T + c*T^2 + d/T
is equivalent to the equation
c*T^3+b*T^2+(a-f(T))*T+d = 0
which can also be solved with ‘roots’ assuming f(T) is known. You will, however, get three solutions to this third degree equation.

Sign in to comment.


Zeeshan Salam
Zeeshan Salam on 1 Dec 2019
how i implement the integral function of this equation sin square 6 theta dtheta

Products

Community Treasure Hunt

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

Start Hunting!