MATLAB 图中图

baqiao

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

%
3 g8 F, c( q, nclc;clear;

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

6 C( S2 U! ]# h% e( @- Ky=sin(1./(t-3));y1=sin(1./(t1-3));

4 A, {) M( y$ s2 ^# ]# L. U# Oplot(t,y);axis('equal');
+ }3 k* A8 ~: Q7 k0 N) `2 K
- N7 B/ l! l/ W( e; \# F0 Z* Taxes('Position',[0.18,0.62,0.28,0.25]); % 生成子图
2 R8 a# R) ], V  z# \plot(t1,y1); % 绘制局部曲线图     值得注意的是上一步生成局部坐标系，此时用plot就是在那个局部坐标系作图

%但是如果在上一步生成坐标系之后，鼠标点了下外面的大坐标系，那么当前坐标系就切换到了外面的打坐标%系，而不是那个局部坐标系，这样plot（t1，y1）就画到外面的大坐标系中了
/ U$ P: x1 Q1 V  O7 K' ixlim([min(t1),max(t1)]); % 设置坐标轴范围
! A. k0 O, X9 o/ w, S+ ~set(gcf,'Color','w');
; _! t7 @% t# H" _' ~! n
8 \, |$ v% F- Q; O& i- e, P( ~( P7 k
3 V% ]5 ^4 L6 n
生成简单的图中图，即在原图上选择放大区域，然后将该区域放大。但是这里不太“智能”，需要手动调所需放大的区域，手动箭头指向。

- }5 _) q/ \& `( ^% H. f2 s
+ F7 y; F$ e4 V, ~/ g$ x
& E% Q) Q0 K9 i) c4 @2 C2 C另一种较为自动的方式是，编程实现放大区域的选择，以及放大的图像展示
4 O" B6 z5 ~7 W7 \3 y  O( j
9 S2 |3 r) I' F5 Shttp://www.ilovematlab.cn/forum.php?mod=viewthread&tid=913&page=1
" Y6 m. S* m; T3 l! ?" S
实现矩形区域的选择，并捕捉该区域的横纵标的范围。mmzoom创建一个缩放坐标轴。



++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
& ]2 b) k: _, w7 u# q$ e0 k2 q
把你要标注的曲线的句柄和字符串包含到legend中
1 }  N* |1 y3 r# c: ^
; {: ?8 r% R9 O1 vfigure(1);
5 [+ `& ~8 W- U. r6 O; u! Ph2=axes('position',[0 0 1 1]); % 获得句柄
) j/ l* s% l# B0 }axis(h2);     %返回 坐标轴范围               为什么要加上这句  
. D; _+ v# r& N: u1 G2 @- b6 M, C
x2=0:pi/50:2*pi;
. G1 X; r6 [$ O" F# X! }8 vy2=sin(x2);
h3=plot(x2,y2,'b-');
h1=axes('position',[0.3 0.2 0.4 0.4]);
2 Y; l1 L8 h* maxis(h1);
x1=0:pi/50:2*pi;
. Z. i0 {( }7 ~) F, u0 g6 zy1=cos(x1);
h4=plot(x1,y1,'r-');
* z1 m- I2 N. c5 Dhold on
h=[h3; h4];
; ]% ?0 M/ l1 t$ k4 S5 Vstr=['大图中的曲线';'小图中的曲线'];
" W" r0 A( b- ]. \3 klegend(h,str);
& J5 P' o* ^. f; f1 n, V% M. m: o6 P3 J

: l/ E  E" i3 A0 J( K6 ?' n' c
3 w, x5 E$ l9 f
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

clear
clc
- d) H7 G8 A6 Aclose all
6 p8 ]: B9 v$ _8 J+ V) ux=-2*pi:0.05*pi:2*pi;
, u- s% M" u( h6 ?y=sin(0.5*x).*cos(2*x);
- i0 \5 g! I' ]fh=figure;
ah=axes('Parent',fh,'Layer','bottom','Position',[0.1 0.1 0.8 0.7],...
    'XLim',[-8 10],'YLim',[-1 1],'XColor','b','YColor','b',...
    'Box','on','XTick',[],'YTick',[]);
$ I0 m3 |8 m3 o4 Wahx=axes('Parent',fh,'Layer','top','Position',[0.1 0.1+0.5*0.7 0.8 eps],...
* N$ j3 p, f; D" k    'XLim',[-8 10],'LineWidth',2.5,'XColor','m','XAxisLocation','top',...
0 ~- ?/ F" Z/ E! @, h; X    'YColor','w','YTick',[]);
! ^& i6 l& v& o4 N. z! Wahy=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',...
    'XColor','w','XTick',[]);
line('XData',x,'YData',y,'Parent',ah,'Color','g')
8 X% N. O0 f8 w: y/ j/ |, X" C; Uahh=axes('Parent',fh,'Layer','top','Position',[0.68 0.12 0.2 0.2]);
: ~& y" J0 d9 J5 ?$ jset(fh,'CurrentAxes',ahh);
lh=ezplot(@(x)(sin(0.5*x).*cos(2*x)),[0.35 0.55]);
set(get(ahh,'Title'),'String','zoomed x:[0.35 0.55]','Color','b');
set(ahh,'XColor','b','YColor','b','Box','on','LineWidth',0.5,...
2 S; H/ P4 V# {0 Z( {3 J* V/ C    'XTick',[0.45],'XTickLabel','0.45','YTick',[0.13],'YTickLabel','0.13',...
. S+ W- k4 R! p  Y    'XGrid','on','YGrid','on')
set(get(ahh,'XLabel'),'String','')
set(get(ahh,'YLabel'),'String','')
, L- R, t- X6 T& g8 x6 w# Lset(lh,'LineStyle','--','Color','r')
set(get(ah,'Title'),'String','Handle-Graphics','FontSize',25,'Color','b')

2 {4 M4 n6 c  s; l# x

7 J5 z$ ]. d* ~6 `# R; {
3 g! x6 K6 v# d1 n3 c5 s

CCxiaom

MATLAB 图中图