MATLAB 图中图

baqiao

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

2 j9 R2 s6 J; w, h5 q4 `" C%
2 t, q9 G9 C4 B5 Qclc;clear;
6 |0 w+ q4 r/ M3 H1 N" Q: o) x
( J3 z' w. k: zt=linspace(0,6,300);t1=linspace(2.8,3.2,300);
3 S/ J+ N8 D3 Z' W: J! {
y=sin(1./(t-3));y1=sin(1./(t1-3));
/ M5 C  H/ G3 g& p3 F' a; }
% {5 s2 H# f$ x# {; b1 ]. i0 Qplot(t,y);axis('equal');
! t+ A- R. G9 K4 w& ]
axes('Position',[0.18,0.62,0.28,0.25]); % 生成子图
0 U$ p: z% f: ~) ]9 G' Kplot(t1,y1); % 绘制局部曲线图     值得注意的是上一步生成局部坐标系，此时用plot就是在那个局部坐标系作图

/ l- {- u7 e) a+ l5 H2 b5 K/ M6 P3 n%但是如果在上一步生成坐标系之后，鼠标点了下外面的大坐标系，那么当前坐标系就切换到了外面的打坐标%系，而不是那个局部坐标系，这样plot（t1，y1）就画到外面的大坐标系中了
2 O' q" S1 R- U' e/ u: exlim([min(t1),max(t1)]); % 设置坐标轴范围
6 b" I7 g0 ?0 \set(gcf,'Color','w');

* Z3 t4 Z8 z9 g, a. M: B
5 ^3 \% [4 K0 {1 u5 ]7 t
& {' M* A. Y- f( b( E& K+ a生成简单的图中图，即在原图上选择放大区域，然后将该区域放大。但是这里不太“智能”，需要手动调所需放大的区域，手动箭头指向。
5 W* K: m0 h2 E/ y* w% B3 Y- L
: l- O# L: `! [. m' J* ]  u$ s
$ s+ n0 A, H- I
/ M7 N" e  d. d( R8 Q. `1 s另一种较为自动的方式是，编程实现放大区域的选择，以及放大的图像展示

/ ?4 ~! t( L1 g3 j% E' yhttp://www.ilovematlab.cn/forum.php?mod=viewthread&tid=913&page=1

实现矩形区域的选择，并捕捉该区域的横纵标的范围。mmzoom创建一个缩放坐标轴。

2 {' A  Y/ Z& G, e4 x

6 d4 E' Y( f; s4 |++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

把你要标注的曲线的句柄和字符串包含到legend中

7 z( v  H& ^% mfigure(1);
h2=axes('position',[0 0 1 1]); % 获得句柄
! t8 N. S! d! V, A1 o, m! R& caxis(h2);     %返回 坐标轴范围               为什么要加上这句  

x2=0:pi/50:2*pi;
  s( S) d. i# ?  r9 Zy2=sin(x2);
h3=plot(x2,y2,'b-');
h1=axes('position',[0.3 0.2 0.4 0.4]);
9 O( K. W4 V. O$ W/ Kaxis(h1);
) l8 j4 J' Z1 @% A  i2 `9 Ox1=0:pi/50:2*pi;
y1=cos(x1);
h4=plot(x1,y1,'r-');
hold on
6 v0 R) g' t1 _h=[h3; h4];
; h3 ]# {" `. W& r* Pstr=['大图中的曲线';'小图中的曲线'];
legend(h,str);


6 Q0 L  k& C# {

+ h8 z7 x- M& z% S# p; i' q++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

$ p4 H1 {5 n; R4 sclear
' n: c6 {, l3 N7 B! mclc
/ D( r, Y6 q% T+ M) a5 H. \close all
, x, d1 T: Z, z( Vx=-2*pi:0.05*pi:2*pi;
y=sin(0.5*x).*cos(2*x);
fh=figure;
ah=axes('Parent',fh,'Layer','bottom','Position',[0.1 0.1 0.8 0.7],...
- x" c( q5 F6 T  H! g0 m    'XLim',[-8 10],'YLim',[-1 1],'XColor','b','YColor','b',...
    'Box','on','XTick',[],'YTick',[]);
3 l/ T9 v7 d* h- [5 iahx=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',[]);
6 \, E4 y9 C% N! L* ^ahy=axes('Parent',fh,'Layer','top','Position',[0.1+0.8*8/18 0.1 eps 0.7],...
    'YLim',[-1 1],'LineWidth',2.5,'YColor','r',...
0 n  D. S, P9 ?( B, \    'XColor','w','XTick',[]);
line('XData',x,'YData',y,'Parent',ah,'Color','g')
ahh=axes('Parent',fh,'Layer','top','Position',[0.68 0.12 0.2 0.2]);
  A- l5 C% b7 p/ C3 ~9 ~; k- Hset(fh,'CurrentAxes',ahh);
+ G9 s: X$ W* Glh=ezplot(@(x)(sin(0.5*x).*cos(2*x)),[0.35 0.55]);
+ n8 C; _' ~( }. h* p, S0 k4 Z5 ]$ Bset(get(ahh,'Title'),'String','zoomed x:[0.35 0.55]','Color','b');
set(ahh,'XColor','b','YColor','b','Box','on','LineWidth',0.5,...
! W" P+ B  X* S% {) p, D& H    'XTick',[0.45],'XTickLabel','0.45','YTick',[0.13],'YTickLabel','0.13',...
    'XGrid','on','YGrid','on')
set(get(ahh,'XLabel'),'String','')
, |4 Y0 R; G' s% pset(get(ahh,'YLabel'),'String','')
set(lh,'LineStyle','--','Color','r')
set(get(ah,'Title'),'String','Handle-Graphics','FontSize',25,'Color','b')
- p3 J* v9 S/ z0 G3 y6 [! r  }
; r+ I$ n; U1 L7 C, {# ]5 V
5 C* {$ ^! n6 S7 l# V2 P/ u
6 I6 Z% q" R8 F8 }- f

CCxiaom

MATLAB 图中图