MATLAB 图中图

baqiao

有时为了在一个较大坐标轴上一个一个小图来说明局部方法的结果，可以调用axes函数实现这个要求。下面通过绘制 y=1/(t-3) 的曲线举例说明该函数的做法。程序如下：

%
clc;clear;

t=linspace(0,6,300);t1=linspace(2.8,3.2,300);

y=sin(1./(t-3));y1=sin(1./(t1-3));

+ L& P  B( o% W" [plot(t,y);axis('equal');

axes('Position',[0.18,0.62,0.28,0.25]); % 生成子图
plot(t1,y1); % 绘制局部曲线图     值得注意的是上一步生成局部坐标系，此时用plot就是在那个局部坐标系作图
) }+ p6 b2 o: B& \, C5 {: \
%但是如果在上一步生成坐标系之后，鼠标点了下外面的大坐标系，那么当前坐标系就切换到了外面的打坐标%系，而不是那个局部坐标系，这样plot（t1，y1）就画到外面的大坐标系中了
xlim([min(t1),max(t1)]); % 设置坐标轴范围
4 u3 K. n; T. p/ h+ M9 Pset(gcf,'Color','w');



生成简单的图中图，即在原图上选择放大区域，然后将该区域放大。但是这里不太“智能”，需要手动调所需放大的区域，手动箭头指向。


$ {( @# m, n# ]! Y! @
  E9 M9 d' x0 P1 J8 y另一种较为自动的方式是，编程实现放大区域的选择，以及放大的图像展示

4 f, S2 [6 C$ n$ ]$ qhttp://www.ilovematlab.cn/forum.php?mod=viewthread&tid=913&page=1

实现矩形区域的选择，并捕捉该区域的横纵标的范围。mmzoom创建一个缩放坐标轴。
0 ]# f8 k4 |9 B8 V! _2 a) d2 ~

! S: E4 f& k# N
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

把你要标注的曲线的句柄和字符串包含到legend中
2 p, X( s* W: Z* {2 c
3 j% o& e6 k/ bfigure(1);
5 g3 J( w* ?' |h2=axes('position',[0 0 1 1]); % 获得句柄
axis(h2);     %返回 坐标轴范围               为什么要加上这句  

2 t, d5 k5 o9 H; z8 N% {x2=0:pi/50:2*pi;
9 l- h& D3 l1 ]2 J6 `! ^& p0 vy2=sin(x2);
h3=plot(x2,y2,'b-');
h1=axes('position',[0.3 0.2 0.4 0.4]);
# z% y% j3 {0 n: ^4 Z  Gaxis(h1);
: `0 a: @/ \- k7 Fx1=0:pi/50:2*pi;
! r- l8 X. Z5 h) [  J8 by1=cos(x1);
h4=plot(x1,y1,'r-');
$ c( K8 i( N; u! S) Y9 Nhold on
h=[h3; h4];
' N& A# w) m, ostr=['大图中的曲线';'小图中的曲线'];
3 @: P2 y% p& klegend(h,str);
+ \7 B$ m- e4 @# y( s% U. o
+ a! y4 t, \9 E3 W2 b8 _; n9 e

- o, r4 p0 w; z! z: J8 \
( M/ c1 {5 H* `6 ^  R% ~$ C; K++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
9 ^; K3 N0 H0 U1 _/ X# a
clear
3 c# P' N3 o: lclc
close all
$ W1 |! S8 c# P7 F3 h2 {8 x7 qx=-2*pi:0.05*pi:2*pi;
y=sin(0.5*x).*cos(2*x);
/ `$ _$ j; S! b# t# V3 |% Ofh=figure;
6 [! l( g/ V2 ^2 \( r% Dah=axes('Parent',fh,'Layer','bottom','Position',[0.1 0.1 0.8 0.7],...
    'XLim',[-8 10],'YLim',[-1 1],'XColor','b','YColor','b',...
& T. X* k, y: {+ m' D    'Box','on','XTick',[],'YTick',[]);
9 E, W& P- ?' Jahx=axes('Parent',fh,'Layer','top','Position',[0.1 0.1+0.5*0.7 0.8 eps],...
    'XLim',[-8 10],'LineWidth',2.5,'XColor','m','XAxisLocation','top',...
    'YColor','w','YTick',[]);
5 b$ F0 S0 {. L6 P+ gahy=axes('Parent',fh,'Layer','top','Position',[0.1+0.8*8/18 0.1 eps 0.7],...
) U7 ]2 O- ~: m7 X% }$ R9 W    'YLim',[-1 1],'LineWidth',2.5,'YColor','r',...
# R2 K  K& N* t& L+ A4 v    'XColor','w','XTick',[]);
2 T& J2 k& x. ?3 Jline('XData',x,'YData',y,'Parent',ah,'Color','g')
* k0 U4 t7 v  k. b& S5 rahh=axes('Parent',fh,'Layer','top','Position',[0.68 0.12 0.2 0.2]);
9 O4 k6 ~: c; G$ t  y$ z  ^set(fh,'CurrentAxes',ahh);
lh=ezplot(@(x)(sin(0.5*x).*cos(2*x)),[0.35 0.55]);
. H- U9 _/ D6 w. z# v5 B+ ]( j: X* nset(get(ahh,'Title'),'String','zoomed x:[0.35 0.55]','Color','b');
  X/ _3 Q2 y3 [set(ahh,'XColor','b','YColor','b','Box','on','LineWidth',0.5,...
/ Q! s- z3 g0 G* k    'XTick',[0.45],'XTickLabel','0.45','YTick',[0.13],'YTickLabel','0.13',...
! q4 o1 E$ F5 Y; K9 q- w* V    'XGrid','on','YGrid','on')
set(get(ahh,'XLabel'),'String','')
: d' \$ ?  M; a. nset(get(ahh,'YLabel'),'String','')
& Q. h5 i+ t2 g6 P. d% `set(lh,'LineStyle','--','Color','r')
0 O* `8 [# T2 q. _; k, Pset(get(ah,'Title'),'String','Handle-Graphics','FontSize',25,'Color','b')
% x# S# ~0 D- q% G; T2 Z8 n6 g$ [8 }

$ d1 H7 b2 C$ F

CCxiaom

MATLAB 图中图