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')
