优化为一个函数输出多个变量有不同的范围。

                             %定义T1, T2, P1, P2范围:
                            
                             T1 = 20:45;
                            
                             T2 = 30:6 0;
                            
                             %值P1和P2
                            
                             P1 =…;
                            
                             P2 =…;
                            
                             %你也需要一个矩阵将你所有的结果
                            
                             %结果等于长度(T1) * (T2) *长度(P1) *长度(P2)。
                            
                             %定义一个计数器来存储你的结果:
                            
                             counter = 1;
                            
                             为i = 1:长度(T1)
                            
                             为j = 1:长度(T2)
                            
                             为k = 1:长度(P1)
                            
                             为m = 1:长度(P2)
                            
                             NetProfit(柜台)= SteamPlant (i, j, k、m);
                            
                             Input_data(柜台)=(计数器,i, j, k、m];
                            
                             counter = counter + 1;
                            
                             结束
                            
                             结束
                            
                             结束
                            
                             结束

这将是非常耗时大值。它将循环通过所有值恒定P1, P2 T1, T2然后改变P1,循环通过P2。它一直持续到结束所有T1值。

2的评论
显示1年长的评论藏1年长的评论

约瑟夫Woloschek 2019年2月13日

当我将这种类型的代码输入脚本,它说我有重复。附加是我的完整代码,前两块是我的尝试,使某种循环功能和你给的答案。下面这是我的主要功能。

当我试着把自己的脚本文件中的代码和下一个函数,它说它不承认原SteamPlant功能我想打电话。

