 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 2011 | 回復(fù): 7 | |||
| 本帖產(chǎn)生 1 個 博學(xué)EPI ，點擊這里進(jìn)行查看 | |||
houbing金蟲 (初入文壇)
|
[交流]
非線性方程組的迭代法（數(shù)值計算高手請進(jìn)）
|
||
| 我在用matlab求解一組非線性方程組的時候遇到了困難，因為初值選擇不合適，迭代幾乎都不收斂，由于數(shù)據(jù)量較大，沒有辦法對每個初值進(jìn)行調(diào)整，有沒有一種迭代算法可以對初值沒有要求，我目前使用的是幾個教科書上的算法，牛頓法，不動點迭代，弦割法。期待有高手可以指點迷津，先行謝過！ |
回復(fù)此樓» 猜你喜歡
296求調(diào)劑
已經(jīng)有4人回復(fù)
-
湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士
已經(jīng)有5人回復(fù)
-
085601材料工程專碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
341求調(diào)劑
已經(jīng)有6人回復(fù)
-
考研化學(xué)學(xué)碩調(diào)劑，一志愿985
已經(jīng)有6人回復(fù)
-
302求調(diào)劑
已經(jīng)有9人回復(fù)
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有6人回復(fù)
-
308求調(diào)劑
已經(jīng)有4人回復(fù)
-
化學(xué)工程321分求調(diào)劑
已經(jīng)有14人回復(fù)
-
材料專碩326求調(diào)劑
已經(jīng)有7人回復(fù)
gaofeng925
版主 (知名作家)
- 博學(xué)EPI: 76
- 應(yīng)助: 504 (博士)
- 貴賓: 5.002
- 金幣: 77551.5
- 散金: 370
- 紅花: 20
- 沙發(fā): 2
- 帖子: 9437
- 在線: 689.2小時
- 蟲號: 957495
- 注冊: 2010-02-22
- 性別: GG
- 專業(yè): 環(huán)境工程
- 管轄: 環(huán)境
2樓2010-05-19 10:02:39
lgm19851116
木蟲 (正式寫手)
清靜的女孩
- 博學(xué)EPI: 40
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1867.6
- 帖子: 433
- 在線: 12.5小時
- 蟲號: 494217
- 注冊: 2008-01-12
- 性別: MM
- 專業(yè): 無機(jī)合成和制備化學(xué)
3樓2010-05-19 10:12:08
houbing
金蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 960.8
- 帖子: 25
- 在線: 3.2小時
- 蟲號: 468800
- 注冊: 2007-11-28
- 性別: MM
- 專業(yè): 固體力學(xué)
4樓2010-05-20 10:20:35
houbing
金蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 960.8
- 帖子: 25
- 在線: 3.2小時
- 蟲號: 468800
- 注冊: 2007-11-28
- 性別: MM
- 專業(yè): 固體力學(xué)
求助
|
為了方便向大家請教，我把我的程序貼了出來，第一次使用matlab，對著手冊編了一周，有不夠簡潔的地方還望見諒：） 基本問題就是求解kesai afa gama J0afa J1afa J0gama J1gama（分別為afa gama的零階和一階bessel函數(shù)）七個變量的非線性方程組；共有5328個數(shù)據(jù)點，每個點都需要求解這樣一個方程組，初值只給了kesai的初值，其它變量有顯式的關(guān)系可以通過kesai求解，實際上是利用迭代法求fkesai＝0； j=1,j=2都是收斂的，j＝3就不收斂了 % 不動點迭代 %define constant clear; E=3000000000; rou=1200; K=2500000000; a=0.015; ita=1000000; sampling_rate=10000000; f=(1:5238)*sampling_rate/5238; im=i; %calculate parameters for j=1:5238 Estar(j)=-im*E*ita*f(j)/(E-im*ita*f(j)); end for j=1:5238 kesai0(j)=sqrt(rou*f(j)^2/Estar(j)); end for j=1:5238 miu(j)=3*K*f(j)*ita*im/(9*K*(1+im*f(j)*ita/E)-im*f(j)*ita); lamda(j)=K-2/3*miu(j); end %initial value of variables for j=1:5238 kesai(j)=kesai0(j); afa(j)=sqrt(rou*f(j)^2/(lamda(j)+2*miu(j))-kesai(j)^2); gama(j)=sqrt(rou*f(j)^2/miu(j)-kesai(j)^2); J0afa(j)=besselj(0,afa(j)*a); J1afa(j)=besselj(1,afa(j)*a); J0gama(j)=besselj(0,gama(j)*a); J1gama(j)=besselj(1,gama(j)*a); fkesai(j)=2*afa(j)/a*(gama(j)^2+kesai(j)^2)*J1afa(j)*J1gama(j)-(gama(j)^2-kesai(j)^2)*J0afa(j)*J1gama(j)-4*kesai(j)*afa(j)*gama(j)*J1afa(j)*J0gama(j); j %iterative n=1; while abs(fkesai(j))>0.0001&(n<=10000) %不動點迭代from fkesai＝0 kesai(j)=(2*afa(j)/a*(gama(j)^2+kesai(j)^2)*J1afa(j)*J1gama(j)-(gama(j)^2-kesai(j)^2)*J0afa(j)*J1gama(j))/(4*afa(j)*gama(j)*J1afa(j)*J0gama(j)); afa(j)=sqrt(rou*f(j)^2/(lamda(j)+2*miu(j))-kesai(j)^2); gama(j)=sqrt(rou*f(j)^2/miu(j)-kesai(j)^2); J0afa(j)=besselj(0,afa(j)); J1afa(j)=besselj(1,afa(j)); J0gama(j)=besselj(0,gama(j)); J1gama(j)=besselj(1,gama(j)); fkesai(j)=2*afa(j)/a*(gama(j)^2+kesai(j)^2)*J1afa(j)*J1gama(j)-(gama(j)^2-kesai(j)^2)*J0afa(j)*J1gama(j)-4*kesai(j)*afa(j)*gama(j)*J1afa(j)*J0gama(j); n=n+1; abs(fkesai) end end |
5樓2010-05-20 11:07:37
hxz0407
金蟲 (小有名氣)
- 博學(xué)EPI: 1
- 應(yīng)助: 0 (幼兒園)
- 金幣: 840.2
- 散金: 60
- 紅花: 2
- 帖子: 193
- 在線: 275.2小時
- 蟲號: 206029
- 注冊: 2006-03-03
- 專業(yè): 環(huán)境工程
6樓2010-05-20 11:28:47
hxz0407
金蟲 (小有名氣)
- 博學(xué)EPI: 1
- 應(yīng)助: 0 (幼兒園)
- 金幣: 840.2
- 散金: 60
- 紅花: 2
- 帖子: 193
- 在線: 275.2小時
- 蟲號: 206029
- 注冊: 2006-03-03
- 專業(yè): 環(huán)境工程
7樓2010-05-20 23:53:16
dingd
鐵桿木蟲 (職業(yè)作家)
- 博學(xué)EPI: 1
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
8樓2011-04-18 10:33:40
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 houbing 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學(xué)工程321分求調(diào)劑 +11 | 大米飯！ 2026-03-15 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +5 | yangfz 2026-03-17 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
|
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 344求調(diào)劑 +3 | knight344 2026-03-16 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 285 求調(diào)劑 資源與環(huán)境 一志愿北京化工大學(xué) +3 | 未名考生 2026-03-10 | 3/150 |
|
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +9 | jiaanl 2026-03-11 | 9/450 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
|
[考研] 283求調(diào)劑，材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|