 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3181 | 回復(fù): 31 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
businiao412金蟲 (正式寫手)
|
[求助]
請(qǐng)教一個(gè)非線性方程非線性項(xiàng)轉(zhuǎn)化為多項(xiàng)式的方法 已有2人參與
|
|
|
RT,樓主近日做振動(dòng)問(wèn)題時(shí)遇到一個(gè)二階非線性微分方程如圖片1，看文獻(xiàn)中其有如下解法:在運(yùn)用攝動(dòng)法之前，先將其中的非線性項(xiàng)（1-q）3/2運(yùn)用三階最小二乘擬合的方法在區(qū)間1-Q<q<1上進(jìn)行展開如圖片2，及圖片3，請(qǐng)問(wèn)這是怎么做到的呢？我查了很多資料就是找不到這個(gè)方法的出處，本人數(shù)學(xué)知識(shí)較為薄弱，望不吝賜教。 圖片1-3如下圖所示： 方程1.png 2.png 3.png [ Last edited by businiao412 on 2014-8-10 at 15:16 ] |
回復(fù)此樓» 猜你喜歡
南京大學(xué)化學(xué)調(diào)劑
已經(jīng)有14人回復(fù)
-
一志愿北京工業(yè)大學(xué)，324分求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑
已經(jīng)有14人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑
已經(jīng)有6人回復(fù)
-
26申博自薦
已經(jīng)有6人回復(fù)
-
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分
已經(jīng)有10人回復(fù)
-
求調(diào)劑 有機(jī)化學(xué)考研356分
已經(jīng)有5人回復(fù)
-
300求調(diào)劑，材料科學(xué)英一數(shù)二
已經(jīng)有10人回復(fù)
-
材料與化工328分調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- matlab解非線性方程組
已經(jīng)有3人回復(fù)
- matlab解非線性方程組程序糾錯(cuò)與正確運(yùn)行
已經(jīng)有6人回復(fù)
- 非線性方程的線性化問(wèn)題，大神們進(jìn)來(lái)指教下！�。�！感謝了！�。。。。。。�！
已經(jīng)有9人回復(fù)
- 1stopt/MATLAB 解非線性方程組問(wèn)題
已經(jīng)有4人回復(fù)
- 一個(gè)關(guān)于有約束非線性方程組的求解問(wèn)題
已經(jīng)有7人回復(fù)
- 求大神幫忙擬合一個(gè)非線性方程，求出模型參數(shù)
已經(jīng)有15人回復(fù)
- 求助解復(fù)雜非線性方程組的好的方法
已經(jīng)有24人回復(fù)
- 用matlab求解一個(gè)非線性方程組的解
已經(jīng)有4人回復(fù)
- 非線性方程求解問(wèn)題，求助各位兄弟姐妹們~~
已經(jīng)有7人回復(fù)
- sos!!1stopt解非線性方程組，一直出現(xiàn)常數(shù)定義錯(cuò)誤的問(wèn)題！求大神幫助！
已經(jīng)有12人回復(fù)
- MATLAB求解非線性方程組
已經(jīng)有5人回復(fù)
- 求Matlab解非線性方程代碼
已經(jīng)有12人回復(fù)
- matlab無(wú)約束非線性方程求最小值的問(wèn)題
已經(jīng)有7人回復(fù)
- 【求助】非線性方程組的求解問(wèn)題
已經(jīng)有6人回復(fù)
- 求助matlab 非線性方程組全部解！
已經(jīng)有6人回復(fù)
- 使用fsolve求解非線性方程問(wèn)題
已經(jīng)有9人回復(fù)
- 需要換迭代方法嗎？
已經(jīng)有7人回復(fù)
- 非線性數(shù)據(jù)擬合的數(shù)學(xué)原理�。。。。�？？？？
已經(jīng)有30人回復(fù)
- 急切求助！一個(gè)mathematica求解一個(gè)簡(jiǎn)單一元非線性方程的問(wèn)題
已經(jīng)有4人回復(fù)
- 1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？
已經(jīng)有13人回復(fù)
- matlab求解非線性方程組，求助！
已經(jīng)有6人回復(fù)
- 【求助】用mathematica 5.0求解一個(gè)非線性方程組失敗，特發(fā)帖求助！
已經(jīng)有5人回復(fù)
- 非線性方程組的迭代法（數(shù)值計(jì)算高手請(qǐng)進(jìn)）
已經(jīng)有7人回復(fù)
- 【求助】求解非線性方程
已經(jīng)有8人回復(fù)
businiao412
金蟲 (正式寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1431.3
- 紅花: 9
- 帖子: 469
- 在線: 58.9小時(shí)
- 蟲號(hào): 1496676
- 注冊(cè): 2011-11-17
- 專業(yè): 石油天然氣開采
9樓2014-08-13 11:34:35
herowolf
木蟲 (正式寫手)
- 應(yīng)助: 25 (小學(xué)生)
- 金幣: 3766.3
- 紅花: 8
- 帖子: 731
- 在線: 381.1小時(shí)
- 蟲號(hào): 3283264
- 注冊(cè): 2014-06-19
- 性別: GG
- 專業(yè): 流體力學(xué)
2樓2014-08-10 16:53:53
guojiashun
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 10 (幼兒園)
- 金幣: 7191.2
- 散金: 82
- 紅花: 12
- 帖子: 940
- 在線: 369.1小時(shí)
- 蟲號(hào): 933003
- 注冊(cè): 2009-12-25
- 性別: GG
- 專業(yè): 機(jī)械動(dòng)力學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
businiao412: 金幣+15, ★★★★★最佳答案, 衷心感謝您細(xì)致透徹的講解，有種豁然開朗的感覺(jué)，我按照您的思路試做一下，有問(wèn)題再向您請(qǐng)教。 2014-08-11 14:57:17
感謝參與，應(yīng)助指數(shù) +1
businiao412: 金幣+15, ★★★★★最佳答案, 衷心感謝您細(xì)致透徹的講解，有種豁然開朗的感覺(jué)，我按照您的思路試做一下，有問(wèn)題再向您請(qǐng)教。 2014-08-11 14:57:17
|
對(duì)你的問(wèn)題理解不是十分準(zhǔn)確。僅僅就我的理解說(shuō)一下： 攝動(dòng)法的基本思想是把響應(yīng)q看做多個(gè)不同量級(jí)的攝動(dòng)的疊加，即：q(t,ε )=u0(t)+ε*u1(t)+ε ^2*u2(t)+......【1】（多尺度法還要把考慮不同的時(shí)間尺度，q(t,ε)=u0(t1,t2,...)+ε*u1(t1,t2,...)+ε ^2*u2(t1,t2,...)+.....），然后把q的近似表達(dá)式【1】代入原方程，通過(guò)整理和推導(dǎo)，根據(jù)方程兩邊 小參數(shù) ε 的同次冪系數(shù)相等的條件，把原方程轉(zhuǎn)化為分別關(guān)于u0、u1、u2.....的微分方程，然后分別逐次求得各攝動(dòng)u0、u1、u2........最后，疊加可得方程q的近似解析解。 這個(gè)非線性微分方程的非線性項(xiàng)是（1-q）^3/2，而在運(yùn)用攝動(dòng)法把式【1】代入原方程中進(jìn)行整理時(shí)，由于（1-q）^3/2這一項(xiàng)帶有分?jǐn)?shù)次冪，故【1】代入后其展開是非常困難的，即使展開其中也帶有 ε 的分?jǐn)?shù)次冪項(xiàng)，也就無(wú)法利用 ε 的同次冪系數(shù)相等的條件來(lái)求解（后面的套路步驟走不下去了）。所以，一般在求解此類非線性微分方程時(shí)，用多項(xiàng)式來(lái)近似表達(dá)諸如（1-q）^3/2這種難于處理的項(xiàng)，這樣當(dāng)式【1】代入后展開、整理就變得很簡(jiǎn)便了。類似的情況如齒輪動(dòng)力學(xué)理論中，一個(gè)主要的非線性因素是齒側(cè)間隙，間隙作用方式用分段函數(shù)來(lái)表達(dá)，而分段函數(shù)的存在也給攝動(dòng)法、多尺度法等的運(yùn)用帶來(lái)困難，數(shù)學(xué)上難于展開處理，故也有學(xué)者用多項(xiàng)式來(lái)近似這個(gè)分段函數(shù)（當(dāng)然還有其它處理方法，這個(gè)還要視造成困難的項(xiàng)的具體函數(shù)性質(zhì)）。 這里提前預(yù)估了一個(gè)擬合區(qū)間q∈[1-Q, 1]，Why？第一，方程首先要有穩(wěn)態(tài)解；第二，也是為了數(shù)學(xué)處理上的便利性。（1）根據(jù)（1-q）^3/2可知，q必須小于等于1，否則出現(xiàn)復(fù)數(shù)項(xiàng)，方程無(wú)穩(wěn)態(tài)解，故已知q的上限；（2）采用多項(xiàng)式來(lái)擬合（1-q）^3/2，取多少階合適？我們當(dāng)然是希望階數(shù)越低越好。階數(shù)越高，式【1】代入后的數(shù)學(xué)展開越復(fù)雜繁瑣。根據(jù)數(shù)學(xué)知識(shí)易知擬合區(qū)間越小，就可以用越低的階數(shù)來(lái)獲得足夠的擬合精度，所以根據(jù)已知條件提前劃定一個(gè)盡量小的區(qū)間十分有益和必要。接下來(lái)要確定q的下限，其中Q是q所能預(yù)期到的最大的峰值（上限）和谷值（下限）之間的差值（peak-to-peak value，我想這里的Q可以通過(guò)數(shù)值求解、實(shí)驗(yàn)或者實(shí)際觀測(cè)得到），那么已知q的上限是1，下限自然就是1-Q了。 最后，文中確定用三次多項(xiàng)式來(lái)擬合，接著就要確定多項(xiàng)式的系數(shù)a0、a1、a2、a3。這里采用“最小二乘法”來(lái)確定，最小二乘法是一種數(shù)值方法，必須基于已知的數(shù)據(jù)進(jìn)行擬合。上面已經(jīng)得到了擬合區(qū)間q∈[1-Q, 1]，那么在這個(gè)區(qū)間內(nèi)取一些點(diǎn)，如1-Q, 1-9*Q/10, 1-8*Q/10, 1-7*Q/10, ......, 1，依次代入（1-q）^3/2中，可得相應(yīng)值Q^3/2, (9*Q/10)^3/2, (8*Q/10)^3/2, (7*Q/10)^3/2, ......, 0，已知了函數(shù)（1-q）^3/2的一組自變量和因變量，套用“最小二乘法”的套路解得多項(xiàng)式的系數(shù)a0、a1、a2、a3也就不難了吧？ 總結(jié)，不論是攝動(dòng)法，還是多項(xiàng)式擬合，都是無(wú)限逼近思想的體現(xiàn)。攝動(dòng)法本身得到的就是方程的無(wú)限近似解，因此用多項(xiàng)式無(wú)限近似方程中復(fù)雜項(xiàng)也就變得無(wú)可厚非了。 PS：周日下午無(wú)事，多敲了些廢話，純屬個(gè)人見(jiàn)解，限于水平，有誤之處請(qǐng)各位指正。 |
3樓2014-08-10 19:43:10
businiao412
金蟲 (正式寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1431.3
- 紅花: 9
- 帖子: 469
- 在線: 58.9小時(shí)
- 蟲號(hào): 1496676
- 注冊(cè): 2011-11-17
- 專業(yè): 石油天然氣開采
4樓2014-08-11 14:59:13
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +3 | 大饞小子 2026-03-29 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 材料與化工272求調(diào)劑 +14 | 阿斯蒂芬2004 2026-03-28 | 14/700 |
|
|
|
[考研] 343求調(diào)劑085601 +3 | 要努力學(xué)習(xí)x 2026-03-29 | 3/150 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5 | hanamiko 2026-03-29 | 5/250 |
|
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
|
[考研] 0703化學(xué) +11 | 妮妮ninicgb 2026-03-27 | 11/550 |
|
|
|
[考研] 一志愿211 初試270分 求調(diào)劑 +8 | 谷雨上岸 2026-03-23 | 9/450 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求導(dǎo)師收 +9 | 天天好運(yùn)來(lái)上岸?/a> 2026-03-24 | 10/500 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
|
[考研] 材料考研求調(diào)劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
|
[考研] 化工專碩求調(diào)劑 +3 | question挽風(fēng) 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒(méi)煩惱 2026-03-23 | 4/200 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|