1、MATLAB如何从硬盘读取文件。

2、如何微调subplot子图的位置。

3、plot命令绘曲线时，曲线上的标志如何调整大小。

4、坐标轴的调整。

6、坐标标题中如何标上标。

7、如何调整图示(legend)的位置。

%----------------------------------

%

这里要画一个2*2共4幅子图。先将第1个子图的位置调整。

h = subplot( 2, 2,

1); % 先让MATLAB默

认绘制第1幅子图，h是子图1的句柄

po = get( h, 'Position'

); % get命令从句柄h中获取'Position'的内容，返回一个含4个元素的一维数组放到po中。这4个元

素分别是子图1的left, bottom, width, height。

subplot( 'Position', [po(1)+0.03, po(2)-0.03, po(3),

po(4)]); 子图1的新位置可以这样调整

%----------------------------------

hold on;

axis([0 13 -3 2]);

set( gca, 'XTick', [1:12]

); gca表示当前对象句柄，set命令分别对当前对象(即子图1)设置坐标轴XTick和YTick属性。这

两个属性分别表示了坐标轴的实际绘值范围。

set( gca, 'YTick', [-3:1:2] );

title( 'The North Hemisphere' );

plot( 1:12, bc, '-r.', 'MarkerSize', 10

); 子图1中第1条曲线用实线绘，带有圆点，红色。MarkerSize属性设

置圆点的大小是10。这样画出来的就是实心圆了。

plot( 1:12, nit, '-b.', 'MarkerSize', 10 );

plot( 1:12, sul, '-g.', 'MarkerSize', 10 );

plot( 1:12, poa, '-m.', 'MarkerSize', 10 );

plot( 1:12, soa, '-k.', 'MarkerSize', 10 );

%zeroArr = zeros( 14 );

%plot( 0:13, zeroArr, '--k' )

xlabel( 'Month' );

ylabel( 'Radiative Effect (Wm^-^2)'

); 单位里有上标，^表示后续一个字符为上标。

下述代码绘子图2、3、4，雷同。

%--------------------------------------------------------------------------

% NH Radiative Forcing

Fut-Mod 子图2

fid_bc=fopen('D:_CurrentPaperRadiativeForcingFutMod_NHbc.dat','r');

bc = fscanf( fid_bc, '%f', [1,12]);

fclose( fid_bc );

fid_nit=fopen('D:_CurrentPaperRadiativeForcingFutMod_NHnit.dat','r');

nit = fscanf( fid_nit, '%f', [1,12]);

fclose( fid_nit );

fid_sul=fopen('D:_CurrentPaperRadiativeForcingFutMod_NHsul.dat','r');

sul = fscanf( fid_sul, '%f', [1,12]);

fclose( fid_sul );

fid_poa=fopen('D:_CurrentPaperRadiativeForcingFutMod_NHpoa.dat','r');

poa = fscanf( fid_poa, '%f', [1,12]);

fclose( fid_poa );

fid_soa=fopen('D:_CurrentPaperRadiativeForcingFutMod_NHsoa.dat','r');

soa = fscanf( fid_soa, '%f', [1,12]);

fclose( fid_soa );

%----------------------------------

h = subplot( 2, 2, 3, 'replace' );

po = get( h, 'Position' );

subplot( 2, 2, 3, 'replace' );

subplot( 'Position', [po(1)+0.03, po(2)+0.03, po(3), po(4)]);

%----------------------------------

box on;

hold on;

axis([0 13 -3 2]);

set( gca, 'XTick', [1:12] );

set( gca, 'YTick', [-3:1:2] );

%title( 'NH Fut-Mod' );

plot( 1:12, bc, '-r.', 'MarkerSize', 10 );

plot( 1:12, nit, '-b.', 'MarkerSize', 10 );

plot( 1:12, sul, '-g.', 'MarkerSize', 10 );

plot( 1:12, poa, '-m.', 'MarkerSize', 10 );

plot( 1:12, soa, '-k.', 'MarkerSize', 10 );

%zeroArr = zeros( 14 );

%plot( 0:13, zeroArr, '--k' )

xlabel( 'Month' );

ylabel( 'Radiative Forcing (Wm^-^2)' );

%--------------------------------------------------------------------------

% SH Radiative Effect

Mod-Noall 子图3

fid_bc=fopen('D:_CurrentPaperRadiativeForcingModNoall_SHbc.dat','r');

bc = fscanf( fid_bc, '%f', [1,12]);

fclose( fid_bc );

fid_nit=fopen('D:_CurrentPaperRadiativeForcingModNoall_SHnit.dat','r');

nit = fscanf( fid_nit, '%f', [1,12]);

fclose( fid_nit );

fid_sul=fopen('D:_CurrentPaperRadiativeForcingModNoall_SHsul.dat','r');

sul = fscanf( fid_sul, '%f', [1,12]);

