 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 913 | 回復(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ù)此樓» 猜你喜歡
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有32人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有240人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有76人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有1人回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)求調(diào)劑 +6 | 苑豆豆 2026-03-20 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 0854電子信息求調(diào)劑 324 +3 | Promise-jyl 2026-03-23 | 3/150 |
|
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +6 | 紙魚ly 2026-03-21 | 6/300 |
|
|
|
[考研]
求調(diào)劑材料學(xué)碩080500，總分289分
 5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
|
[考研] 0703 調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-16 | 6/300 |
|
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +11 | 冷笙123 2026-03-17 | 11/550 |
|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
|
[考研] 材料工程（專）一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|