分享

注释版本的 matlab 画图程序

 永恒_世界 2014-11-13

注释版本的 matlab 画图程序

(2010-08-11 20:02:25)
标签:

杂谈

分类: matlab

%%%%%常用选项和小技巧%%%%%

%画等值线
[cc hh]=contour(peaks(30),'LINESPEC','b-')
clabel(cc,hh,'manual')

%写文本
text(5,10,'bf math sl math it math rm math alpha','color',[0.1 0.1 0.9],'fontsize',24)

%设置线宽
set(gca,'linewidth',2)

%写标题并设置字体的大小
ti=title('Title of My Figure','color','blue')
set(ti,'fontsize',24)
drawnow

%输出文件
print -dpsc plotE.ps
print -append -dpsc plotE.ps
print -djpeg100 plotE.jpg

%给定图窗口标题
figure('Name','My Figure1')

%设置使用调色板
map=hsv(32);
colormap(map)

%设定等值线的范围
caxis([-6 6])

%设定colorbar的方向和位置
hc=colorbar('hori');
po=get(hc,'position');
%set(hc,'position',[po(1) po(2)+0.14 po(3) po(4)-0.01]);
%set(hc,'XLim',[0 300]);

%设置绘图的缺省值
set(0,'DefaultLineLineWidth',2)
set(0,'DefaultAxesFontSize',18)
set(0,'DefaultAxesLineWidth',2)
set(0,'DefaultAxesTickLength',[0.01 0.025])
set(0,'DefaultPatchLineWidth',2)
set(0,'DefaultSurfaceLineWidth',2)
set(0,'DefaultRectangleLineWidth',2)

set(0,'DefaultLineLineWidth','remove')
set(0,'DefaultAxesFontSize','remove')
set(0,'DefaultAxesLineWidth','remove')
set(0,'DefaultAxesTickLength','remove')
set(0,'DefaultPatchLineWidth','remove')
set(0,'DefaultSurfaceLineWidth','remove')
set(0,'DefaultRectangleLineWidth','remove')

%%set(0,'DefaultLineLineWidth','factory')
%%get(0,'factory')

%设置坐标轴的间隔和显示
set(gca,'XTick',[0:1:25],'YTick',[0:300:6000])
set(gca,'XTickLabelMode','manual')
set(gca,'XTickLabel','0|||||5|||||10|||||15|||||20|||||25|')

%分别设置等值线和等值线标注的间隔
[cc hh]=contour(peaks(30),[-6:2:8])
c(cc,hh,[-6:4:8],'fontsize',15,'color','r','rotation',0,'labelspacing',200)

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%几个实例%%%%%

%给定月份名字
name='JanFebMarAprMayJunJulAugSepOctNovDec'
%给定每层深度
depth=[5 15 25 35 45 55 65 75 85 95 110 130 155 185 220 260 305 365 450 575 755 1115 1735 2615 3645 4830]

%给定经纬度
x=31:1:288;
y=-61:1:61;
[lo la]=meshgrid(x,y);

%等值线间距
vt=2:2:30;
vs=30:0.5:36;
vc=10:5:200;

%打开文件读取数据
h1=fopen(['plotE' '.dat']);

u=fscanf(h1,'%f',[258 123]);
v=fscanf(h1,'%f',[258 123]);
t=fscanf(h1,'%f',[258 123]);
s=fscanf(h1,'%f',[258 123]);
ccc=fscanf(h1,'%f',[258 123]);

%设置流速标尺
u(20,220)=0.5;

%%%%%%%%figure1%%%%%%%%等值线
figure;
m_proj('Lambert','lon',[31 288],'lat',[-61 61]);
hold on;
[cc hh]=m_contour(lo,la,ccc);
c(cc,hh,'fontsize',6);
title('carbon concentration');
m_coast('patch',[0.1 0.1 0.8])
m_grid('linestyle','none','box','fancy','linewidth',2,'XaxisLocation','bottom')
print -dpsc mmap.ps

%%%%%%%%figure2%%%%%%%%矢量图
figure;
m_proj('Equidistant','lon',[31 288],'lat',[-61 61]);
hold on;
m_coast('patch',[0.3 0.5 0.3]);
m_grid('linestyle','none','tickdir','out','linewidth',2,'XaxisLocation','top')
hq=m_quiver(lo,la,u',v',2)
set(hq,'color','b')
m_text(240,54,'0.5 cm s^{-1}','color','r','fontsize',10)
print -dpsc -append mmap.ps

%%%%%%%%figure3%%%%%%%%使用高分辨率岸线
figure;
m_proj('Equidistant','lon',[31 288],'lat',[-61 61]);
hold on;
m_contour(lo,la,t)
hq=m_quiver(lo,la,u,v,2)
set(hq,'color','b')
m_text(240,54,'50 cm s^{-1}','color','r','fontsize',10)
m_coast('patch',[0.1 0.1 0.8]);
m_grid('linestyle','none','tickdir','out','linewidth',2,'XaxisLocation','top')
%m_gshhs_i('patch',[0.1 0.1 0.8]);
%m_gshhs_h('patch',[0.1 0.1 0.6]);
%m_gshhs_h('save','gumby');
%m_usercoast('gumby','patch',[0.1 0.1 0.6]);
print -dpsc -append mmap.ps

%%%%%%%%figure4%%%%%%%%直线图
figure;
m_proj('Equidistant','lon',[31 288],'lat',[-61 61]);
hold on;m_coast('patch',[0.1 0.1 0.8]);
lo1=[180:5:270]
la1=[-30:5:60]
m_line(lo1,la1,'marker','.','markersize',20,'linestyle','none')
m_grid('linestyle','none','tickdir','out','linewidth',2,'XaxisLocation','top','xtick',12)
print -dpsc -append mmap.ps



    本站是提供个人知识管理的网络存储空间,所有内容均由用户发布,不代表本站观点。请注意甄别内容中的联系方式、诱导购买等信息,谨防诈骗。如发现有害或侵权内容,请点击一键举报。
    转藏 分享 献花(0

    0条评论

    发表

    请遵守用户 评论公约