 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2028 | 回復(fù): 7 | |||
| 本帖產(chǎn)生 1 個 博學(xué)EPI ，點擊這里進行查看 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
houbing金蟲 (初入文壇)
|
[交流]
非線性方程組的迭代法（數(shù)值計算高手請進）
|
||
| 我在用matlab求解一組非線性方程組的時候遇到了困難，因為初值選擇不合適，迭代幾乎都不收斂，由于數(shù)據(jù)量較大，沒有辦法對每個初值進行調(diào)整，有沒有一種迭代算法可以對初值沒有要求，我目前使用的是幾個教科書上的算法，牛頓法，不動點迭代，弦割法。期待有高手可以指點迷津，先行謝過！ |
回復(fù)此樓» 猜你喜歡
310求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有17人回復(fù)
-
材料工程專碩求調(diào)劑
已經(jīng)有3人回復(fù)
-
337求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工328分調(diào)劑
已經(jīng)有9人回復(fù)
-
318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有6人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有4人回復(fù)
-
327求調(diào)劑
已經(jīng)有4人回復(fù)
lgm19851116
木蟲 (正式寫手)
清靜的女孩
- 博學(xué)EPI: 40
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1867.6
- 帖子: 433
- 在線: 12.5小時
- 蟲號: 494217
- 注冊: 2008-01-12
- 性別: MM
- 專業(yè): 無機合成和制備化學(xué)
3樓2010-05-19 10:12:08
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
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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +4 | 兔兔兔111223314 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +6 | 家佳佳佳佳佳 2026-03-29 | 6/300 |
|
|
|
[考研] 290求調(diào)劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +7 | 26要上岸 2026-03-28 | 7/350 |
|
|
|
[考研] 環(huán)境工程 085701，267求調(diào)劑 +6 | minht 2026-03-29 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +9 | 張zz111 2026-03-27 | 10/500 |
|
|
|
[考研] 求收留 +5 | 1943443204 2026-03-28 | 5/250 |
|
|
|
[考研] 343求調(diào)劑 +6 | 愛羈絆 2026-03-29 | 6/300 |
|
|
|
[考研] 348求調(diào)劑 +5 | 小懶蟲不懶了 2026-03-28 | 5/250 |
|
|
|
[考研] 調(diào)劑求院校招收 +6 | 鶴鯨鴿 2026-03-28 | 6/300 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
|
[考研] 材料考研求調(diào)劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機 2026-03-25 | 4/200 |
|
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|