matlab 画图双坐标画图问题

未解决问题：1、若解决双坐标与子图并和问题，右端曲线问题

                        2、相反 。。。

因为plotyy画图曲线右端没能到坐标轴，以下程序实现，不过subplot用不了

解决右端曲线问题：

%import data
x1=linspace(1,226235,226235);
bx=importdata('bx.txt');
%by=importdata('by.txt');
%bz=importdata('bz.txt');

x2=linspace(1,222017,222017);
vpx=importdata('vx.txt');
%vpy=importdata('vy.txt');
%vpz=importdata('vz.txt');

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% fit

figure(1);
subplot(2,1,1);
subplot(2,1,1,'replace'); 
ax=gca;
subplot(2,1,1,ax);
%ax = subplot();
subplot(ax);

clf reset,h_ap=axes('Position',[0.13,0.13,0.7,0.75]);
set(gca,'xlim',[0 226235]);
set(gca,'xTick',[0:22623.5:226235]);
%set(gca,'xtick',[],'xticklabel',[])
set(gca,'yLim',[-15 15]);
set(gca,'yTick',[-15:5:15]);
ylabel('Bx(nT)');
line(x1,bx,'Color','b');
 
h_at=axes('Position',get(h_ap,'Position')); 
set(h_at,'Color','none','Xcolor','k','Ycolor','r'); 
set(h_at,'Xaxislocation','top')；
set(h_at,'Yaxislocation','right','Ydir','rev') %<15>  % 以上四句语句实现上X坐标，右Y坐标

set(gca,'xlim',[0 222017]);
set(gca,'xTick',[0:22201.7:222017]);
set(gca,'xtick',[],'xticklabel',[]);
set(gca,'yLim',[-900 -300]);
set(gca,'yTick',[-900:100:-300]);
set(gca,'YDir','normal');
line(x2,vpx,'Color','r','Parent',h_at) 
ylabel('Vpx(km/s)','Rotation',270);

%subplot(2,1,1)

解决双坐标与子图问题

x=linspace(1,226235,226235);
y=importdata('bx.txt');
x1=linspace(1,222017,222017);
y1=importdata('vx.txt');

%hold on;
figure(1);
subplot(2,1,1);
[AX,H1,H2] =plotyy(x,y,x1,y1);% 获取坐标轴、图像句柄
%set(get(AX(1),'ylabel'),'string', '峰值透过率','fontsize',16);
%set(get(AX(2),'ylabel'),'string', '半高宽/nm','fontsize',16);
set(gca,'xtick',[],'xticklabel',[]);
%xlabel('倾斜角度/°','fontsize',16);
set(H1,'Linestyle','-');
set(H2,'Linestyle','-');
set(gcf,'color','white')
set(gca,'linewidth',1.5) %设置边框宽度  

x=linspace(1,226235,226235);
by=importdata('by.txt');
x1=linspace(1,222017,222017);
vpy=importdata('vy.txt');

hold on
figure(1);
subplot(2,1,2);
[AX,H1,H2] =plotyy(x,by,x1,vpy);% 获取坐标轴、图像句柄
%set(get(AX(1),'ylabel'),'string', '峰值透过率','fontsize',16);
%set(get(AX(2),'ylabel'),'string', '半高宽/nm','fontsize',16);
xlabel('倾斜角度/°','fontsize',4);
set(H1,'Linestyle','-');
set(H2,'Linestyle','-');
set(gcf,'color','white')
set(gca,'linewidth',0.5) %设置边框宽度

参考：

哲的博文：http://bilu.blog.163.com/blog/static/176178129201182943032724/

baidujingying: https://jingyan.baidu.com/article/37bce2be63d4bf1002f3a23e.html