MATLAB 图中图

baqiao

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

%
+ k. L' C$ Z( H' Pclc;clear;

- X4 H' L/ z1 H3 X0 ^t=linspace(0,6,300);t1=linspace(2.8,3.2,300);
- c$ W5 ^3 z3 k: a5 F3 B
& F# P- v# W' @, j3 I1 zy=sin(1./(t-3));y1=sin(1./(t1-3));
+ w. w& c+ f& B7 J& j* T
. ^5 o; e- ^$ ]! b, J# [, S1 [' E! Zplot(t,y);axis('equal');

8 S1 I" m: [* }) zaxes('Position',[0.18,0.62,0.28,0.25]); % 生成子图
& L3 P! H7 B5 Cplot(t1,y1); % 绘制局部曲线图     值得注意的是上一步生成局部坐标系，此时用plot就是在那个局部坐标系作图

# L; e0 e% R2 _) `9 s/ n0 [%但是如果在上一步生成坐标系之后，鼠标点了下外面的大坐标系，那么当前坐标系就切换到了外面的打坐标%系，而不是那个局部坐标系，这样plot（t1，y1）就画到外面的大坐标系中了
xlim([min(t1),max(t1)]); % 设置坐标轴范围
5 P" T" {9 \3 T) oset(gcf,'Color','w');
9 W2 o  E: N* k; ]" U& `

3 ?; M+ ^% m/ q: I1 j( U
. h  d: s6 j7 j9 i' Z# z生成简单的图中图，即在原图上选择放大区域，然后将该区域放大。但是这里不太“智能”，需要手动调所需放大的区域，手动箭头指向。

' D# u7 `' N8 E6 `; ~: \
$ l  K0 g" m) M
另一种较为自动的方式是，编程实现放大区域的选择，以及放大的图像展示

5 ?3 ]# b- l/ g( ~4 S1 {http://www.ilovematlab.cn/forum.php?mod=viewthread&tid=913&page=1

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

0 A$ w+ e. z4 G

6 U+ o6 \# K% J4 A++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
4 F6 o( C' `# o  t7 d% R1 p
把你要标注的曲线的句柄和字符串包含到legend中
3 X0 v5 ?4 F- X+ |& W
figure(1);
h2=axes('position',[0 0 1 1]); % 获得句柄
axis(h2);     %返回 坐标轴范围               为什么要加上这句  

x2=0:pi/50:2*pi;
2 q- ^3 X0 f# my2=sin(x2);
h3=plot(x2,y2,'b-');
h1=axes('position',[0.3 0.2 0.4 0.4]);
axis(h1);
x1=0:pi/50:2*pi;
1 R8 V" D, K7 D; f3 t# Iy1=cos(x1);
+ C4 t% s! q- X) ~% }( k5 O% Ih4=plot(x1,y1,'r-');
4 U  H: k  S0 Ghold on
h=[h3; h4];
str=['大图中的曲线';'小图中的曲线'];
legend(h,str);




++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
4 |4 C- }8 t5 i
clear
) {5 {9 B& g0 D/ i: Gclc
close all
, n* c/ K( ]8 }" Z+ \; C5 P) Vx=-2*pi:0.05*pi:2*pi;
y=sin(0.5*x).*cos(2*x);
# P+ n) l- E$ Z$ }5 r$ Y+ Yfh=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',[]);
ahx=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',...
6 ^: H# g% s  t  b    'YColor','w','YTick',[]);
ahy=axes('Parent',fh,'Layer','top','Position',[0.1+0.8*8/18 0.1 eps 0.7],...
+ R% Y6 l9 K. l* [    'YLim',[-1 1],'LineWidth',2.5,'YColor','r',...
    'XColor','w','XTick',[]);
line('XData',x,'YData',y,'Parent',ah,'Color','g')
, W3 J9 e! _: ~ahh=axes('Parent',fh,'Layer','top','Position',[0.68 0.12 0.2 0.2]);
set(fh,'CurrentAxes',ahh);
. @$ ?! d; M/ B7 d- N: E0 mlh=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');
" ]* \! f; ]* I: E& r1 N  ?set(ahh,'XColor','b','YColor','b','Box','on','LineWidth',0.5,...
    'XTick',[0.45],'XTickLabel','0.45','YTick',[0.13],'YTickLabel','0.13',...
    'XGrid','on','YGrid','on')
set(get(ahh,'XLabel'),'String','')
; x2 |; L# ~) J1 M# V* ]set(get(ahh,'YLabel'),'String','')
set(lh,'LineStyle','--','Color','r')
set(get(ah,'Title'),'String','Handle-Graphics','FontSize',25,'Color','b')



! k1 h7 G# c) E
2 |5 V4 n+ H6 h

CCxiaom

MATLAB 图中图