 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 3184 | 回復(fù): 3 | |||
sunarising銅蟲 (小有名氣)
|
[交流]
【求助】非線性方程組求解 已有2人參與
|
用matlab自帶的fsolve函數(shù)和擬牛頓法等一系列優(yōu)化算法求解結(jié)果都不理想。方程組形式和系數(shù)矩陣算法如下。請(qǐng)教如何求解這個(gè)復(fù)雜的非線性方程組。 系數(shù)矩陣算法如下： function g=g_matrix() global rate n topsize=rate*n*32e-9; basesize=n*32e-9; dtop=rate*32e-9; dbase=32e-9; toploc=[-topsize/2+dtop/2:dtop:topsize/2-dtop/2]; baseloc=[-basesize/2+dbase/2:dbase:basesize/2-dbase/2]; E1=2e4; a0=5e-9; xinf1=rate*2000*32e-9; xinf2=2000*32e-9; c1=-(xinf1+a0)*log((xinf1+a0)^2/(a0)^2)+(xinf1-a0)*log((xinf1-a0)^2/(a0)^2); c2=-(xinf2+a0)*log((xinf2+a0)^2/(a0)^2)+(xinf2-a0)*log((xinf2-a0)^2/(a0)^2); g=zeros(n,n); for i=1:n for j=1:n if (j==i) g(i,j)=-(2*a0*log(4)+c1)/(2*pi*E1*a0*dtop)-(2*a0*log(4)+c2)/(2*pi*E1*(1+0*baseloc(i))*a0*dbase); else g(i,j)=2*(log(abs(xinf1))-log(abs(toploc(i)-toploc(j))))/(pi*E1*dtop)+... 2*(log(abs(xinf2))-log(abs(baseloc(i)-baseloc(j))))/(pi*E1*dbase*(1+0*baseloc(i))); end end end function x=x_vector() global n rate topsize=rate*n*32e-9; basesize=n*32e-9; dtop=rate*32e-9; dbase=32e-9; toploc=[-topsize/2+dtop/2:dtop:topsize/2-dtop/2]; baseloc=[-basesize/2+dbase/2:dbase:basesize/2-dbase/2]; x=zeros(n,1); for i=1:n x(i)=toploc(i)-baseloc(i); x(i)=1e9*x(i); end 請(qǐng)教各位大蟲幫幫忙。謝謝啦。 [ Last edited by sunarising on 2010-7-5 at 10:42 ] |
回復(fù)此樓» 猜你喜歡
266分，求材料相關(guān)專業(yè)調(diào)劑
已經(jīng)有8人回復(fù)
-
材料求調(diào)劑 一志愿哈工大總分298分，前三科223分
已經(jīng)有10人回復(fù)
-
本科211總分289，08工學(xué)真心求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有13人回復(fù)
-
食品工程專碩一志愿中海洋309求調(diào)劑
已經(jīng)有9人回復(fù)
-
一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料科學(xué)與工程調(diào)劑
已經(jīng)有9人回復(fù)
-
材料與化工調(diào)劑一志愿大連海事085600，349
已經(jīng)有5人回復(fù)
-
理學(xué)07化學(xué) 303求調(diào)劑
已經(jīng)有4人回復(fù)
-
甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機(jī)理
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- Matlab求解偏微分方程組
已經(jīng)有22人回復(fù)
- matlab求解非線性方程組
已經(jīng)有16人回復(fù)
- 求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；
已經(jīng)有11人回復(fù)
- matlab的fsove 命令求解非線性方程組
已經(jīng)有6人回復(fù)
- 【求助】多元方程組求解及擬合
已經(jīng)有9人回復(fù)
- 【求助】用mathematica 5.0求解一個(gè)非線性方程組失敗，特發(fā)帖求助！
已經(jīng)有5人回復(fù)
- 【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實(shí)數(shù)解！
已經(jīng)有18人回復(fù)
- 非線性方程組的迭代法（數(shù)值計(jì)算高手請(qǐng)進(jìn)）
已經(jīng)有7人回復(fù)
- 【求助】求教matlab解非線性方程組
已經(jīng)有9人回復(fù)
- 【求助】求解非線性方程
已經(jīng)有8人回復(fù)
- 【分享】PDECOL-經(jīng)典數(shù)值求解非線性偏微分方程組源程序
已經(jīng)有142人回復(fù)
2樓2010-07-05 16:34:19
sunarising
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 408.1
- 帖子: 73
- 在線: 32小時(shí)
- 蟲號(hào): 572239
- 注冊(cè): 2008-06-11
- 專業(yè): 力學(xué)中的基本問題和方法
3樓2010-07-05 18:30:49
sunarising(金幣+30): 2010-07-06 08:41:52
|
增加遺傳代數(shù)試試 或許是精度還不夠 繼續(xù)算很多代以后看看結(jié)果是否合適 另外 f=f+((li(ii)^2+(x(ii)+u)^2)*(li(ii)-Fi(ii)/k)^2-li(ii)^2*l^2)^2*100; f=f+sum((x+u)./li.*Fi)^2*100; 后面都乘了100 可以把現(xiàn)在的程序下求出的解代入這兩組表達(dá)式 如果對(duì)這倆組約束方程的重視程度相當(dāng)?shù)脑? 應(yīng)該在后面乘上合適的權(quán)重 使他們的值大致相當(dāng) 當(dāng)然也可以根據(jù)實(shí)際情況調(diào)整權(quán)重，側(cè)重某組約束 |
4樓2010-07-05 18:47:12
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 sunarising 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +7 | DLkz1314. 2026-03-30 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +5 | min3 2026-03-24 | 6/300 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +6 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
|
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過六級(jí)。 +9 | XX小鄧 2026-03-29 | 9/450 |
|
|
|
[考研] 材料專碩調(diào)劑 +11 | 椰椰。 2026-03-29 | 11/550 |
|
|
|
[考研] 348求調(diào)劑 +6 | 小懶蟲不懶了 2026-03-28 | 6/300 |
|
|
|
[考博] 26申博自薦 +6 | whh869393 2026-03-24 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-24 | 6/300 |
|
|
|
[考研] 調(diào)劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
|
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|
|
[考研] 266求調(diào)劑 +11 | 陽陽哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級(jí)，六級(jí)521，求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
|
[考研] 各位老師您好：本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|