matlab Code 电极电势转换

clc
clear
close all
%transform between Hg/HgO---RHE---Ag/AgCl
RHE=-0.1:0.1:1.4;
PH1=14;
PH2=-log10(10.^(-14)/0.5);
PH3=0;

% E1_AgCl=RHE-0.059*PH1-0.1976;
% E1_HgO=RHE-0.059*PH1-0.098;
% E1_SCE=RHE-0.059*PH1-0.241;
% E1_HgCl2=RHE-0.059*PH2-0.245;
% E2_AgCl=RHE-0.059*PH2-0.1976;
% E2_HgO=RHE-0.059*PH2-0.098;
% E2_SCE=RHE-0.059*PH2-0.241;
% E2_HgCl2=RHE-0.059*PH2-0.245;

E3_AgCl=RHE-0.059*PH3-0.1976;
E3_HgO=RHE-0.059*PH3-0.098;
E3_SCE=RHE-0.059*PH3-0.241;
E3_HgCl2=RHE-0.059*PH3-0.245;
%E(RHE) = E Ag/AgCl + 0.059 pH + E0 Ag/AgCl
% Where E0Ag/AgCl = 0.1976 V at 25oC and E Ag/AgCl is your working potential.
% E (vs. RHE) = E (vs. Hg/Hg2Cl2) + 0.245 + 0.0591 V × pH
%E(RHE) = E(SCE) + 0.241 + 0.059 pH = E(Hg/HgO) +0.098 + 0.059 pH.
%chui xian
Sizes=size(RHE);
aa=-0.2*ones(Sizes(2),Sizes(1));
bb=-0.4*ones(Sizes(2),Sizes(1));
for ii=1:1:Sizes(2)
    line([-0.2,RHE(ii)],[E3_SCE(ii),E3_SCE(ii)]','linestyle','--', 'Color','r', 'LineWidth', 2)
    line([RHE(ii),RHE(ii)],[E3_SCE(ii),bb(ii)],'linestyle','--', 'Color','r', 'LineWidth',2)
end

% hold on
% plot(RHE,E3_AgCl,'.b-');
hold on
plot(RHE,E3_SCE,'.b-');
% hold on
% plot(RHE,E3_HgCl2,'.g-');
hold off
xlabel('RHE');
axis equal
axis([-0.2 1.6 -0.4 1.2])
set(gca,'ytick',-0.2:0.05:1.2)
E=[RHE;E3_AgCl;E3_SCE];
title('0.5M H_2SO_4     E3_S_C_E=RHE-0.059*PH-0.241')