 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 922 | 回復(fù): 2 | ||
wtisunday木蟲 (正式寫手)
|
[求助]
MATLAB 非線性擬合 fminsearch中包含非線性方程的求解問題
|
|
format short global a1 a2 a3 a4 a5 a6 a7 a8 b1 b2 b3 b4 b5 b6 b7 b8 k0=[[3 3 3 2 2 2 1 1];[1 1 1 1 1 1 1 1]]'; %參數(shù)的初始值 x0=[1,0]'; %拉伸應(yīng)力初始值 t=[0.00423,0.01621,0.06216,0.23841,0.91437,3.50692,14.6856,101.343]'; G=[336370,64246.4,44849.7,10649.7,5151.99,1224.43,63.4524,2.67362]'; %relaxation time,relaxation modulus實(shí)驗(yàn)數(shù)據(jù)[t,G] xexp=[0.03162 0.1 0.31623 1 3.16228 10.5635 20.395]'; yexp=[84834.3 84834.3 79215.4 62322.2 52509.5 32501.6 25570.4]'; %xexp,yexp實(shí)驗(yàn)數(shù)據(jù)[頻率x、拉伸粘度y] [k,fval]=fminsearch(@ObjFunc,k0,x0,xexp,yexp);%非線性擬合 function f=ObjFunc(x0,xexp,yexp) G=[336370,64246.4,44849.7,10649.7,5151.99,1224.43,63.4524,2.67362]'; t=[0.00423,0.01621,0.06216,0.23841,0.91437,3.50692,14.6856,101.343]'; [x]=fsolve(@kineticseqs103,x0); ff=sum(x(:,2)-x(:,1))./xexp; f=(sum(ff'-yexp))^2; function pp=rheolomodel(k,x,g,xexp) global a1 a2 a3 a4 a5 a6 a7 a8 b1 b2 b3 b4 b5 b6 b7 b8 k=[[a1,a2,a3,a4,a5,a6,a7,a8];[b1,b2,b3,b4,b5,b6,b7,b8]]'; q=1+(x(:,1)+2*x(:,2))/3./G.^0.5; u=2.*k(:,1)*exp(2/k(:,2)-1)*(1-1/q.^2)+1/q.^2*(1-0.3/k(:,2).*... ((x(:,1)).^2+2*(x(:,2)).^2)/3/g.^2); pp1=u.*x(:,1)-2*t.*xexp.*(x(:,1)+G)+G.*(u-1)+... 0.3./k(:,2)./G.*(x(:,1))^2; pp2=u.*x(:,2)+t.*xexp.*x(:,2)+G.*(u-1)+... 0.3./k(:,2)./G.*(x(:,2))^2+t.*G.*xexp; pp=[pp1,pp2]'; 以上是我依照前人經(jīng)驗(yàn)編的程序，由于是新手，在程序出錯(cuò)和修改上解決不了，還望各位蟲友幫幫忙�。� |
回復(fù)此樓» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有72人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博士后入站、德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國(guó)內(nèi)博后入站德國(guó)兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團(tuán)隊(duì)招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 合同和編制的啥區(qū)別 求解
已經(jīng)有10人回復(fù)
- 關(guān)于一級(jí)A標(biāo)準(zhǔn)執(zhí)行問題的一點(diǎn)看法
已經(jīng)有8人回復(fù)
- 2012標(biāo)書格式問題
已經(jīng)有3人回復(fù)
- 本科畢業(yè)生求解
已經(jīng)有9人回復(fù)
- 關(guān)于測(cè)ee值的問題
已經(jīng)有20人回復(fù)
- matlab非線性參數(shù)擬合問題
已經(jīng)有7人回復(fù)
- 請(qǐng)教一個(gè)數(shù)據(jù)非線性擬合的問題
已經(jīng)有11人回復(fù)
- 求助：非線性方程系數(shù)回歸問題？
已經(jīng)有10人回復(fù)
- 關(guān)于origin里面的非線性擬合的問題
已經(jīng)有6人回復(fù)
- 求助origin高手 非線性擬合問題
已經(jīng)有6人回復(fù)
- 【求助】求解非線性方程
已經(jīng)有8人回復(fù)
- 請(qǐng)教關(guān)于origin中自定義函數(shù)非線性回歸擬合問題！
已經(jīng)有10人回復(fù)
dingd
鐵桿木蟲 (職業(yè)作家)
- 計(jì)算強(qiáng)帖: 4
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
2樓2012-02-16 08:23:48
wtisunday
木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1742.3
- 散金: 30
- 紅花: 1
- 帖子: 937
- 在線: 568.8小時(shí)
- 蟲號(hào): 444484
- 注冊(cè): 2007-10-27
- 性別: GG
- 專業(yè): 高分子材料的加工與成型
|
我的問題是多模態(tài)問題 就是已知t的i個(gè)分量和G的i個(gè)分量，通過兩個(gè)方程求出x的i個(gè)分量和y的i個(gè)分量，再通過x和y分量和建立優(yōu)化函數(shù)，通過優(yōu)化實(shí)現(xiàn)a的2i個(gè)分量的確定。 比如 i=3 u(i)=a(i)*exp(a(i+4)-t(i))+1/(x(i)+y(i))/G(i).....中間變量表達(dá)式 f1=u(i)*x(i)-t(i)*(x(i)+y(i))+g(i)*((u(i)-1)+expdata; f2=u(i)*y(i)+t(i)*y(i)+y(i)^(2*xexpdata); f=[f1,f2]';......分量方程組（每一分量均有此關(guān)系） Y=sum(x(i)-y(i))/xexpdata.....計(jì)算yexpdata Op=sum(Y-yexpdata)^2......優(yōu)化函數(shù) |
3樓2012-02-16 08:51:17
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wtisunday 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 275求調(diào)劑 +10 | Micky11223 2026-03-25 | 13/650 |
|
|---|---|---|---|---|
|
|
[考研] 289求調(diào)劑 +7 | 新時(shí)代材料 2026-03-27 | 7/350 |
|
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑 +3 | Micky11223 2026-03-25 | 3/150 |
|
|
|
[考研] 0856材料化工調(diào)劑 總分330 +8 | zhubinhao 2026-03-27 | 8/400 |
|
|
|
[考研] 334求調(diào)劑 +3 | 雨清天晴 2026-03-21 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 劉柯@ 2026-03-24 | 4/200 |
|
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 10/500 |
|
|
|
[考研] 材料考研求調(diào)劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 10/500 |
|
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +4 | 15188958825 2026-03-25 | 4/200 |
|
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
|
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601） +5 | WW.' 2026-03-23 | 7/350 |
|
|
|
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|