[运算] 解三角函数的定积分,解出来还是积分?

楼主: JorDy (乔弟)   2015-05-11 19:33:51
http://i.imgur.com/1bjJcuN.jpg
想解这个积分
以下是我打的指令
syms r t theta p lcu lsio2 lsi
cost= cos(theta)
rt=(r+t)
a= (2*rt*cost-(2*t/cost))/lcu
b=2*t/(lsio2*cost)
c= (p-2*rt*cost)/lsi
A=rt*cost
f=A/(a+b+c)
int(f,theta,0,pi/2)
小弟本来默认解出来后会是除了theta之外的变量所组成的函数
但是解出来后,却还是一个定积分
ans =
int((cos(theta)*(r + t))/((p - cos(theta)*(2*r + 2*t))/lsi -
((2*t)/cos(theta) - cos(theta)*(2*r + 2*t))/lcu + (2*t)/(lsio2*cos(theta))),
theta, 0, pi/2)
再怎么解都是这个积分
请问各位前辈,是我哪里有打错吗?
还是需要其他的方式求解?
谢谢
作者: celestialgod (天)   2015-05-12 15:52:00
因为解不出来
楼主: JorDy (乔弟)   2015-05-12 19:55:00
但这是从paper上看到的,应该是可以解!?
作者: celestialgod (天)   2015-05-12 20:33:00
因为有奇点吧
楼主: JorDy (乔弟)   2015-05-12 20:35:00
但是用quad就可以解,不过解出来是个值,其他变量也不见了
作者: rex0707 (我就是我)   2015-05-12 23:19:00
我曾经遇过某些积分matlab积不出来 但是手算是可以积的
楼主: JorDy (乔弟)   2015-05-12 23:30:00
这个手算我试过,但……似乎不太可能耶……
作者: rex0707 (我就是我)   2015-05-13 00:03:00
那就用数值积分吧...matlab的ans不就是你的积分式吗?
楼主: JorDy (乔弟)   2015-05-13 10:44:00
对…可是我想变成其他变量的函数
作者: celestialgod (天)   2015-05-13 10:59:00
这个论文如果没给closed form,很有可能只是利用数值积分去算出答案而已毕竟有奇异点,matlab不保证一定能解

Links booklink

Contact Us: admin [ a t ] ucptt.com