matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');6 a9 u1 x+ E' E% E. n. B) ]2 i) _, V9 H
   
2.   J5 z$ Z0 p6 {% _2 Y# R# c
   
3. hnb=fgetl(fphn);
   3 ?/ \0 o9 T; u3 m/ j5 [- N6 v\" r
4. 5 \# ?* c. i/ k9 V  d
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
   3 c  ?8 Y2 W  A1 s, H1 H0 f' s9 [
6. 5 N3 P% m) K. O5 J) u/ B
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';( a8 t. ~: r4 S- Y: \4 _9 a
   
2. 
   # S% o% \' c\" W2 W3 J' \$ i/ a, C
3. xa=hnmap(:,[1]);4 M8 r* y1 W, G) S
   
4. , l5 t4 o0 w5 |$ B+ C4 A
   
5. ya=hnmap(:,[2]);
   \" u5 g% T) `& R- l
6. 
   7 t+ D: i) K1 m1 ?7 a8 W

复制代码

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

1. fp=fopen('LATLON57.txt','r');
     ]( d4 `; C4 t
2. 
   # l; \9 W\" s* l2 p+ I7 B' E* R\" `
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
     z9 C6 a% T0 b# s& j: D0 O5 g( ^
4. 6 N- q7 D2 w' R' N9 w9 W
   
5. fclose(fp);
     E% |5 f7 l; d4 W9 Z( H
6. ' B$ o# h( w0 K, O
   

复制代码

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

1. LL57=LL57';
   9 m8 v  S/ j9 g7 H\" i
2. & }3 l9 T/ ?/ N0 q8 Q* Y+ q, o
   
3. x=LL57(:,[3])/10;+ H0 _' T8 c\" I: Z0 Y; ^- w
   
4. ! N# o- A9 L0 s& |# l5 d# @+ [3 f
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');( _1 W1 i, d6 a% U4 [- `# h3 K$ N
   
2. 
   ! @  ^& \* {, `- e. q% r
3. ymd57=fscanf(fpy,'%d',[3,1]);
   ) X1 L' w3 T. ]/ N5 h
4. 
   * X; Y$ o. `1 A. g6 D1 T3 E, J
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);
   % c2 @! t: C; D# ~# l6 p: l
6. 
   4 S/ s- h2 D  S
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   1 L9 F% Z) _% P\" ~
2. . t- T* D+ G( i
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   - [& F+ X. N3 ?
2. , P! i$ y\" a' w1 H$ l7 [
   
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));$ P7 l) x% Z2 G$ l$ A' o5 ?' f
   
2. 
   & K+ R5 n4 F: C+ E- W. B7 I
3. zi=griddata(x,y,z,xi,yi,'cubic');
   1 v. g* @. m+ z6 S0 @
4. 
   8 V3 j, {+ ~$ |7 }4 F$ P- n

复制代码

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

1. hold on;+ v& R9 d7 G. _* j\" o! v  T
   
2. 
   ; ?. B  Z- z4 n: h9 C5 E! C
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);# C+ N: @2 V6 |2 Z
   
2. & ~4 y; s$ H9 m) Z6 e
   
3. hold off;

复制代码

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

最后结果图如下：

具体附件代码如下：

1 U- i; a# |9 h; q' C! [



6 [1 i0 q7 A6 j/ W6 w