matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');& g6 N/ w( M7 m; ~4 p1 f! I0 m
   
2. ! V+ w. j/ o9 u+ o
   
3. hnb=fgetl(fphn);9 D7 \+ \# x  ~# P% w2 d- l
   
4. 
   ; T( i% ]+ F' X
5. hnmap=fscanf(fphn,'%f %f',[2,59]);' q& @9 q  x. n( R5 c) j; {
   
6. 
   7 n- V\" k$ X. L3 H: a5 d# o
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';
   ) e3 p0 g9 }7 _3 z3 \
2. . a# j. y# H5 q8 d
   
3. xa=hnmap(:,[1]);
   & F$ ~. v: _' x; n2 x
4. - H( @3 P, D2 e/ T  u
   
5. ya=hnmap(:,[2]);
   & `( r: d; V! E- }5 d9 w% ?
6. 1 Z9 C  B0 F# j, p; _3 P
   

复制代码

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

1. fp=fopen('LATLON57.txt','r');1 O$ Y* E  j0 C  P6 N2 l\" x: D/ m
   
2. ' q% _\" g3 F\" n; U  T
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   ; x5 Q! W3 ]' `3 E1 ?8 c
4. - E( n$ w1 L- S* ]\" G% x
   
5. fclose(fp);
   \" ~; u& s, s0 Z+ u' W1 F4 F
6. 8 }3 c1 _\" B8 g' z+ i\" `. K
   

复制代码

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

1. LL57=LL57';3 o' T0 S/ X+ T3 s$ z! n
   
2. 
   ! e; C3 V( V8 \$ k0 a
3. x=LL57(:,[3])/10;$ \9 m6 ^6 \5 G\" C- W
   
4. 
   : A2 V/ i- L  N4 S0 H' o) l
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');
   9 s4 q& W\" L# c
2. 
   & p% \9 V0 E, U9 m3 B6 p- A, p: _/ Q
3. ymd57=fscanf(fpy,'%d',[3,1]);2 ?3 D( V3 ?7 y/ J4 k
   
4.   H! z# }% d3 K8 v9 C) L* P
   
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);  W6 Z& J- f9 I9 Z7 |
   
6. 
   1 w% Z0 I3 ]3 \; E: Q
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   7 N' q3 [9 p3 P% H- P
2. ( J+ b5 z3 R- c! U7 J9 O
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   7 Z5 o6 r/ Q0 ~' B; O
2. 
   ( o3 L5 k8 b9 L' y
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));
   7 g: F) D' _+ d7 V
2. 
   3 l: h7 t3 @% X
3. zi=griddata(x,y,z,xi,yi,'cubic');
   $ f% E% \% ~\" s- i2 I
4. - g' \) X: I& n( U; M7 L0 {
   

复制代码

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

1. hold on;
   ( H0 A8 j) ]% n  x3 D1 ?8 k
2. 
     c( @8 ~5 ?\" n% y1 u7 |
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   , @: y! Z2 F3 \' P1 ?+ s# L( ]
2. + |3 z$ A' p# U& D3 j+ H. Q8 O. a
   
3. hold off;

复制代码

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

9 Z/ N+ W4 p  d, o/ [

最后结果图如下：

* _, J* ?1 h' s5 k

具体附件代码如下：

( U; m! {5 N" r


7 W0 x6 {; f2 B; ]
7 n* E' O& F# [3 V: @; }8 d
6 x2 T8 k( |' X0 A( `