matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');( r3 Q  ]2 d* _6 t- q# |4 s
   
2. 
   7 q* b6 F* _$ X0 K( ]+ x6 Q
3. hnb=fgetl(fphn);\" S% L* F) G& Q8 E# o
   
4. 4 O) \) n' E! g, S( e/ c% u& a
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
     W, L# y+ c/ x* R
6. / f) c, w- @* {& j
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';
     {' z# v; }  p+ l
2. 2 C4 C) X/ w5 i  ~: Q3 R% m6 e
   
3. xa=hnmap(:,[1]);
     x  U3 T( l$ C; T7 ?/ [& w
4. 
   & Z' x, S% r! f; a$ y; v( y
5. ya=hnmap(:,[2]);4 x$ F& y! |; r  O0 c  p. T
   
6. \" M1 }: K6 @8 F, w
   

复制代码

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

1. fp=fopen('LATLON57.txt','r');# k# g. g\" c0 e/ X( M' ?
   
2. % j6 _( R6 d6 ^$ U# M5 x; Z6 A* F
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);0 ~& a% |4 W& @' i: Y
   
4. ) i, \' D. }# N, R$ k
   
5. fclose(fp);6 D, i\" E& o# P9 H9 Q5 d- }/ _
   
6. 6 s4 W7 U+ j( l
   

复制代码

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

1. LL57=LL57';8 d) M( {8 K) T6 ^% A8 K9 z
   
2. 
   ' B' j3 \& L% i( z+ y1 p8 O
3. x=LL57(:,[3])/10;. X! ]\" g1 g6 d0 T4 s& `
   
4. 
   1 T$ ~- M\" C& e; C
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');
   ) p* o& P4 n# N4 i' {7 a
2. ! `8 H7 I. L7 x7 T( c1 s+ O' ^4 }
   
3. ymd57=fscanf(fpy,'%d',[3,1]);9 s4 E- @+ I* I* t. P, I2 T% K
   
4. 
   3 q: ]# v# s. u; @+ j9 e' I
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);
   , M; I6 d8 e1 Q/ l
6. , Z4 ~( b# b, F3 T2 {9 l
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   8 \0 p\" R5 N$ p4 V% m5 L: D
2. 
   ) K* l, {, W- P% S3 G
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   5 G% R8 v2 y' M8 f( G
2. 
   ( ?* ^; }+ Z. b9 U. n8 ^
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));( A- K8 {6 L# X\" ~. Z7 f/ ?# T
   
2. * ?( N\" B) }1 \8 i4 W- W! Q
   
3. zi=griddata(x,y,z,xi,yi,'cubic');
   9 @+ ~' |- n8 T; |
4. 
   \" y$ t6 n* G7 h. b2 c% H

复制代码

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

1. hold on;
   5 P! `1 g. `4 K' [& x/ O6 [. ~% Y
2. 
   . ]+ B/ [4 b2 k- Q  \& Z
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   2 K) Y1 `) |1 F7 f2 x' z6 M7 X
2. $ W4 G4 E* s% k, h
   
3. hold off;

复制代码

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

9 C3 i+ H& F) a6 N- B

最后结果图如下：

8 w1 ^0 Y: M  ^: d' ]
/ E# P( @4 Y6 |6 \% T

具体附件代码如下：

: M8 Q8 T; O( Y2 n1 P* A& o, M
# ^! ?& [2 [7 _; k$ g8 U
: w# c1 E5 k8 v* f! U- j! i/ h