Python _Mix*9

  • 1. 函数
    • 函数是对功能的封装
  • 语法:
    • def 函数名(形参列表):
      • 函数体(代码块) 代码块中有可能包含return
    • 调用:
      • 函数名(实参列表)
  • 1 def mix(a,b):   #def 函数名(a和b代表的是形参)
    2     return a if a > b else b    #return是返回结果
    3 mix(1,2)        #这里就是函数的调用
    4 print(mix(3,2))

     

  • 2. 返回值*(给函数的使用者一个访问结果)
    • return : 在函数执行的时候. 如果遇到return. 直接返回,
      • 1.如果函数什么都不写, 不写return, 没有返回值. 得到的是None
      • 2.在函数中间或者末尾写return, 返回的是None
  • def mix(a,b):  #def 函数名(a和b代表的是形参)
        return  #在函数中间或者末尾写return,返回的结果是None
    
    mix(1,2)        #这里就是函数的调用  括号里的1和2就是实参
    print(mix(3,2))  
    #打印的结果是None

     

      • 3.在函数中写return 值. 返回一个值.
      • 4.在函数中可以返回多个返回值, return 值1, 值2, 值3...., 接收的是元组
  • 3. 参数(函数执行的时候给形参传递信息)
    • *形参:函数声明的位置的变量
    • *实参:函数调用的时候给的具体的值
    • 传参:把实参交给形参的过程
      • 1. 实参:
        • 1. 位置参数, 按照形参的参数位置, 给形参传值
        • 2. 关键字参数, 按照形参的名字给形参传值
        • 3. 混合参数. 即用位置参数, 也用关键参数
      • 顺序是:先位置后关键字
      • 2. 形参:
        • 1. 位置参数
        • 2. 默认值参数 先位置后默认值

转载于:https://www.cnblogs.com/mixw/p/9445710.html

