matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   5 [  M. W4 y0 Q( h9 }
2. - `\" S/ G\" e  d7 d8 Y; {/ I
   
3. hnb=fgetl(fphn);) X, ^\" W7 v) e' e9 s9 R
   
4. . P! d8 b# S/ T5 s/ }- R
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
   + T1 V  ]! |1 s% w% F- M
6. 1 c2 f$ x) o1 P4 D1 i
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';: N9 u# f* y- f, ~
   
2. 
   8 p5 J' ~2 K7 h0 U( y; g3 n
3. xa=hnmap(:,[1]);4 d$ E6 R& p( @! S
   
4. 
   3 l( z6 ?( D\" o9 v; C3 {4 k7 r
5. ya=hnmap(:,[2]);  P8 X$ x, e4 P) q4 N1 @0 G
   
6. 
   4 X+ C- [; ?\" {; U

复制代码

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

1. fp=fopen('LATLON57.txt','r');( d. [) g8 T! g! [; H) C\" }
   
2. 
     l; S; f- X\" b& Z
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   9 E- G\" s- E5 w+ K$ c
4. 
   ' A6 L/ c$ S9 {- \/ |
5. fclose(fp);
   + ~+ {: Z* z9 C4 ~) m; @, X
6. 0 q1 M/ R& H/ N4 k
   

复制代码

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

1. LL57=LL57';
   + y/ V3 t( V4 u& A4 Q$ g
2. 
   : ^- [% L- l3 H% e+ t0 d
3. x=LL57(:,[3])/10;
   & u6 f\" O  @4 }* G3 N' }3 v
4. 
   . y\" w5 Q4 o  l2 W2 l( g/ ~
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');, _5 L% r, F' {: n
   
2. 4 a% n( p5 U8 e( l
   
3. ymd57=fscanf(fpy,'%d',[3,1]);
   / ~7 w' n5 |$ Z) H- L
4. 
   / n0 v3 Y4 j7 g\" Q1 j# ~  B% q
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);5 D) R\" k7 A' x( B  O1 W  g
   
6. . b: M! R# d* N7 `
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';( b& n& C8 a9 b0 d3 }$ m9 c$ n7 b
   
2. 
   ) w- I8 k$ U7 o1 W7 F9 ^
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   % ?8 u8 F- k5 c1 F
2. - `! D% F5 ^5 q1 c\" P+ F+ \
   
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));
   & r4 P0 g% [; {5 |. }9 I9 z, Y8 y4 \
2. 3 k% N0 r/ t' {+ Q8 [0 G7 T
   
3. zi=griddata(x,y,z,xi,yi,'cubic');
   9 f7 {* l/ |9 ^% Z* K+ \
4. # \3 |3 n5 d& {1 S- _
   

复制代码

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

1. hold on;, l/ V* T  a# U8 C
   
2. 
   8 y\" _3 S% @7 t: L; M) K) l
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);/ C0 q! C9 z; W
   
2. 
   6 ~- U: a' |4 v' e& Y' k4 f
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。
1 F7 ~4 R$ m. `) o! l, R  t
/ n0 x) f9 [1 [1 K
/ U4 ~4 b# ~* s" t

最后结果图如下：

, |$ f$ Y: T+ K6 A+ A9 d

具体附件代码如下：

' g) r2 e! f; {