 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 914 | 回復: 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]]'; %參數的初始值 x0=[1,0]'; %拉伸應力初始值 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實驗數據[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實驗數據[頻率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]'; 以上是我依照前人經驗編的程序，由于是新手，在程序出錯和修改上解決不了，還望各位蟲友幫幫忙��！ |
回復此樓» 猜你喜歡
請問還有沒有用Latex寫文章的小伙伴們？
已經有0人回復
-
光學工程學碩調劑信息
已經有32人回復
-
物理學I論文潤色/翻譯怎么收費?
已經有204人回復
-
歡迎加入課題組
已經有0人回復
-
散金幣，求好運，祝面上順利！
已經有76人回復
-
華東師范大學芯片設計徐瓏真實水平咋樣
已經有1人回復
-
調劑
已經有0人回復
-
【新加坡】納米電子器件項目組有“聯合培養(yǎng)博士生”名額
已經有0人回復
-
中德博士后2026項目5月30日截止- 國內博士后入站、德國兩年全職
已經有0人回復
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內博后入站德國兩年全職
已經有0人回復
-
【新加坡】納米電子器件項目組有“聯合培養(yǎng)博士生”名額
已經有1人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 合同和編制的啥區(qū)別 求解
已經有10人回復
- 關于一級A標準執(zhí)行問題的一點看法
已經有8人回復
- 2012標書格式問題
已經有3人回復
- 本科畢業(yè)生求解
已經有9人回復
- 關于測ee值的問題
已經有20人回復
- matlab非線性參數擬合問題
已經有7人回復
- 請教一個數據非線性擬合的問題
已經有11人回復
- 求助：非線性方程系數回歸問題？
已經有10人回復
- 關于origin里面的非線性擬合的問題
已經有6人回復
- 求助origin高手 非線性擬合問題
已經有6人回復
- 【求助】求解非線性方程
已經有8人回復
- 請教關于origin中自定義函數非線性回歸擬合問題！
已經有10人回復
dingd
鐵桿木蟲 (職業(yè)作家)
- 計算強帖: 4
- 應助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
2樓2012-02-16 08:23:48
wtisunday
木蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 1742.3
- 散金: 30
- 紅花: 1
- 帖子: 937
- 在線: 568.8小時
- 蟲號: 444484
- 注冊: 2007-10-27
- 性別: GG
- 專業(yè): 高分子材料的加工與成型
|
我的問題是多模態(tài)問題 就是已知t的i個分量和G的i個分量，通過兩個方程求出x的i個分量和y的i個分量，再通過x和y分量和建立優(yōu)化函數，通過優(yōu)化實現a的2i個分量的確定。 比如 i=3 u(i)=a(i)*exp(a(i+4)-t(i))+1/(x(i)+y(i))/G(i).....中間變量表達式 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]';......分量方程組（每一分量均有此關系） Y=sum(x(i)-y(i))/xexpdata.....計算yexpdata Op=sum(Y-yexpdata)^2......優(yōu)化函數 |
3樓2012-02-16 08:51:17
相關版塊跳轉
我要訂閱樓主 wtisunday 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學求調劑 +4 | 奶油草莓. 2026-03-22 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 333求調劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
|
[考研] 291 求調劑 +4 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 5/250 |
|
|
|
[考研] 0703 調劑 +3 | 我可以上岸的對?/a> 2026-03-16 | 6/300 |
|
|
|
[考研] 一志愿西安交通大學材料工程專業(yè) 282分求調劑 +11 | 楓橋ZL 2026-03-18 | 13/650 |
|
|
|
[考研] 308求調劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
|
[考研] 298求調劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
|
[考研] 311求調劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿重慶大學085700資源與環(huán)境總分308求調劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
|
[考研] 南昌大學材料專碩311分求調劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
|
[考研] 307求調劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿蘇州大學材料求調劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] A區(qū)線材料學調劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
|
[考研] 295材料求調劑，一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
|
[考研] 廣西大學家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調劑） +3 | 123阿標 2026-03-17 | 3/150 |
|
|
|
[考研] 0856調劑，是學校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
|
[考研] 材料學碩318求調劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
|
[考研] 085600材料與化工調劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 0703化學調劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[碩博家園] 湖北工業(yè)大學 生命科學與健康學院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|