fclose( fid_sul );

fid_poa=fopen('D:_CurrentPaperRadiativeForcingModNoall_SHpoa.dat','r');

poa = fscanf( fid_poa, '%f', [1,12]);

fclose( fid_poa );

fid_soa=fopen('D:_CurrentPaperRadiativeForcingModNoall_SHsoa.dat','r');

soa = fscanf( fid_soa, '%f', [1,12]);

fclose( fid_soa );

%----------------------------------

h = subplot( 2, 2, 2, 'replace' );

po = get( h, 'Position' );

subplot( 2, 2, 2, 'replace' );

subplot( 'Position', [po(1)-0.03, po(2)-0.03, po(3), po(4)]);

%----------------------------------

box on;

hold on;

axis([0 13 -1.2 0.8]);

set( gca, 'XTick', [1:12] );

set( gca, 'YTick', [-1.2:0.4:0.8] );

title( 'The South Hemisphere' );

plot( 1:12, bc, '-r.', 'MarkerSize', 10 );

plot( 1:12, nit, '-b.', 'MarkerSize', 10 );

plot( 1:12, sul, '-g.', 'MarkerSize', 10 );

plot( 1:12, poa, '-m.', 'MarkerSize', 10 );

plot( 1:12, soa, '-k.', 'MarkerSize', 10 );

%zeroArr = zeros( 14 );

%plot( 0:13, zeroArr, '--k' )

xlabel( 'Month' );

%ylabel( 'Radiative Effect (Wm^-^2)' );

%--------------------------------------------------------------------------

% SH Radiative Forcing

Fut-Mod 子图4

fid_bc=fopen('D:_CurrentPaperRadiativeForcingFutMod_SHbc.dat','r');

bc = fscanf( fid_bc, '%f', [1,12]);

fclose( fid_bc );

fid_nit=fopen('D:_CurrentPaperRadiativeForcingFutMod_SHnit.dat','r');

nit = fscanf( fid_nit, '%f', [1,12]);

fclose( fid_nit );

fid_sul=fopen('D:_CurrentPaperRadiativeForcingFutMod_SHsul.dat','r');

sul = fscanf( fid_sul, '%f', [1,12]);

fclose( fid_sul );

fid_poa=fopen('D:_CurrentPaperRadiativeForcingFutMod_SHpoa.dat','r');

poa = fscanf( fid_poa, '%f', [1,12]);

fclose( fid_poa );

fid_soa=fopen('D:_CurrentPaperRadiativeForcingFutMod_SHsoa.dat','r');

soa = fscanf( fid_soa, '%f', [1,12]);

fclose( fid_soa );

%----------------------------------

h = subplot( 2, 2, 4, 'replace' );

po = get( h, 'Position' );

subplot( 2, 2, 4, 'replace' );

subplot( 'Position', [po(1)-0.03, po(2)+0.03, po(3), po(4)]);

%----------------------------------

box on;

hold on;

axis([0 13 -1.2 0.8]);

set( gca, 'XTick', [1:12] );

set( gca, 'YTick', [-1.2:0.4:0.8] );

%title( 'SH Fut-Mod' );

plot( 1:12, bc, '-r.', 'MarkerSize', 10 );

plot( 1:12, nit, '-b.', 'MarkerSize', 10 );

plot( 1:12, sul, '-g.', 'MarkerSize', 10 );

plot( 1:12, poa, '-m.', 'MarkerSize', 10 );

plot( 1:12, soa, '-k.', 'MarkerSize', 10 );

%zeroArr = zeros( 14 );

%plot( 0:13, zeroArr, '--k' )

xlabel( 'Month' );

%ylabel( 'Radiative Forcing (Wm^-^2)' );

我将legend放在了子图4上。

gca=legend( 'BC', 'Nitrate', 'Sulfate', 'POA',

'SOA', 4

); 4表示把legend放在子图的右下角，还有几个数字的含义是：

0

= Automatic "best" placement (least conflict with data)

1

= Upper right-hand corner (default)

2

= Upper left-hand corner

3

= Lower left-hand corner

4

= Lower right-hand corner

-1

= To the right of the plot

po=get( gca, 'Position'

); 发现这样放置后legend要挡住图，因此需要再微调一下。获得legend的'Position'值。

set( gca, 'FontSize', 8, 'Position', [po(1)-0.01, po(2)+0.01,

po(3), po(4)]

); 重新设置legend的位置，同时设置legend里面的字体为8号。

legend('boxoff'); 不画legend的外框。

强调的是上述调整legend的值要不断地试。因为legend相对子图的位置还要随画图窗口大小变

化而变化。如果你看不懂这句，试试就知道了。

我一般是将MATLAB画出的图打印成PDF，再用Acrobat打开截屏，贴到WORD中，这样图

像质量似乎比较好。谁还有更好的将MATLAB图转贴到WORD的方法，欢迎赐教。