 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 593 | 回復(fù): 0 | ||
thooo1985新蟲 (初入文壇)
|
[求助]
MATLAB二元非線性優(yōu)化問(wèn)題
|
|
麻煩老師幫忙看一下程序的問(wèn)題。中間函數(shù)相互關(guān)聯(lián)，不知道該怎么編寫函數(shù)了，程序如下。謝謝 function [w,x,n]=tansuo3(W,x1,x2,x4) %初始輸入值% % w=20; % x(1,1)=180; % x(1,2)=4; % x(1,4)=3; x(1,1)=x1; x(1,2)=x2; x(1,4)=x4; %第一次方程求解得到的值% y(1,1)=2.2171+0.021226.*x(1,1)-0.47545.*x(1,2); y(1,2)=-18.09883+0.32313.*x(1,1)-1.00499.*x(1,2)-7.84734*10.^(-4).*x(1,1).^2; y(1,3)=y(1,1); %方程組（y(1,1),y(1,2)）與方程組(y(1,3),y(1,4))的關(guān)系% x(1,3)=(%y(1,3)-4.12488+1.5333.*x(1,4)-0.13043.*x(1,4).^2)./0.031; y(1,4)=0.61721+0.058066.*x(1,3)+0.01073.*x(1,4)-5.5*10.^(-3).*x(1,3).*x(1,4); %第一次計(jì)算滿足條件的最小n值% w(1,1)=W; n(1,1)=(3*w(1,1)+3*y(1,2)/2-y(1,4)*3/2)/(2*y(1,2)); n(1,1)=ceil(n(1,1)); %返回計(jì)算最終的w值% if n>2 w(1,1)=y(1,4)/2+y(1,2)/3+(n(1,1)-2)*y(1,2)*2/3+y(1,2)/2; elseif n==2 w(1,1)=y(1,4)/2+y(1,2)/3+y(1,2)/2; elseif n<2 w(1,1)=y(1,4)/2; end %以第一次計(jì)算的w為條件，計(jì)算第二次至第i次的x(i,1),x(i,2),x(i,3),x(i,4)和最小值n% k=y(1,1); i=1; while k<=w %循環(huán)終止條件% i=i+1; %10是假定的值，循環(huán)可能在10次內(nèi)停止% %第i次方程組的關(guān)系，與第一次的關(guān)系相同% % y(i,1)=2.2171+0.021226.*x(i,1)-0.47545.*x(i,2); % y(i,2)=-18.09883+0.32313.*x(i,1)-1.00499.*x(i,2)-7.84734*10.^(-4).*x(i,1).^2; % y(i,3)=y(i,1); % x(i,3)=(y(i,3)-4.12488+1.5333.*x(i,4)-0.13043.*x(i,4).^2)./0.031; % y(i,4)=0.61721+0.058066.*x(i,3)+0.01073.*x(i,4)-5.5*10.^(-3).*x(i,3).*x(i,4); %%%%%%%%%%%%%%%%%%%%%%%第二次循環(huán)時(shí)的電流未知 w(i,1)=w(i-1,1)-y(i-1,2)/2; %第i次w的值與第i-1次的關(guān)系% [X fval]=fmincon(@GAsovle1,[100,3,95,3.2],[164,3.6,94,2.6],[206,6.4,136,5.4]) %求滿足x(i,1),x(i,2),x(i,3),x(i,4)范圍要求的最小的n的值% n(i,1)=ceil(n(i,1)) %對(duì)n取整% %取整后的n帶入w的計(jì)算公式，求解新的w值% y(1)=; y(2)=(-18.09883+0.32313.*X(1)-1.00499.*X(2)-7.84734*10.^(-4).*X(1).^2); y(3)=; y(4)=(0.61721+0.058066.*X(3)+0.01073.*X(4)-5.5*10.^(-3).*X(3).*X(4)); if n>2 w(i,1)=y(i,4)/2+y(i,2)/3+(n(i,1)-2)*y(i,2)*2/3+y(i,2)/2; elseif n==2 w(i,1)=y(i,4)/2+y(i,2)/3+y(i,2)/2; elseif n<2 w(i,1)=y(i,4)/2; end %輸出第i次計(jì)算的w,n,x的值 w(i,1) n(i,1) x(i,1)=X(1); x(i,2)=X(2); x(i,3)=X(3); x(i,4)=X(4); k=k+y(i,1) end |
回復(fù)此樓» 猜你喜歡
267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿華中科技大學(xué)，080502，354分求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿武理材料305分求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料專碩英一數(shù)二306
已經(jīng)有3人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有7人回復(fù)
-
311求調(diào)劑
已經(jīng)有3人回復(fù)
-
26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線
已經(jīng)有7人回復(fù)
-
295求調(diào)劑
已經(jīng)有5人回復(fù)
-
354求調(diào)劑
已經(jīng)有5人回復(fù)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 thooo1985 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +7 | llllkkkhh 2026-03-18 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-18 | 3/150 |
|
|
|
[考研] 354求調(diào)劑 +3 | Tyoumou 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿天大材料與化工（085600）總分338 +3 | 蔡大美女 2026-03-13 | 3/150 |
|
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +9 | yangfz 2026-03-17 | 9/450 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|