我真的新MATLAB抱歉如果我犯了一个很明显的错误。我非常感谢帮助!

                                  % T1 = (40、45);
                                 
                                  % boiler_T_range = [400450500550600];
                                 
                                  % boiler_P_range = (5、6、7、8、9、10、11、12、13、14、15);
                                 
                                  % turbine_P_range = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12、13、14、15);
                                 
                                  % turbine1_eff = [.80、.85 .90 .95];
                                 
                                  % turbine2_eff = [.80、.85 .90 .95];
                                 
                                  %
                                 
                                  % = 1:2
                                 
                                  % condenser_T = T1 (a);
                                 
                                  % (Net_profit (a)] = SteamPlant (T1,600 15 8,1 13 .85, .85);
                                 
                                  %结束
                                 
                                  %总(,,)= (Net_profit);
                                 
                                  %定义T1, T2, P1, P2范围:
                                 
                                  T1 = 40:45;
                                 
                                  T3 = 550:600;
                                 
                                  %值P1和P2
                                 
                                  P3 = 5:15;
                                 
                                  P4 = 1:15;
                                 
                                  %你也需要一个矩阵将你所有的结果
                                 
                                  %结果等于长度(T1) * (T2) *长度(P1) *长度(P2)。
                                 
                                  %定义一个计数器来存储你的结果:
                                 
                                  counter = 1;
                                 
                                  为i = 1:长度(T1)
                                 
                                  为j = 1:长度(T3)
                                 
                                  为k = 1:长度(P3)
                                 
                                  为m = 1:长度(P4)
                                 
                                  Net_profit(柜台)= SteamPlant (i, j, k、m);
                                 
                                  Input_data(柜台)=(计数器,i, j, k、m];
                                 
                                  counter = counter + 1;
                                 
                                  结束
                                 
                                  结束
                                 
                                  结束
                                 
                                  结束
                                 
                                  函数(表、Specific_Work Net_profit) = SteamPlant (T1, T3, P3, P4、Reheat1 Reheat2, Turbine1, Turbine2)
                                 
                                  %函数使用输入参数给周期属性
                                 
                                  state3 = steam_properties (“P”P3,“T”T3 + 273.15,“如果”);
                                 
                                  state4ideal = steam_properties (“P”P4,“年代”state3.s,“如果”);
                                 
                                  h4a = 0.15 * state3.h + Turbine1 * state4ideal.h;
                                 
                                  state1 = steam_properties (“T”T1 + 273.15,”[],“如果”);
                                 
                                  state4actual = steam_properties (“P”P4,“h”h4a,“如果”);
                                 
                                  state2ideal = steam_properties (“P”P3,“年代”state1.s_l,“如果”);
                                 
                                  h2a = ((state2ideal.h-state1.h_l) /炮)+ state1.h_l;
                                 
                                  state2actual = steam_properties (“P”P3,“h”h2a,“如果”);
                                 
                                  state5 = steam_properties (“P”P4,“T”T3 + 273.15,“如果”);
                                 
                                  state6ideal = steam_properties (“T”T1 + 273.15,“年代”state5.s,“如果”);
                                 
                                  h6a = 0.15 * state5.h + Turbine2 * state6ideal.h;
                                 
                                  state6actual = steam_properties (“P”state1.P,“h”h6a,“如果”);
                                 
                                  state2b = steam_properties (“P”state2ideal.P,”0,“如果”);
                                 
                                  %表向量
                                 
                                  州= {“状态1”;“理想状态2”;“实际状态2”;“状态3”;“理想状态4”;“实际状态4”;“国家5”;“国家6理想”;“国家实际6”};
                                 
                                  温度= [state1.T;state2ideal.T;state2actual.T;state3.T;state4ideal.T;state4actual.T;state5.T;state6ideal.T;state6actual.T];
                                 
                                  压力= [state1.P;state2ideal.P;state2actual.P;state3.P;state4ideal.P;state4actual.P;state5.P;state6ideal.P;state6actual.P];
                                 
                                  焓= [state1.h_l;state2ideal.h;state2actual.h;state3.h;state4ideal.h;state4actual.h;state5.h;state6ideal.h;state6actual.h];
                                 
                                  熵= [state1.s_l;state2ideal.s;state2actual.s;state3.s;state4ideal.s;state4actual.s;state5.s;state6ideal.s;state6actual.s];
                                 
                                  质量= [0;state2ideal.x;state2actual.x;state3.x;state4ideal.x;state4actual.x;state5.x;state6ideal.x;state6actual.x];
                                 
                                  %的表
                                 
                                  表=表(状态、温度、压力、焓、熵、质量);
                                 
                                  %计算工作
                                 
                                  Net_work = 20000
                                 
                                  Specific_Work_1to2 = h2a - state1.h_l
                                 
                                  Specific_HT_2to3 = state3.h-h2a
                                 
                                  Specific_Work_3to4 = state3。h - h4a
                                 
                                  Specific_HT_4to5 = state5。h - h4a
                                 
                                  Specific_Work_5to6 = state5。h - h6a
                                 
                                  MassFlow = Net_work / (Specific_Work_3to4 + Specific_Work_5to6-Specific_Work_1to2)
                                 
                                  Total_Work_1to2 = MassFlow * Specific_Work_1to2
                                 
                                  Total_HT_2to3 = MassFlow * Specific_HT_2to3
                                 
                                  Total_Work_3to4 = MassFlow * Specific_Work_3to4
                                 
                                  Total_HT_4to5 = MassFlow * Specific_HT_4to5
                                 
                                  Total_Work_5to6 = MassFlow * Specific_Work_5to6
                                 
                                  % t蒸汽圆顶
                                 
                                  Ts = linspace (20 + 273.15, 373 + 273.15, 500);
                                 
                                  为我= 1:1:长度(Ts)
                                 
                                  satt = steam_properties (“T”Ts (i),”[],“如果”);
                                 
                                  x1t (i) = satts.s_l;
                                 
                                  x2t (i) = satts.s_v;
                                 
                                  走(我)= Ts(我);
                                 
                                  结束
                                 
                                  Ts_pointx = [x1t(结束),x2t(结束)];
                                 
                                  Ts_pointy =[走(结束),走(结束)];
                                 
                                  % p - v蒸汽圆顶
                                 
                                  Pv = linspace (5 e - 3, 22, 10000);
                                 
                                  为我= 1:1:长度(Pv)
                                 
                                  satpv = steam_properties (“P”光伏(i),”[],“如果”);
                                 
                                  x1p (i) = satpv.v_l;
                                 
                                  x2p (i) = satpv.v_v;
                                 
                                  ypv (i) = Pv (i);
                                 
                                  结束
                                 
                                  Pv_pointx = [x1p(结束),x2p(结束)];
                                 
                                  Pv_pointy = [ypv(结束),ypv(结束)];
                                 
                                  %的图表
                                 
                                  图(1)
                                 
                                  xT = [state1。s_l state2ideal。年代,state2b。s_l state2b。s_v state3。年代,state4ideal。年代,state5。年代,state6ideal。年代,state1.s_l];
                                 
                                  欧美= [state1。T, state2ideal。T, state2b。T, state2b。T, state3。T, state4ideal。T, state5。T, state6ideal。T, state1.T]
                                 
                                  股上扬= [state1。s_l state2actual。年代,state2b。s_l state2b。s_v state3。年代,state4actual。年代,state5。年代,state6actual。年代,state1.s_l];
                                 
                                  (state1 yTa =。T, state2actual。T, state2b。T, state2b。T, state3。T, state4actual。T, state5。T, state6actual。T, state1.T];
                                 
                                  TsPlot =情节(xT,欧美,“b - o ')
                                 
                                  持有在
                                 
                                  情节(x1t,走,“r”)
                                 
                                  持有在
                                 
                                  情节(x2t,走,“r”)
                                 
                                  持有在
                                 
                                  情节(yTa股上扬,“米以下方式”)
                                 
                                  持有在
                                 
                                  情节(Ts_pointx Ts_pointy,“r”)
                                 
                                  包含(“熵(kJ /公斤* k)”)
                                 
                                  ylabel (“温度(k)”)
                                 
                                  图(2)
                                 
                                  xP = [state1。v_l state2ideal。v, state2b。v_l state2b。v_v state3。v, state4ideal。v, state5。v, state6ideal。v, state1.v_l];
                                 
                                  yP = [state1。P state2ideal.P state2b。P, state2b。P, state3。P, state4ideal。P, state5。P, state6ideal。P, state1.P];
                                 
                                  xPa = [state1。v_l state2actual。v, state2b。v_l state2b。v_v state3。v, state4actual。v, state5。v, state6actual。v, state1.v_l];
                                 
                                  (state1 yPa =。P state2actual.P state2b。P, state2b。P, state3。P, state4actual。P, state5。P, state6actual。P, state1.P];
                                 
                                  PvPlot =情节(xP, yP,“b - o ')
                                 
                                  持有在
                                 
                                  情节(xPa yPa,“米以下方式”)
                                 
                                  持有在
                                 
                                  情节(x1p ypv,“r”)
                                 
                                  持有在
                                 
                                  情节(x2p ypv,“r”)
                                 
                                  持有在
                                 
                                  情节(Pv_pointx Pv_pointy,“r”)
                                 
                                  集(gca),“Xscale”,“日志”);
                                 
                                  集(gca),“Yscale”,“日志”);
                                 
                                  包含(“具体体积(m ^ 3 /公斤))
                                 
                                  ylabel (“压力(MPa)”)
                                 
                                  %循环效率
                                 
                                  % % Ideal_efficieny = (1 - ((state3。h - state2ideal.h_l) + (state5) / (state3.h-state2ideal.h)
                                 
                                  Actual_efficiency = ((state3.h-state2actual.h) + (state5.h-state4actual.h)——(state6actual.h-state1.h_l)) / ((state3.h-state2actual.h) + (state5.h-state4actual.h))
                                 
                                  %天然气成本
                                 
                                  BTU = (Total_HT_2to3 + Total_HT_4to5) / 1.06
                                 
                                  mmBTU = BTU / (10 ^ 6)
                                 
                                  NG_cost_year = mmBTU * 3 * 86400 * 365
                                 
                                  %动力费用
                                 
                                  MW = ((Total_Work_5to6 + Total_Work_3to4) - Total_Work_1to2) / 1000
                                 
                                  Power_revenue_year = 25 MW * * 24 * 365
                                 
                                  Net_profit = Power_revenue_year - NG_cost_year
                                 
                                  %权衡曲线
                                 
                                  Trade_off = linspace (Reheat1 Reheat2 500);
                                 
                                  为我= 1:1:长度(Trade_off)
                                 
                                  Reheat_prop = steam_properties (“P”Trade_off(我),”[],“如果”);
                                 
                                  % x1p (i) = satpv.s_l;
                                 
                                  % x2p (i) = satpv.s_v;
                                 
                                  % ypv (i) = Pv (i);
                                 
                                  结束
                                 
                                  结束

沃尔特·罗伯森 2019年2月13日

我预测你试图在SteamPlant存储所有这一切。m。当你定义一个函数在脚本函数不能具有相同的名称作为脚本。

登录置评。

类别

工程机械工程热力学和热传递

找到更多的在热力学和热传递在帮助中心和文件交换

下载188bet金宝搏

MATLAB

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!