[運算] 解三角函數的定積分,解出來還是積分?

看板MATLAB作者 (喬弟)時間10年前 (2015/05/11 19:33), 編輯推噓0(0012)
留言12則, 3人參與, 最新討論串1/1
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://www.ptt.cc/bbs/MATLAB/M.1431344033.A.E4A.html

05/12 15:52, , 1F
因為解不出來
05/12 15:52, 1F

05/12 19:55, , 2F
但這是從paper上看到的,應該是可以解!?
05/12 19:55, 2F

05/12 20:33, , 3F
因為有奇點吧
05/12 20:33, 3F

05/12 20:35, , 4F
但是用quad就可以解,不過解出來是個值,其他變數也不見了
05/12 20:35, 4F

05/12 23:19, , 5F
我曾經遇過某些積分matlab積不出來 但是手算是可以積的
05/12 23:19, 5F

05/12 23:30, , 6F
這個手算我試過,但……似乎不太可能耶……
05/12 23:30, 6F

05/13 00:03, , 7F
那就用數值積分吧...
05/13 00:03, 7F

05/13 00:10, , 8F
matlab的ans不就是你的積分式嗎?
05/13 00:10, 8F

05/13 10:44, , 9F
對…可是我想變成其他變數的函數
05/13 10:44, 9F

05/13 10:59, , 10F
這個論文如果沒給closed form,很有可能只是利用數
05/13 10:59, 10F

05/13 10:59, , 11F
值積分去算出答案而已
05/13 10:59, 11F

05/13 11:00, , 12F
畢竟有奇異點,matlab不保證一定能解
05/13 11:00, 12F
文章代碼(AID): #1LK9EXvA (MATLAB)
文章代碼(AID): #1LK9EXvA (MATLAB)