作者JorDy (乔弟)
看板MATLAB
标题[运算] 解三角函数的定积分，解出来还是积分？
时间Mon May 11 19:33:51 2015
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)
再怎麽解都是这个积分
请问各位前辈，是我哪里有打错吗？
还是需要其他的方式求解？
谢谢
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 123.195.4.89
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/MATLAB/M.1431344033.A.E4A.html
1F：→ celestialgod: 因为解不出来 05/12 15:52
2F：→ JorDy: 但这是从paper上看到的，应该是可以解!? 05/12 19:55
3F：→ celestialgod: 因为有奇点吧 05/12 20:33
4F：→ JorDy: 但是用quad就可以解，不过解出来是个值，其他变数也不见了 05/12 20:35
5F：→ rex0707: 我曾经遇过某些积分matlab积不出来 但是手算是可以积的 05/12 23:19
6F：→ JorDy: 这个手算我试过，但……似乎不太可能耶…… 05/12 23:30
7F：→ rex0707: 那就用数值积分吧... 05/13 00:03
8F：→ rex0707: matlab的ans不就是你的积分式吗? 05/13 00:10
9F：→ JorDy: 对…可是我想变成其他变数的函数 05/13 10:44
10F：→ celestialgod: 这个论文如果没给closed form，很有可能只是利用数 05/13 10:59
11F：→ celestialgod: 值积分去算出答案而已 05/13 10:59
12F：→ celestialgod: 毕竟有奇异点，matlab不保证一定能解 05/13 11:00