matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

这两行清空MATLAB的工作区，并关闭所有图形窗口。

1. fphn=fopen('hunan.txt','r');
   9 z, E) J8 K3 V$ b
2. 
   # d0 [* `2 q: \. T8 G
3. hnb=fgetl(fphn);. o5 x$ G8 |' w3 z; C
   
4. 
   & i( a1 ]; j: ]/ p6 \0 o
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
   : V! f( c8 B( {: g9 W  C: h
6. 9 F/ |. {% a! `9 o, @
   
7. fclose(fphn);

复制代码

这部分打开名为 'hunan.txt' 的文件，读取文件中的湖南省的经纬度数据。fgetl 用于读取文件的第一行，而 fscanf 用于读取文件的剩余部分。hnmap 存储了湖南省59个位置的经纬度，每一列对应一个位置。文件的每一行应包含两个浮点数。

1. hnmap=hnmap';
   \" ^; M: ]& _' n/ v0 i0 a
2. 
   2 E0 M! T' {% D5 K
3. xa=hnmap(:,[1]);
   8 i8 n( c1 T# G\" t5 W. E
4. 
   1 t* L2 R9 b% F- @
5. ya=hnmap(:,[2]);( F' L& \' R9 d) b. z- w& X  \
   
6. 8 w: Z( r5 x9 I' Q
   

复制代码

这里对 hnmap 进行了转置，然后将经度和纬度分别存储在 xa 和 ya 中。

1. fp=fopen('LATLON57.txt','r');3 H# Q: H- J$ @! l. B6 s
   
2. * w' z( g# `! P0 t/ c5 ]) {7 p
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);5 }! b7 a0 _- W: S2 d) }) s
   
4. , ~1 n; D2 {5 M- I
   
5. fclose(fp);5 _) O& X\" k, J, S
   
6. * V2 H- w. m\" E. ~% ~/ m& h7 `\" z: B
   

复制代码

这部分打开名为 'LATLON57.txt' 的文件，读取文件中的湖南省97个县的编号和经纬度数据。LL57 存储了97个县的相关信息，每一列对应一个县。文件的每一行应包含一个整数和两个浮点数。

1. LL57=LL57';  E5 P9 q\" Q, C3 |\" V
   
2. 
   % s- m1 j3 }, P0 N3 y4 ]0 t0 \( X
3. x=LL57(:,[3])/10;% D/ w9 B% I& ~4 W
   
4. 
   7 h. r/ S8 l4 Y
5. y=LL57(:,[2])/10;

复制代码

这里对 LL57 进行了转置，然后将经度和纬度分别存储在 x 和 y 中。在这里，经度和纬度都除以10，可能是由于数据单位不同而进行的缩放。

1. fpy=fopen('etw00100.txt','r');$ u4 F' u3 @8 R' a+ n% I
   
2. 
   ) i% i8 C0 r$ C1 W$ I! x
3. ymd57=fscanf(fpy,'%d',[3,1]);
   4 N. s- U. m0 a9 z! k\" P  ~7 ]8 ]) x
4. 
   + S# H+ v4 p5 L8 {
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);; L& F2 v' ~+ w8 r
   
6. . M0 `, g. S7 B$ p
   
7. fclose(fpy);

复制代码

这部分打开名为 'etw00100.txt' 的文件，读取文件中的湖南省97个县的温度数据。yu97 存储了97个县的温度信息，每一列对应一个县。文件的每一行应包含一个整数和两个浮点数。

1. yu97=yu97';
   , q& g% @4 T, _/ o7 r
2. 
   / R1 P$ _/ N3 y0 X; c8 j2 v$ a
3. z=yu97(:,[2]);

复制代码

这里对 yu97 进行了转置，然后将温度数据存储在 z 中。

1. hold on;
   3 [, X( s1 R3 y) h; }- j
2. 2 }) F& F\" w# f( T: d1 H8 e
   
3. plot(xa,ya,'.','markersize',5,'color','red');

复制代码

这部分开始绘制图形，用红色的小点表示湖南省的边界。

1. plot(x,y,'.','markersize',6);

复制代码

这部分继续绘制图形，用蓝色的小点表示湖南省97个县的位置。

1. [xi,yi]=meshgrid(linspace(min(x),max(x),25),linspace(min(y),max(y),25));
   / I5 d, @( _% I0 r, I4 F: m' H
2. 
   8 S: x6 j* s3 }0 ?/ n, ^
3. zi=griddata(x,y,z,xi,yi,'cubic');' i! V$ W* R& }9 C\" |: b
   
4. ! p# q6 r0 x$ B7 p& i- U+ f
   

复制代码

这里创建了一个网格，然后使用 griddata 进行三维数据插值，得到 zi 表示的湖南省97个县的温度分布。linspace 用于生成指定范围内的等间隔数据。

1. hold on;
   % g0 m# J0 j2 Z# c3 l\" k; l
2. & e1 e# f  K* l9 o8 L' z4 C
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

这部分绘制等值线图，表示湖南省97个县的温度等值线。contour 用于画等值线，xi 和 yi 是插值后的网格，zi 是温度数据，'b-' 表示蓝色实线。

1. clabel(c,h);$ ~' ?% g) b4 p' ~\" |( ?
   
2. 
   2 Z. [- G; x. s0 ?7 x3 ~
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。

$ a- R4 Q) o- f

最后结果图如下：

5 a" h8 ~; d& P& j; h0 X1 r

具体附件代码如下：

% {# i# }9 W1 a1 D2 h' h9 n  }
: L; `: S  o9 W
9 U; m* Y9 |/ o+ W" K: S6 S, Y6 {% C