matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   8 b) P/ ^9 [6 w. Y9 W/ z
2. * b: G2 f2 A) Y: C0 H* u( d1 t
   
3. hnb=fgetl(fphn);& K) Z4 m. s$ ]; A# }/ s  }4 {+ A; v
   
4. 
   % B5 C$ t5 {8 O6 Y0 g1 u/ W% E! v$ ?
5. hnmap=fscanf(fphn,'%f %f',[2,59]);% {# m% m1 {* E! o8 }9 j( v
   
6. 7 o\" H6 q1 _' C# L
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';
   ! H) P- }/ `, k; U
2. % q0 C, P* D. ~8 s9 ^5 z7 m
   
3. xa=hnmap(:,[1]);; t; M# |6 x7 w) |$ w2 |2 C3 O  R- H
   
4. 
   ! t; |8 d7 y: j( R& n2 J/ a0 y
5. ya=hnmap(:,[2]);1 K  Z; M; _) A( H5 `- F
   
6. 
   - D  v0 k5 ^) b& W& e8 Y

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   3 L( t9 L8 M2 b- {
2.   e, X) d. E! m
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);. k: ^/ {3 R\" t7 _+ U2 \! d
   
4.   g  c4 v& j* q# z# i\" b
   
5. fclose(fp);
   0 O& m3 d& m0 U- a0 P; T( V+ [' g- K
6. , o, S: O! N/ I! Q6 B1 E
   

复制代码

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

1. LL57=LL57';1 f* O5 H2 D) T! ~, ]
   
2. 
   ) ?6 Y$ a+ V: ~
3. x=LL57(:,[3])/10;
   7 }* ~\" I: |$ ~: E0 ]4 |# ]) Q
4. 6 w! P3 X  O4 l, R
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');
   & w, i( T; L$ W  C
2. 
   ; D: H7 ~! Z3 R& X. o' a
3. ymd57=fscanf(fpy,'%d',[3,1]);
   + @3 ^9 ~5 b+ u- z
4. 
   6 ]* n% o( `& J\" v' D1 l
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);
   3 s/ d( B+ R! G/ ?8 [
6. 
     f; |3 O2 a% ~+ Q
7. fclose(fpy);

复制代码

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

1. yu97=yu97';. z0 j' c. V  o4 Y
   
2. - }3 U6 y4 D* o/ d0 V; V1 u
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   * J! G# U# ^4 w9 \4 S8 h: `
2. 
   ! m4 C. c, n+ c) D% y\" D$ J! |
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));3 j) e* C; x: {# ]/ k5 z3 H\" F
   
2. : d( r. ~6 @. j+ O0 Q
   
3. zi=griddata(x,y,z,xi,yi,'cubic');% E/ u  ^- d3 P# V/ |% h
   
4. + z% D  g& h) F4 N( q) E
   

复制代码

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

1. hold on;# J: E$ F. {: P3 a
   
2. 9 a7 D5 E( C5 X% y1 P0 q4 W
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   0 q9 L) |( a3 s' n  n\" I
2. 
   / j/ y9 _) {2 r2 p
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
8 w4 Y: S' W! x9 g- j1 S3 n请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。
: J/ P! `- W  h; v+ z! W  Y

  y: G( M( u$ L- K. H5 ^" P# k

最后结果图如下：

具体附件代码如下：

! E8 ~' z+ f. t9 y  i% U( _) _
/ I0 q- {+ X/ H8 S7 F1 V
; G- e$ Q  c3 o! I: `" y

& h8 i  h) j6 I" h