matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   ( v1 v. J7 N1 X6 Y7 }1 o
2. 
   4 ?8 W$ x; w4 G1 a+ o1 l% u
3. hnb=fgetl(fphn);\" S+ K! |7 [' B. T) {0 A4 o) J
   
4. 
   , |) u$ ]7 b7 n0 V
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
   * Y3 Y4 I1 f5 V; o: D
6. ( X5 L7 _, `\" h0 _0 @, H( [
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';1 _& O0 v! m$ h( f3 Q3 }
   
2. # k* }+ x9 ]# g7 b, P\" T+ h9 t: T
   
3. xa=hnmap(:,[1]);
   7 z2 R4 u# i$ a3 j6 W7 C, \
4. 9 Z1 b$ m' w2 O3 P4 m4 ]. ]
   
5. ya=hnmap(:,[2]);
   ' Q7 Z! F4 R# d0 t; a7 b\" Z8 F! q
6. 
   5 ^, B# _8 D* |1 d\" Z

复制代码

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

1. fp=fopen('LATLON57.txt','r');1 x8 |$ @2 F% z/ y2 q
   
2. 
   3 J# P2 Y( Y6 C! ?* y5 s
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   7 D' o5 D/ p, N; K) i* V
4. 
   & }/ l0 |- ^) p
5. fclose(fp);1 h* Y0 D! W/ U# E' @) x; w
   
6. 
   $ o0 O$ I! }: Q2 r# F

复制代码

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

1. LL57=LL57';
   ; P) B5 X4 v0 K9 t! o; R$ O
2. 
   3 Z& W8 n3 b& A$ b
3. x=LL57(:,[3])/10;; M' {1 |! x$ p) C
   
4. / B7 Z+ H\" \# ]' p. Y
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');5 m3 A: V* E( t6 g9 H
   
2. 
   # U$ Z) @' L; B\" J+ s- k
3. ymd57=fscanf(fpy,'%d',[3,1]);\" M2 N1 q, \7 K$ K' I
   
4. 
   3 a% X( `5 c) l& Q
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);5 _9 O8 T& `1 _% ^+ Q
   
6. 7 u: J/ `' a8 V! c* J
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';- ^6 z- c& M$ X  V. J
   
2. ) C3 C1 {+ V9 O% p  U) e. r7 c' n
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;# V1 s  Y2 k/ i\" d7 H' R
   
2. % J4 f\" T3 l0 l
   
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));
   % v4 {; c/ t8 I  t1 ]& T$ @
2. ! s) |0 g4 ]+ a$ z( u
   
3. zi=griddata(x,y,z,xi,yi,'cubic');
   ! d) y8 [# o& F
4. \" n; G' |  T6 A
   

复制代码

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

1. hold on;
   ; o2 B+ A5 T% M; E. V3 t
2. + G  H: `6 }3 u4 m& A& B
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);, p. n: J, {/ F( N/ R. I
   
2. 
   3 z+ L: G# u1 }2 f! v1 m/ O
3. hold off;

复制代码

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

+ u$ e+ h# f/ e3 U) x2 g
7 s& v( q& K6 ^- z  u

最后结果图如下：

  k" P$ D# n. f7 r# W

具体附件代码如下：

, i% G4 P  E' t5 [0 [* v9 {& |; N


6 F1 H% l8 \: ~1 b) |3 s) `

3 u% @% b  x; s/ y