 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 2015 | 回復: 7 | |||
| 本帖產(chǎn)生 1 個 博學EPI ，點擊這里進行查看 | |||
houbing金蟲 (初入文壇)
|
[交流]
非線性方程組的迭代法（數(shù)值計算高手請進）
|
||
| 我在用matlab求解一組非線性方程組的時候遇到了困難，因為初值選擇不合適，迭代幾乎都不收斂，由于數(shù)據(jù)量較大，沒有辦法對每個初值進行調(diào)整，有沒有一種迭代算法可以對初值沒有要求，我目前使用的是幾個教科書上的算法，牛頓法，不動點迭代，弦割法。期待有高手可以指點迷津，先行謝過！ |
回復此樓» 猜你喜歡
一志愿中國海洋大學，生物學，301分，求調(diào)劑
已經(jīng)有4人回復
-
一志愿武理材料305分求調(diào)劑
已經(jīng)有6人回復
-
材料專碩英一數(shù)二306
已經(jīng)有3人回復
-
085600材料與化工調(diào)劑 324分
已經(jīng)有7人回復
-
311求調(diào)劑
已經(jīng)有3人回復
-
26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線
已經(jīng)有7人回復
-
295求調(diào)劑
已經(jīng)有5人回復
-
354求調(diào)劑
已經(jīng)有5人回復
-
材料專業(yè)求調(diào)劑
已經(jīng)有3人回復
-
求材料調(diào)劑
已經(jīng)有7人回復
gaofeng925
版主 (知名作家)
- 博學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
木蟲 (正式寫手)
清靜的女孩
- 博學EPI: 40
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1867.6
- 帖子: 433
- 在線: 12.5小時
- 蟲號: 494217
- 注冊: 2008-01-12
- 性別: MM
- 專業(yè): 無機合成和制備化學
3樓2010-05-19 10:12:08
houbing
金蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 960.8
- 帖子: 25
- 在線: 3.2小時
- 蟲號: 468800
- 注冊: 2007-11-28
- 性別: MM
- 專業(yè): 固體力學
4樓2010-05-20 10:20:35
houbing
金蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 960.8
- 帖子: 25
- 在線: 3.2小時
- 蟲號: 468800
- 注冊: 2007-11-28
- 性別: MM
- 專業(yè): 固體力學
求助
|
為了方便向大家請教，我把我的程序貼了出來，第一次使用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
6樓2010-05-20 11:28:47
7樓2010-05-20 23:53:16
dingd
鐵桿木蟲 (職業(yè)作家)
- 博學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ā)表 | |
|---|---|---|---|---|
|
[考研] 354求調(diào)劑 +3 | Tyoumou 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學上- 2026-03-14 | 9/450 |
|
|
[碩博家園] 湖北工業(yè)大學 生命科學與健康學院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
|
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習 +6 | Losir 2026-03-12 | 7/350 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學獎覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。，。，。，。i 2026-03-12 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
|
[考研] 0817化學工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
|
[考研] 070303一志愿西北大學學碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|