|
回复 6楼 zebtra_stripe 的帖子
大约就是这个样, 一些参数慢慢调至LZ要的, 苦功!- load Z; x=Z(:,1); y=Z(:,2); z=Z(:,3);
- axes('Position',[0.1 0.1 0.7 0.8]);
- bb=0:0.001:1; [X,Y]=meshgrid(bb,bb);
- L=griddata(x,y,z,X,Y,'cubic');
- pcolor(X,Y,L); shading interp; caxis([-1 1.5]); % colorbar; colormapeditor
- title('lyapunov exponent-InventoryRate-SLRate?')
- xlabel('InventoryRate'),ylabel('SLRate'),zlabel('lyapunov_exponent')
- cmap=colormap;
- hh=axes('Position',[0.85 0.1 0.1 0.8],'box','off', 'Color','none');
- fill([0.1 0.4 0.4 0.1],[0.1 0.1 0.2 0.2],cmap(end,:)); text(0.5,0.15,'High');
- axis([0 1 0 1]); hold on
- fill([0.1 0.4 0.4 0.1],[0.3 0.3 0.4 0.4],cmap(1,:)); text(0.5,0.35,'Low');
- set(hh,'Xticklabel',[],'yticklabel',[],'Xtick',[],'ytick',[]);
复制代码
[ 本帖最后由 ChaChing 于 2009-12-1 11:58 编辑 ] |
|