解决使用牛顿拉富生method-Continuosly运行

3视图(30天)

显示旧的评论

一边抚摸Chandrasekar 2017年8月31日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/354737-solve-using-newton-raphson-method-continuosly-running

评论道: 克丽丝Fedorenko 2017年9月7日

                         函数电解()
                        
                         关闭所有
                        
                         clc
                        
                         清晰的所有
                        
                         % %定义常量
                        
                         stoic_coeff_H2 = 2;
                        
                         stoic_coeff_O2 = 1;
                        
                         F = 96485.3365;% amppere-sencond /摩尔库仑/秒
                        
                         R = 8.314;%通用气体常数,j.mol - 1. - k - 1
                        
                         Power_fluc = 2 * 10 ^ 6;
                        
                         P = 100 * 10 ^ 6;% Pa
                        
                         P_H = 1.013 * 10 ^ 6;% Pa
                        
                         nfar_eff = 1;%法拉第效率,通常假设为1
                        
                         I_0 = 0.005;%安培/平方米
                        
                         Area_elec = 0.05;% m2,通常提供了术语当前/ apecific区域,因此使用这个术语
                        
                         α= 0.5;%的co-rffici ent分裂,通常0.5(假设等量离子的阳极和阴极)
                        
                         delta_thickness = 0.0178 * 10 ^ (2);%厚度的膜电解电池,土地
                        
                         λ= 15;% PEM膜中的含水量
                        
                         m = 1000;%的摩尔浓度elctrolyte AE,摩尔/ m3
                        
                         sigma_PEM = 0.14;
                        
                         T = 150;
                        
                         T = 273.15 + T;

% %牛顿拉富生方法i = 1:尺寸(Power_fluc 2);Y = @(我)(I * (1.229 - (8.5 * 10 ^ (4) * (T - 298)) + (R * T *日志(1.5 * (P-P_H) / P_H)) / ((2 * F)) + (R * T * asin (I / (2 * I_0)) /(α2 * * F)) + (delta_thickness。*我/ sigma_PEM))) f_prime = @ (I)(-((即* ((R * T / (sqrt ((4 * (I_0 ^ 2)) - (I ^ 2)) * 2 *α* F)) + (delta_thickness / sigma_PEM))) + (1.229 - 8.5 * 10 ^ (4) * (T - 298)) + (R * T *日志(1.5 * (P-P_H) / P_H) / (2 * F)) + (R * T * asin (I / (2 * I_0)) /(α2 * * F)) + (delta_thickness。*我/ sigma_PEM))));% I = 0(大小(Power_fluc, 2));iter = 0;% sin_value = (I / (2 * I_0));%如果sin_value > -π/ 2%如果sin_value <π/如果我2% < (2 * I_0), abs (Y (I)) > (10 ^ 3) iter = iter + 1;x1 = I - (Y(我)/ (f_prime(我)));我= x1;Y(我)结束结束% % % iter ndot_H2 = nfar_eff结束。*我(我)* Area_elec / (stoic_coeff_H2 * F) %的氢气流ndot_O2 = nfar_eff。*我(我)* Area_elec / (stoic_coeff_O2 * F) %的氧气流图绘制(我Power_fluc)图绘制(我ndot_H2)图散射(Power_fluc、ndot_H2 Power_fluc, ndot_O2)结束