| 6 | 1/1 | 返回列表 |
| 查看: 1165 | 回復: 5 | ||
三星2010木蟲 (小有名氣)
|
[求助]
請高人指點一個matlab編程問題,謝謝
|
|
如題:主要是通過求解一個非線性方程組來求最值,一些基本的參數(shù)及方程如下: taof=7.2; Gf=0.576; delta1=0.034; deltaf=0.16; b1=25; E1=500000; t1=5.12; b2=400; E2=10000; t2=20; eql=600; beta=b2*E2*t2/(b1*E1*t1); lamda=((taof^2/2/Gf*(1/E1/t1+b1/b2/E2/t2)))^0.5; lamda1=(lamda^2*2*Gf/(delta1*taof))^0.5; equation1:lamda2/lamda1*tanh(0.5*lamda1*(L-x(1)-x(2)))-(sin(lamda2*x(1))-beta*sin(lamda2*x(2)))/((beta*cos(lamda2*x(2)))-cos(lamda2*x(1)))=0; equation2:lamda2*sinh(L-x(1)-x(2))*tan(lamda2*x(1))+1-lamda1*(cosh(L-x(1)-x(2))+1)=0; 其中x(1)+x(2)<L, L從400變化到1200,L的步長為20,求每一L對應(yīng)的x(1)和x(2),最好通過循環(huán)輸出吧,我用fsolve嘗試了一下,并不理想,請高人給以指點,先謝謝啦! |

木蟲 (小有名氣)

木蟲 (小有名氣)

|
樓主可嘗試編寫m文件如下: function xfsolve clear all;clc; global L; for L=400:20:1200 x0=[L/2,L/2]; k=1; xa(k, =fsolve(@abc,x0)k=k+1; end function f=abc(x) global L; taof=7.2; Gf=0.576; delta1=0.034; deltaf=0.16; b1=25; E1=500000; t1=5.12; b2=400; E2=10000; t2=20; eql=600; beta=b2*E2*t2/(b1*E1*t1); lamda=((taof^2/2/Gf*(1/E1/t1+b1/b2/E2/t2)))^0.5; lamda1=(lamda^2*2*Gf/(delta1*taof))^0.5; lamda2=(lamda^2*2*Gf/(deltaf-delta1)/taof)^0.5; f(1)=lamda2/lamda1*tanh(0.5*lamda1*(L-x(1)-x(2)))-(sin(lamda2*x(1))-beta*sin(lamda2*x(2)))/((beta*cos(lamda2*x(2)))-cos(lamda2*x(1))); f(2)=lamda2*sinh(L-x(1)-x(2))*tan(lamda2*x(1))+1-lamda1*(cosh(L-x(1)-x(2))+1); 其中x0是迭代初值,樓主可依情況適當更改,xa(k, 即為要求的x(1),x(2)。 |
木蟲 (小有名氣)

| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
[考研] 一志愿天津大學化學工藝專業(yè)(081702)315分求調(diào)劑 +10 | yangfz 2026-03-17 | 10/500 |
|
|
[考研] 070303一志愿西北大學學碩310找調(diào)劑 +6 | d如愿上岸 2026-03-12 | 9/450 |
|
|
[考研] 085600材料與化工 +5 | 安全上岸! 2026-03-16 | 5/250 |
|
|
[教師之家] 焦慮 +8 | 水冰月月野兔 2026-03-13 | 12/600 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風 2026-03-17 | 5/250 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 材料專碩274一志愿陜西師范大學求調(diào)劑 +6 | 薛云鵬 2026-03-13 | 6/300 |
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 297求調(diào)劑 +4 | 學海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
[考研] 一志愿山大07化學 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|