function V=pemfc(x)
% load simudata
% I=simudata(:,1);
V=zeros(200,1);
Copyright By PowCoder代写 加微信 powcoder
W=zeros(200,1);
load testdata
I=testdata(:,1);
U=testdata(:,2);
Pa=testdata(:,3);
Pc=testdata(:,4);
T=testdata(:,5);
for i=1:200
% load pemfcdata
% I=pemfcdata(:,2);
% U=pemfcdata(:,1);
% x1=-0.9512;
% x2=4.3e-5;
% x3=7.4e-5;
% x4=-1.87e-4;
% x5=3.8337e-6;
% x6=0.0467;
% x7=0.1855;
% r= 23;%为膜含水量
% Rc=1e-4; %欧姆阻抗中的contact resistances
% b=0.016;
% x=[x1,x2,x3,x4,x5,x6,x7,r,Rc,b]’;
r= x(8);%为膜含水量
Rc=x(9); %欧姆阻抗中的contact resistances
b=x(10);%为浓差极化的系数
RHa=1;%RHa是阳极中水蒸气的相对湿度
RHc=1;%RHc是阳极中水蒸气的相对湿度
% T=343.15;%设定电堆运行温度为70摄氏度，即343.15K
% Pa=3;%阳极气体的压力为1bar
% Pc=5;%阴极气体的压力为1bar
A=27;%活化面积为27cm-2
J=I(i)/A;%J为电流密度，单位是A·cm-2
Jmax=0.86;%Jmax是最大电流密度，单位是A·cm-2
L=1.27e-2;%L是质子交换膜的厚度，单位为cm
P_H20=10^(2.95e-2*(T(i)-273.15)-9.18e-5*((T(i)-273.15)^2)+1.44e-7*((T(i)-273.15)^3)-2.18);%P_H2O是水蒸汽分压
u=[RHa,P_H20,T(i),0.9869*Pa(i),0.9869*Pc(i),RHc]’;
P_H2=0.5*u(1)*u(2)*((exp(1.635*J/(u(3)^1.334))*u(1)*u(2)/u(4)).^(-1)-1);%计算氢气分压
P_O2=u(6)*u(2)*((exp(4.192*J/(u(3)^1.334))*u(6)*u(2)/u(5)).^(-1)-1);%计算氧气分压
E=1.229+(log(P_H2)+0.5*log(P_O2))*u(3)*(4.308e-5)-(u(3)-298.15)*(8.5e-4);
CH2=P_H2/(1.09e6*exp(77/u(3)));%计算氢气浓度
CO2=P_O2/(5.08e6*exp(-498/u(3)));%计算氧气浓度
Vact=(x(1)+0.00286*u(3)+0.0002*u(3)*log(A)+x(2)*u(3)*log(CH2)+x(3)*u(3)*log(CO2)+x(4)*u(3)*log(I(i)))*(-1);
Vohm=((J.*0.03+1+(u(3)/303)^2*(J.^2.5)*0.062)*181.6./((r-0.634-J.*3)*exp(4.18*(1-303/u(3))))*L/A+Rc).*I(i);
Vcon=-b*log(1-J./Jmax)+x(5)*(exp(x(6)*I(i)-x(7)/u(5)))/u(5);
Vcell=E-Vact-Vohm-Vcon;%单个燃料电池的输出电压
Vout=24*Vcell;
W(i)=Vout;
K(i)=(Vout-U(i));
V=(sum(K.^2));
程序代写 CS代考 加微信: powcoder QQ: 1823890830 Email: powcoder@163.com