| 9 | 1/1 | 返回列表 |
| 查看: 1711 | 回復(fù): 8 | ||
tarcy金蟲 (正式寫手)
|
[求助]
求助高手 解一個積分方程
|
|
|
求解磁球的震動方程,里面有一個積分表達(dá)式不太明白怎么得出的。 如圖,第一行方程為被積項表達(dá)式,左端為φ關(guān)于時間t的倒數(shù)項,右側(cè)表達(dá)式包含變量t和φ的表達(dá)式。 第二行方程為一個積分項,同時求極限。 請高手幫解方程,謝謝。 積分公式.jpg |
至尊木蟲 (著名寫手)

銀蟲 (小有名氣)

金蟲 (正式寫手)
至尊木蟲 (著名寫手)
|
首先做變量替換 x=2(Bt- phi), a=A/B. 那么 dx /dt = 2B -2 (A* Sin(x) ), 即 dx/ ( 1-a*Sin(x) ) = 2B* dt. 這是分離變量x, t. 積分兩邊, 右端是2Bt +C, 其中C是任意常數(shù). 左端用萬能公式 sin(x)= 2*tan(x/2) / (1+tan^2(x/2)), 被積函數(shù)可寫成 dx/ ( 1-a*Sin(x) ) = sec^2(x/2) dx / [ 1+tan^2(x/2) -2a*tan(x/2) ] = -2 * d(a-tan(x/2)) / [(1-a^2) + (a-tan(x/2))^2] 這是一個標(biāo)準(zhǔn)積分, 關(guān)于arctan函數(shù)的, 所以左端是 -2/ sqrt(1-a^2) *arctan[ (a-tan(x/2)) / sqrt(1-a^2)]. 因此 tan(x/2) -a = sqrt(1-a^2)* tan(B*sqrt(1-a^2)*t +C). 忐忑地問一下樓主, 是不是 有|A| < |B|的限制條件? 當(dāng)t趨于無窮大時, sqrt(1-a^2)* tan(B*sqrt(1-a^2)*t +C)是周期性的無界的振蕩函數(shù), 因此, tan(x/2)也是振蕩的, 沒有極限. 因此不能用L'Hospital 法則求極限 int_{0}^T sin(x) dt / T. 其實, sin(x(t) )基本上是個周期函數(shù), 所以樓主要求的極限等于一個周期上的積分值, 具體來說就是, int_{- pi/2}^{pi/2} 2*( a+sqrt(1-a^2)*tan(theta) ) / [1 + ( a+sqrt(1-a^2)*tan(theta) )^2 ] d(theta). 經(jīng)過一些雜七雜八的化簡后, 要求計算Integrate[ y / [(1+y^2)(1-a^2+(y-a)^2)], {y, -infinity, infinity} ] Mathematica軟件計算結(jié)果為 -Pi/ (2a), 所以極限值等于 -Pi* sqrt( (B/A)^2-1). 答案有些詭異, 樓主看看符號有沒有 問題, 是否與預(yù)期的相符, 等等. |

至尊木蟲 (著名寫手)
|
哎呀好象Mathematica軟件計算結(jié)果有問題, 那個負(fù)的答案是在某種條件下得到的, http://www.wolframalpha.com/inpu ... 2C+infinity%7D+%5D+ 我看了一下原函數(shù)表達(dá)式, 覺得還是這個答案更合理: Pi* [ B/A - sqrt( (B/A)^2 -1 ) ]. 樓主可以代入一些參數(shù), 驗證合理性. |

木蟲之王 (文學(xué)泰斗)
peterflyer
金蟲 (正式寫手)
至尊木蟲 (著名寫手)

| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求助 +3 | 夢里的無言 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 286求調(diào)劑 +5 | Faune 2026-03-21 | 5/250 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +4 | 葵梓衛(wèi)隊 2026-03-18 | 6/300 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[論文投稿]
申請回稿延期一個月,編輯同意了。但系統(tǒng)上的時間沒變,給編輯又寫郵件了,沒回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|