| 24小時(shí)熱門(mén)版塊排行榜 |
| 3 | 1/1 | 返回列表 |
| 查看: 838 | 回復(fù): 2 | |||
周亦良金蟲(chóng) (小有名氣)
|
[求助]
matlab求兩個(gè)自變量的最值 已有1人參與
|
|
clear clc T=0.3; Vsw=2500; Vvw=3900; Vs=100; Vv=187; yc=25; yd=16; yst=19; yw=10; kh1=0.3; H=10; ru=0.2; b=0.2*H; h=0.75*H; B=30*pi/180; D=B/2; C=10*pi/180; A=0:0.01*pi/180:pi/2; t=0:0.001:0.3; yy=((h/H)^2*yst+(1-(h/H)^2)*yd) kh=yst*kh1/(yy*(1-ru)) kv=0.5*kh; m1=T*Vs*kh*(1+tan(A)*tan(C))/(2*pi^2*H^2*tan(A))*(2*pi*H*cos(2*pi*(t/T-H/(T*Vs)))+T*Vs*(sin(2*pi*(t/T-H/(T*Vs)))-sin(2*pi*t/T))); m2=T*Vv*kv*(1+tan(A)*tan(C))/(2*pi^2*H^2*tan(A))*(2*pi*H*cos(2*pi*(t/T-H/(T*Vv)))+T*Vv*(sin(2*pi*(t/T-H/(T*Vv)))-sin(2*pi*t/T))); Ke=(1+tan(A)*tan(C))*sin((A+B))/tan(A)/cos((A+B+D-C))-m1*cos((A+B))/cos((A+B+D-C))-m2*sin((A+B))/cos((A+B+D-C)) 麻煩用matlab寫(xiě)出求得Ke的最小值(Ke為正的最小值)程序,以及相應(yīng)A,t的值 發(fā)自小木蟲(chóng)Android客戶(hù)端 |
金蟲(chóng) (小有名氣)
|
clear; clc format long g T=0.3; Vsw=2500; Vvw=3900; Vs=100; Vv=187; yc=25; yd=16; yst=19; yw=10; kh1=0.3; H=10; ru=0.2; b=0.2*H; h=0.75*H; B=30*pi/180; D=B/2; C=10*pi/180; % A=0:0.01*pi/180:pi/2; % t=0:0.001:0.3; syms A t yy=((h/H)^2*yst+(1-(h/H)^2)*yd); kh=yst*kh1/(yy*(1-ru)); kv=0.5*kh; m1=T*Vs*kh*(1+tan(A)*tan(C))/(2*pi^2*H^2*tan(A))*(2*pi*H*cos(2*pi*(t/T-H/(T*Vs)))+T*Vs*(sin(2*pi*(t/T-H/(T*Vs)))-sin(2*pi*t/T))); m2=T*Vv*kv*(1+tan(A)*tan(C))/(2*pi^2*H^2*tan(A))*(2*pi*H*cos(2*pi*(t/T-H/(T*Vv)))+T*Vv*(sin(2*pi*(t/T-H/(T*Vv)))-sin(2*pi*t/T))); Ke = (1+tan(A)*tan(C))*sin((A+B))/tan(A)/cos((A+B+D-C))-m1*cos((A+B))/cos((A+B+D-C))-m2*sin((A+B))/cos((A+B+D-C)); Kem = matlabFunction(Ke); [x,fval,exitflag] = fminsearch(@(x) Kem(x(1),x(2)),[0.5, 1.3]); A0 = x(1) t0 = x(2) Ke_min = fval 計(jì)算結(jié)果: A0 = 0.281986095843263 t0 = 1.33097106429981 Ke_min 2.04515149602023 求最值,可根據(jù)具體問(wèn)題,選用 fmincon 、fminunc 、fminsearch 等函數(shù)。 |
金蟲(chóng) (小有名氣)
送紅花一朵 |
2樓: Originally posted by zyh700600 at 2022-05-05 11:58:26 clear; clc format long g T=0.3; Vsw=2500; Vvw=3900; Vs=100; Vv=187; yc=25; yd=16; yst=19; yw=10; kh1=0.3; H=10; ru=0.2; b=0.2*H; h=0.75*H; B=30*pi/180; D=B/2; C=10*pi/180; ... 你好,很高興你的解答,但我改變kh1=0.6時(shí),會(huì)顯示超過(guò)了函數(shù)計(jì)算的最大數(shù)目- 請(qǐng)?jiān)龃?MaxFunEvals 選項(xiàng)。請(qǐng)問(wèn)應(yīng)該如何調(diào)整?非常感謝 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化工學(xué)碩306求調(diào)劑 +7 | 42838695 2026-03-12 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +5 | yfihxh 2026-03-09 | 5/250 |
|
|
[考研] 0857環(huán)境調(diào)劑 +5 | 熠熠_11 2026-03-10 | 5/250 |
|
|
[考研] 2026考研求調(diào)劑-材料類(lèi)-本科211一志愿985-初試301分 +10 | 蟲(chóng)友233 2026-03-07 | 10/500 |
|
|
[考研] 復(fù)試調(diào)劑 +6 | 呼呼?~+123456 2026-03-08 | 8/400 |
|
|
[考研] 物理學(xué)求調(diào)劑 +4 | sx宋 2026-03-05 | 4/200 |
|
|
[考研] 0856求調(diào)劑 +3 | squirtle11 2026-03-07 | 3/150 |
|
|
[考研] 293一志愿華東理工 0817化學(xué)工程與技術(shù) 調(diào)劑 +5 | fjj0912 2026-03-07 | 5/250 |
|
|
[考研] 0701-322 求調(diào)劑 +3 | jiliuxian 2026-03-06 | 8/400 |
|
|
[考研] 288求調(diào)劑085600材料與化工 +13 | Daunrin 2026-03-07 | 15/750 |
|
|
[考研] 求調(diào)劑 +3 | 拾柒12。 2026-03-08 | 3/150 |
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
[考研] 求調(diào)劑,不管什么專(zhuān)業(yè),我是可塑造的人才一枚,希望遇到知己老師撈撈我 +4 | 13102137290 2026-03-06 | 5/250 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)0856材料與化工,前三科206,總分283,求調(diào)劑 +7 | 26考研求調(diào)劑 2026-03-06 | 7/350 |
|
|
[考研] 306求調(diào)劑 +7 | Bahati 2026-03-05 | 7/350 |
|
|
[考研] 哈爾濱理工大學(xué)2026年研究生調(diào)劑,材料科學(xué)與化學(xué)工程學(xué)院研究生調(diào)劑 +3 | xinliu866 2026-03-06 | 3/150 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-06 | 4/200 |
|
|
[考研] 287求調(diào)劑 +3 | 看看我. 2026-03-05 | 6/300 |
|
|
[考研] 085600,一志愿鄭州大學(xué),280分求調(diào)劑 +7 | Wuqi725 2026-03-05 | 7/350 |
|
|
[考研] 復(fù)試調(diào)劑 +5 | 呼呼?~+123456 2026-03-05 | 5/250 |
|