| 5 | 1/1 | 返回列表 |
| 查看: 1171 | 回復: 5 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
三星2010木蟲 (小有名氣)
|
[求助]
請高人指點一個matlab編程問題,謝謝
|
|
|
如題:主要是通過求解一個非線性方程組來求最值,一些基本的參數及方程如下: 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對應的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)。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311(085601)求調劑 +6 | liziyeyeye 2026-03-28 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 安徽大學專碩生物與醫(yī)藥專業(yè)(086000)324分,英語已過四六級,六級521,求調劑 +9 | 美味可樂雞翅 2026-03-26 | 10/500 |
|
|
[考研] 張芳銘-中國農業(yè)大學-環(huán)境工程專碩-298 +9 | 手機用戶 2026-03-26 | 9/450 |
|
|
[考研] 322求調劑:一志愿湖南大學 材料與化工(085600),已過六級。 +10 | XX小鄧 2026-03-29 | 10/500 |
|
|
[考研] 一志愿西交大080500材料學碩349 +5 | jqx1258 2026-03-31 | 6/300 |
|
|
[考研] 343求調劑 +8 | 愛羈絆 2026-03-28 | 8/400 |
|
|
[考研] 化學工程085602 305分求調劑 +28 | RichLi_ 2026-03-25 | 36/1800 |
|
|
[考研] 一志愿211,335分,0856,求調劑院校和導師 +10 | 傾____蕭 2026-03-27 | 11/550 |
|
|
[考研] 08工科,295,接受跨專業(yè)調劑 +6 | lmnlzy 2026-03-30 | 6/300 |
|
|
[考研] 085600材料與化工329分求調劑 +5 | 這是有太陽哇 2026-03-25 | 5/250 |
|
|
[考研] 285求調劑 +6 | AZMK 2026-03-29 | 9/450 |
|
|
[考研] 297求調劑 +17 | 田洪有 2026-03-26 | 18/900 |
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
[考研] 086502化學工程342求調劑 +6 | 阿姨復古不過 2026-03-27 | 6/300 |
|
|
[考研] 0703化學求調劑,各位老師看看我!! +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[有機交流]
高溫高壓反應求助
10+4
|
chibby 2026-03-25 | 4/200 |
|
|
[考研] 285求調劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
[考研] 085600,材料與化工321分,求調劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 【2026考研調劑】制藥工程 284分 求相關專業(yè)調劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 347求調劑 +4 | L when 2026-03-25 | 4/200 |
|