MATLAB 板


LINE

※ 引述《va1kyrie (123)》之銘言: : 各位先進大家好 : 小弟想請問一個問題 : 我現在要解一個積分: : ∫ (1/(exp((x-F-V)/T))+1) * (R)/((x-E-U)^2+R^2) dx : 我要對x積分 : 積分區間是從-∞到F : 而F=0.2 E=0.3 U=0.06 R=0.001 T=0.002都是已知定數 : V則是一個變數 V=0:0.001:1 後面還會用到 : 請問這樣能積嗎? : 請知道的人教我 : 或是提點我一下該用哪些指令~ 先把方程式寫容易看一些 V 1 * R ∫--------------------------------------------------------dx -∞ x - F - V 2 2 ( exp( ----------- ) + 1 ) * ( ( x - E - U ) + R ) T 邊界是無限大沒辦法直接做數值積分 所以先進行變數變換 1 1 - t -1 令t = --------- => x = ------- => dx = -------dt x + 1 t t * t 1 下界-∞ => ---------- = 0 -∞ + 1 1 上界 V => ---------- V + 1 接著把x替換掉 1 ------- V + 1 - R ∫ ---------------------------------------------------------------dt 0 1 - t / ------- -F -V \ | t | 1 - t 2 2 2 ( exp| -------------- | + 1 ) * ( ( ------- - E - U ) + R ) * t \ T / t 變成有限區間就可以直接做數值積分了 (那直接求-∞開始區間的H開頭高斯積分法我不太信任) (而且H開頭高斯積分法的高斯點與weighting也太難求) 接下來有兩種作法 1 一種是把0到-------的區間再做變數變換處理成-1到1之間 V + 1 然後就可以去找網路上面很多的-1到1高斯積分點資料 直接使用高斯積分法便可得到積分值 因為經過變數變換後知道變數會轉變成怎樣的方程式 所以不管V怎麼變,其區間都會維持在-1到1之間 因此V不管怎麼變都可以使用完全相同的一套最基本的高斯積分點來計算 第二種作法就是不管變數變換法處理成-1到1區間 1 改成求0到-------這個區間內的高斯積分點再來做高斯積分法 V + 1 因為V會變化,所以每次變化都要再求一次新區間內的高斯積分點 不過這個Matlab Center裡面有別人已經寫好的funtion了 只要輸入上下區間以及需要幾個高斯點就會把高斯點和其weighting算出來給你 撿別人寫好的來用會簡單一些,可以自己少算一次變數變換 最後要注意的一點就是高斯點數目請請使用偶數,以免出現0的高斯點來 至於要取幾個高斯點才會轉就麻煩請進行收斂測試來決定 -- Deserves death! I daresay he does. Many that live deserve death. And some die that deserve life. Can you give that to them? Then be not too eager to deal out death in the name of justice, fearing for your own safty. Even the wise cannot see all ends. Gandalf to Frodo --



※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 125.233.241.253 ※ 編輯: Gwaewluin 來自: 125.233.241.253 (10/17 22:48)
1F:推 gamer:沒記錯的話,Gauss point好像沒有0 10/17 23:09
2F:→ gamer:不過這個題目出現0應該也不會有什麼問題 10/17 23:09
3F:→ gamer:喔,轉換後的t代0會爆掉,我看成轉換前的 10/17 23:10
4F:推 Gwaewluin:奇數個又是在-1到1區間的話,中間的高斯點會是0,因為在 10/17 23:11
5F:→ Gwaewluin:-1到1區間高斯點是左右對稱的,所以中間那點會為0 10/17 23:12
6F:推 va1kyrie:我是原問者~ 謝謝你的解答 10/18 00:26
7F:→ va1kyrie:不過積分上限是定值(F=0.2) 不是變數 10/18 00:30
8F:推 Gwaewluin:這不重要,算出變數變換後區間內高斯點即可簡單求出數值 10/18 00:41
9F:→ Gwaewluin:積分,不需要想太複雜,高斯積分很好用低 XD 10/18 00:41
10F:推 zupo:推!! 厲害 10/18 00:44
11F:推 zupo:不過.Matlab Center裡面人寫的函數名稱是什麼? 10/18 00:48
12F:推 Gwaewluin:http://0rz.tw/b93bS 10/18 01:11
13F:推 gamer:喔,是我記錯了。 10/18 02:41
14F:→ gamer:L-G polynomial的解查表就可以了,自己解太浪費時間了 10/18 02:42
15F:→ gamer:通常函數不要太非線性,取4點5點結果就非常好了 10/18 02:43
16F:→ gamer:http://0rz.tw/d23f2 10/18 02:44
17F:推 zupo:謝謝 ^^ 10/18 10:37
18F:推 guteres:這篇值366銀 10/18 21:16







like.gif 您可能會有興趣的文章
icon.png[問題/行為] 貓晚上進房間會不會有憋尿問題
icon.pngRe: [閒聊] 選了錯誤的女孩成為魔法少女 XDDDDDDDDDD
icon.png[正妹] 瑞典 一張
icon.png[心得] EMS高領長版毛衣.墨小樓MC1002
icon.png[分享] 丹龍隔熱紙GE55+33+22
icon.png[問題] 清洗洗衣機
icon.png[尋物] 窗台下的空間
icon.png[閒聊] 双極の女神1 木魔爵
icon.png[售車] 新竹 1997 march 1297cc 白色 四門
icon.png[討論] 能從照片感受到攝影者心情嗎
icon.png[狂賀] 賀賀賀賀 賀!島村卯月!總選舉NO.1
icon.png[難過] 羨慕白皮膚的女生
icon.png閱讀文章
icon.png[黑特]
icon.png[問題] SBK S1安裝於安全帽位置
icon.png[分享] 舊woo100絕版開箱!!
icon.pngRe: [無言] 關於小包衛生紙
icon.png[開箱] E5-2683V3 RX480Strix 快睿C1 簡單測試
icon.png[心得] 蒼の海賊龍 地獄 執行者16PT
icon.png[售車] 1999年Virage iO 1.8EXi
icon.png[心得] 挑戰33 LV10 獅子座pt solo
icon.png[閒聊] 手把手教你不被桶之新手主購教學
icon.png[分享] Civic Type R 量產版官方照無預警流出
icon.png[售車] Golf 4 2.0 銀色 自排
icon.png[出售] Graco提籃汽座(有底座)2000元誠可議
icon.png[問題] 請問補牙材質掉了還能再補嗎?(台中半年內
icon.png[問題] 44th 單曲 生寫竟然都給重複的啊啊!
icon.png[心得] 華南紅卡/icash 核卡
icon.png[問題] 拔牙矯正這樣正常嗎
icon.png[贈送] 老莫高業 初業 102年版
icon.png[情報] 三大行動支付 本季掀戰火
icon.png[寶寶] 博客來Amos水蠟筆5/1特價五折
icon.pngRe: [心得] 新鮮人一些面試分享
icon.png[心得] 蒼の海賊龍 地獄 麒麟25PT
icon.pngRe: [閒聊] (君の名は。雷慎入) 君名二創漫畫翻譯
icon.pngRe: [閒聊] OGN中場影片:失蹤人口局 (英文字幕)
icon.png[問題] 台灣大哥大4G訊號差
icon.png[出售] [全國]全新千尋侘草LED燈, 水草

請輸入看板名稱,例如:Tech_Job站內搜尋

TOP