| 6 | 1/1 | 返回列表 |
| 查看: 1164 | 回復: 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對應的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ā)表 | |
|---|---|---|---|---|
|
[考研] 281求調劑(0805) +6 | 煙汐憶海 2026-03-16 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 344求調劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 286求調劑 +6 | lemonzzn 2026-03-16 | 9/450 |
|
|
[考研] 311求調劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 0817 化學工程 299分求調劑 有科研經歷 有二區(qū)文章 +7 | rare12345 2026-03-18 | 7/350 |
|
|
[考研] 0703化學調劑 ,六級已過,有科研經歷 +10 | 曦熙兮 2026-03-15 | 10/500 |
|
|
[考研] 0703化學調劑 +4 | pupcoco 2026-03-17 | 7/350 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考博] 環(huán)境領域全國重點實驗室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考研] 334求調劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
[考研] 本人考085602 化學工程 專碩 +16 | 不知道叫什么! 2026-03-15 | 18/900 |
|
|
[考研] 332求調劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 一志愿,福州大學材料專碩339分求調劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] 藥學383 求調劑 +3 | 藥學chy 2026-03-15 | 4/200 |
|
|
[考研] 085600調劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 中科大材料與化工319求調劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[考研] 材料080500調劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
[考研] 308 085701 四六級已過求調劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
[考研] 考研調劑 +4 | 芬達46 2026-03-12 | 4/200 |
|
|
[考研] 求調劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|