matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');$ |' a6 N# G: W2 z
   
2. ! K9 W0 P/ E5 Z  y& Q' m, p$ T
   
3. hnb=fgetl(fphn);1 R6 O  h0 O, Q1 V8 [! \- `
   
4. \" l2 O; j) J; ^, y& o+ X, q
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);! n2 k! x. L; l4 f
   
6. 2 Q3 E6 i$ v5 }: S
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';- ?# C0 K  n2 ]& ?4 G. [
   
2. ) [  H( t\" x7 x7 V( d1 ]; o, x8 m1 j/ A
   
3. xa=hnmap(:,[1]);: X9 B. B* q3 Z; @' u
   
4. 
   ' F4 G% N\" p+ e3 q\" N
5. ya=hnmap(:,[2]);\" k, `' y+ E% z2 w: l' O/ E\" e
   
6. 
   3 p3 D# h5 |0 A

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   3 V* k. B; c% U; y
2. 
   , y: y- @0 L$ e2 z  S- w\" D
3. LL57=fscanf(fp,'%d %f %f',[3,97]);' M! z+ i1 o9 s+ `9 r2 H2 ?% _
   
4. 6 d( d' P4 }+ n, O
   
5. fclose(fp);
   / e4 c3 Q6 J9 i$ J3 V! R
6. % ~; R1 U1 o6 f/ E0 u! N# I
   

复制代码

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

1. LL57=LL57';
   % F4 d$ l' v3 p; L\" {0 q0 q( D8 D
2. ) U  o$ T6 {  X0 k$ J
   
3. x=LL57(:,[3])/10;
   # J9 k' X\" F2 v6 M! G
4. 
     [, A# u+ p$ J* U- L6 C6 P! Q; ^
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');
   0 ^: I9 F2 O5 Z1 d- G8 O
2. : M: m6 |  u\" a  k1 @
   
3. ymd57=fscanf(fpy,'%d',[3,1]);
   9 ^5 K# u5 `  |) F6 K
4. 
   \" w: I3 z' @\" r) P
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);; h  K! l9 x3 M! x; u: l
   
6. 0 N+ O$ Y' I/ n& Q  N9 X
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';  s/ V# ~9 c! M/ M* y\" o9 z, g
   
2. 1 U, Y' {5 p$ ~1 |5 k  E
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   % N( I3 Y7 C4 f
2. ' D* O4 w\" R, V9 D' Y7 d: R
   
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 m6 d5 F' b( L$ L
   
2. 
   7 k& \  k+ {/ E' E3 Z( }
3. zi=griddata(x,y,z,xi,yi,'cubic');8 u3 h\" A7 b; i0 X$ q5 V  z2 Z
   
4. 
   / v) \/ R& j\" _9 c0 N: K

复制代码

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

1. hold on;
   6 F' ~0 W  T2 P' a! d$ e
2. & G! ~3 A/ z8 A, D& [
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);: |, b( m$ Q# {
   
2. 
   - k4 ~$ R\" e+ {; P) i
3. hold off;

复制代码

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

最后结果图如下：

# U2 y0 V. x0 R+ `! k

* L( H6 W0 ?) E8 e
# y  z0 o8 T0 ~' F7 J# y

具体附件代码如下：

# ?1 d8 [! f, ?5 u: U* H' d

, {/ a1 c$ ]$ F/ y" Q3 |
0 v& b# x" p0 J2 X3 w