How can I solve this integral equation?

조회 수: 3 (최근 30일)
Roberto
Roberto 2017년 4월 11일
답변: Zeeshan Salam 2019년 12월 1일
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!

답변 (3개)

Roger Stafford
Roger Stafford 2017년 4월 12일
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
Jose Lopez
Jose Lopez 2017년 4월 12일
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 2017년 4월 12일
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.

댓글을 달려면 로그인하십시오.


Roberto
Roberto 2017년 4월 29일
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
Roger Stafford
Roger Stafford 2017년 4월 30일
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.

댓글을 달려면 로그인하십시오.


Zeeshan Salam
Zeeshan Salam 2019년 12월 1일
how i implement the integral function of this equation sin square 6 theta dtheta

카테고리

Help CenterFile Exchange에서 Numerical Integration and Differentiation에 대해 자세히 알아보기

제품

Community Treasure Hunt

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

Start Hunting!

Translated by