MATLAB画双轴-倒立降水柱状图

% P=fix(rand(1,400)*500);   %产生所需数据

% Qobs=fix((rand(1,400)+10)*5);

% Qsim=Qobs-1;

function [NC,WBE]=Nash(Qobs,Qsim,P)

[m,n]=size(Qsim); 

if nargin==3  %如果输入参数为3个

    figure;  

    [AX,H1,H2]=plotyy(1:n,Qsim,1:n,P,'plot','bar');%画双轴，AX(1)左轴，AX（2）右轴，H为曲线本身

    hold on    

    set(AX(2),'YDir','reverse','Ylim',[0,max(P)*2],'YTick',[0:100:max(P)*2],'FontSize',12);  %设置右边轴为倒立

    set(AX(1),'YLim',[min(Qsim),1.3*max(Qsim)],'YTick',[min(Qsim):10:1.3*max(Qsim)],'Fontsize',12,'YColor','k');

    H3=plot(1:n,Qobs,'r-','LineWidth',1.1);

    %设置坐标轴的标题

    set(get(AX(1),'Xlabel'),'String','时间(d)');

    set(get(AX(1),'Ylabel'),'String','流量（m^3/s)');

    set(get(AX(2),'Ylabel'),'string','降雨量/mm');

    %设置figuer中线，柱的属性

    set(H2,'BarWidth',0.1);

    set(H1,'LineWidth',2,'Color','g');  %设置H1的曲线宽为4，颜色为g绿色,

    set(H3,'LineWidth',2,'Color','b')

    %设置背景色和图例

    set(gcf,'color','w');  %设置当前figure的背景颜色

    legend('模拟值','实测值','降水量');%设置图例

    WBE=sum(Qobs-Qsim)/sum(Qobs);

    NC=1-sum((Qobs-Qsim).^2)/sum((Qobs-mean(Qobs)).^2);

elseif nargin<2

     errordlg('less input！','Error','on');

else    

 WBE=sum(Qobs-Qsim)/sum(Qobs);

 NC=1-sum((Qobs-Qsim).^2)/sum((Qobs-mean(Qobs)).^2);
end