 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 3186 | 回復(fù): 3 | |||
sunarising銅蟲 (小有名氣)
|
[交流]
【求助】非線性方程組求解 已有2人參與
|
用matlab自帶的fsolve函數(shù)和擬牛頓法等一系列優(yōu)化算法求解結(jié)果都不理想。方程組形式和系數(shù)矩陣算法如下。請教如何求解這個復(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 請教各位大蟲幫幫忙。謝謝啦。 [ Last edited by sunarising on 2010-7-5 at 10:42 ] |
回復(fù)此樓» 猜你喜歡
070300化學求調(diào)劑
已經(jīng)有9人回復(fù)
-
085404 22408 315分
已經(jīng)有6人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有22人回復(fù)
-
347求調(diào)劑
已經(jīng)有10人回復(fù)
-
354求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿華東理工大學，080500學碩，317分，求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿a區(qū)211，085601-307分求調(diào)劑
已經(jīng)有11人回復(fù)
-
284求調(diào)劑
已經(jīng)有7人回復(fù)
-
求調(diào)劑：085600材料與化工，考材科基，總分319
已經(jīng)有18人回復(fù)
-
362求調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- Matlab求解偏微分方程組
已經(jīng)有22人回復(fù)
- matlab求解非線性方程組
已經(jīng)有16人回復(fù)
- 求高人指點用matlab求解非線性方程組，解決了追加100金幣；
已經(jīng)有11人回復(fù)
- matlab的fsove 命令求解非線性方程組
已經(jīng)有6人回復(fù)
- 【求助】多元方程組求解及擬合
已經(jīng)有9人回復(fù)
- 【求助】用mathematica 5.0求解一個非線性方程組失敗，特發(fā)帖求助！
已經(jīng)有5人回復(fù)
- 【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實數(shù)解！
已經(jīng)有18人回復(fù)
- 非線性方程組的迭代法（數(shù)值計算高手請進）
已經(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小時
- 蟲號: 572239
- 注冊: 2008-06-11
- 專業(yè): 力學中的基本問題和方法
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)在的程序下求出的解代入這兩組表達式 如果對這倆組約束方程的重視程度相當?shù)脑? 應(yīng)該在后面乘上合適的權(quán)重 使他們的值大致相當 當然也可以根據(jù)實際情況調(diào)整權(quán)重，側(cè)重某組約束 |
4樓2010-07-05 18:47:12
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 sunarising 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 考研調(diào)劑求助 +4 | 13287130938 2026-03-31 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 272求調(diào)劑，接受跨專業(yè)調(diào)劑！ +3 | 閑魚盧 2026-03-31 | 3/150 |
|
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑 +15 | Micky11223 2026-03-25 | 16/800 |
|
|
|
[考研] 本2一志愿C9-333分，材料科學與工程，求調(diào)劑 +5 | 升升不降 2026-03-31 | 5/250 |
|
|
|
[考研] 370求調(diào)劑 +3 | 080700調(diào)劑 2026-03-30 | 3/150 |
|
|
|
[考研] 317求調(diào)劑 +8 | 十閑wx 2026-03-24 | 8/400 |
|
|
|
[考研] 求收留 +6 | 1943443204 2026-03-28 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +16 | 新時代材料 2026-03-27 | 16/800 |
|
|
|
[考研] 295材料工程專碩求調(diào)劑 +10 | 1428151015 2026-03-27 | 10/500 |
|
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +16 | 李李不服輸 2026-03-26 | 16/800 |
|
|
|
[考研] 環(huán)境科學與工程334分求調(diào)劑 +6 | 王一一依依 2026-03-30 | 8/400 |
|
|
|
[考研] 332求調(diào)劑 +6 | Lyy930824@ 2026-03-29 | 6/300 |
|
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
|
[考研] 311求調(diào)劑 +10 | lin0039 2026-03-26 | 10/500 |
|
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
|
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
|
[有機交流]
高溫高壓反應(yīng)求助
 10+4
|
chibby 2026-03-25 | 4/200 |
|
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|