``` function fo = gaCAES1(x) % Caution: 杈撳叆鐨勫帇鍔涙槸kPa锛屽昂瀵告槸mm %% 绯荤粺鎬诲拰 %鍒濆鍖栧弬鏁? ambient_pressure=101.32; ambient_temperature=293.15; P_start=x(6); P_end=7500; compression_ratio=(P_end/ambient_pressure)^(1/3); [h_ambient s_ambient]=refpropm(&#39;HS&#39;,&#39;T&#39;,ambient_temperature,&#39;P&#39;,ambient_pressure,&#39;air.mix&#39;); %#ok<*NCOMMA> delta_m=370000; % vol=4527.23;% m^3 D_end=refpropm(&#39;D&#39;,&#39;T&#39;,ambient_temperature,&#39;P&#39;,P_end,&#39;air.mix&#39;); % m_end=D_end*vol; D_start=refpropm(&#39;D&#39;,&#39;T&#39;,ambient_temperature,&#39;P&#39;,P_start,&#39;air.mix&#39;); % m_start=D_start*vol; vol=delta_m/(D_end-D_start); %% 涓?绾у帇缂╂満 COMP1 %浣庡帇鍘嬬缉鏈? adiabatic_efficiency = 0.85; % 缁濈儹鏁堢巼 % compression_ratio = 7.87; % 鍘嬬缉姣?(宸茶鎬讳綋鎺у埗) m_rate=20; % 鍏呰兘娈电殑娴侀噺 s_low_comp_in=refpropm(&#39;S&#39;,&#39;T&#39;,ambient_temperature,&#39;P&#39;,ambient_pressure,&#39;air.mix&#39;); inlet_entropy=s_low_comp_in;% 杈撳叆鍏ュ彛鐔靛?硷紝鍗曚綅锛欽/(kg*K) P_low_comp_in=ambient_pressure;% 杈撳叆鍏ュ彛鍘嬪己锛屽崟浣嶏細kPa inlet_pressure=P_low_comp_in; outlet_pressure=P_low_comp_in*compression_ratio; % 璁$畻鍑哄彛鐒撳?? h2_adiabatic = refpropm(&#39;H&#39;,&#39;P&#39;,outlet_pressure,&#39;S&#39;,inlet_entropy,&#39;air.mix&#39;); % 缁濈儹杩囩▼涓嬬殑鍑哄彛鐒撳?? h1 = refpropm(&#39;H&#39;,&#39;P&#39;,inlet_pressure,&#39;S&#39;,inlet_entropy,&#39;air.mix&#39;); % 鍏ュ彛鐒撳?? h2_actual = h1 + 1/adiabatic_efficiency*(h2_adiabatic - h1); % 瀹為檯鍑哄彛鐒撳?硷紝鑰冭檻缁濈儹鏁堢巼 % 璁$畻鍘嬬缉鏈哄姛鑰? h_low_comp_work=h2_actual - h1; low_comp_work = h_low_comp_work*m_rate; % 鍘嬬缉鏈烘墍鍋氱殑鍔熺巼 % 璁$畻鍑哄彛娓╁害 low_comp_Temp_out=refpropm(&#39;T&#39;,&#39;H&#39;,h2_actual,&#39;P&#39;,outlet_pressure,&#39;air.mix&#39;); s_low_comp_out=refpropm(&#39;S&#39;,&#39;T&#39;,low_comp_Temp_out,&#39;P&#39;,outlet_pressure,&#39;air.mix&#39;); % fprintf(&#39;浣庡帇鍘嬬缉鏈虹殑杈撳叆鍔熶负%.0fkW\n&#39;,low_comp_work/1000); ex_low_comp_in=(h1-h_ambient)-ambient_temperature*(s_low_comp_in-s_ambient); ex_low_comp_out=(h2_actual-h_ambient)-ambient_temperature*(s_low_comp_out-s_ambient); ix_low_comp=ex_low_comp_in+low_comp_work/m_rate-ex_low_comp_out; %% 浣庡帇鎹㈢儹鍣? mc_low = 0.4*m_rate; % 鍐锋祦浣撹川閲忔祦閲忥紝鍗曚綅 kg/s T_hi = low_comp_Temp_out; % 鐑祦浣撳叆鍙f俯搴︼紝鍗曚綅 K cp_h = refpropm(&#39;C&#39;,&#39;T&#39;,T_hi,&#39;P&#39;,outlet_pressure,&#39;air.mix&#39;); % 鐑祦浣撳畾鍘嬫瘮鐑锛屽崟浣? J/(kg.K) while 1 %浣庡帇鎹㈢儹鍣? cond=0; % 鍒ゆ柇椤洪?嗘祦 cp_c = 2100; % 鍐锋祦浣撳畾鍘嬫瘮鐑锛屽崟浣? J/(kg.K) T_ci = 293.15; % 鍐锋祦浣撳叆鍙f俯搴︼紝鍗曚綅 K mh = m_rate; % 鐑祦浣撹川閲忔祦閲忥紝鍗曚綅 kg/s pp=5; % 鑺傜偣娓╁樊 P=outlet_pressure; omega=0; % 婕忔皵鐜? % 瑙f柟绋嬪緱鍒板喎娴佷綋鍑哄彛娓╁害 syms T_co T_ho %#ok<*NASGU> % 瀹氫箟绗﹀彿鍙橀噺 T_co,T_ho if cond==0 T_ho=T_ci+pp; %閫嗘祦 eqn=mc_low*cp_c*(T_ci-T_co)-mh*cp_h*(T_ho-T_hi)== 0; T_co_sol=solve(eqn,T_co); T_co_low=double(T_co_sol); P_out=P*(1-omega); else T_ho=T_co+pp; %椤烘祦 H_s=refpropm(&#39;H&#39;,&#39;T&#39;,T_ho,&#39;P&#39;,P,&#39;air.mix&#39;); H_a=refpropm(&#39;H&#39;,&#39;T&#39;,T_ho,&#39;P&#39;,P*(1-omega),&#39;air.mix&#39;); de_H=H_s-H_a; eqn = mc_low * cp_c * (T_ci - T_co) - mh * cp_h * (T_ho - T_hi) + de_H == 0; % 鐑钩琛℃柟绋? T_co_sol = solve(eqn, T_co); % 姹傝В鏂圭▼ T_co_low = double(T_co_sol); % 灏嗚В杞崲涓烘暟鍊? P_out=P*(1-omega); end if T_co_low<=T_hi-15 break else mc_low=mc_low+0.05; end end Q_ch1=mh*cp_h*(T_hi-T_ho); T_mean_ch1=(T_hi-T_ho+T_ci-T_co_low)/log((T_hi-T_co_low)/(T_ho-T_ci)); A_ch1=Q_ch1/0.8/1100/T_mean_ch1; h_air_charging_low_in=refpropm(&#39;H&#39;,&#39;T&#39;,T_hi,&#39;P&#39;,P,&#39;air.mix&#39;); h_air_charging_low_out=refpropm(&#39;H&#39;,&#39;T&#39;,T_ho,&#39;P&#39;,P_out,&#39;air.mix&#39;); s_air_charging_low_in=refpropm(&#39;S&#39;,&#39;T&#39;,T_hi,&#39;P&#39;,P,&#39;air.mix&#39;); s_air_charging_low_out=refpropm(&#39;S&#39;,&#39;T&#39;,T_ho,&#39;P&#39;,P_out,&#39;air.mix&#39;); % fprintf(&#39;浣庡帇鎹㈢儹鍣ㄧ┖姘旂殑鍑哄彛娓╁害涓?%fK\n&#39;,T_ho);```显示完整的修改后的代码
最新发布
03-26
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值