matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');, l; f( I0 [8 F
   
2. 7 `- U! w/ x; B2 }\" l\" L
   
3. hnb=fgetl(fphn);
   \" s+ {, v. M- G; B* j
4. 
   ! ^: `+ l- M, r) [* C8 t2 k2 Y
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
   5 a$ A8 K& Q5 B  P! g
6. + f6 q& X; Q0 s: D
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';8 y; z$ G' p0 x7 [9 z2 c3 T
   
2. - h: t3 J/ T- i
   
3. xa=hnmap(:,[1]);0 P# _8 k) q- V0 x
   
4. 
   : I2 Y3 O$ ?0 f( z0 F+ L: L
5. ya=hnmap(:,[2]);% R- B* m0 u1 ?1 o1 ~! ]* x
   
6. # C. k' k- E/ b' y4 ~& s
   

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   3 f\" _5 B, x# C5 |
2. - x8 h6 @2 f1 F! Q: y
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   ' C4 \2 H: v% F) A3 r
4. 
   2 ^/ D4 \8 ]1 j; d8 @/ U/ ]0 V- O
5. fclose(fp);+ m- I6 E- R: Q# N1 N4 \0 k
   
6. % @# z+ I+ m/ G0 F\" S5 [+ Z5 ?
   

复制代码

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

1. LL57=LL57';
   * c$ E$ z: M( G$ \8 q3 E* v3 s: N
2. ) `# W8 W( B* P% e8 M
   
3. x=LL57(:,[3])/10;
   9 X% P3 `! g- q$ P' I$ ^) K
4. 
   ! F+ c& v, C* l5 s  D, y5 {  e& D) {
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');1 M. Y- P( |8 d! ]/ X) o4 w$ s0 z
   
2. * z\" P# Q  B; h* S  P5 g
   
3. ymd57=fscanf(fpy,'%d',[3,1]);
   - A6 a  _7 U- P5 x& n; A6 t) X; b
4. 
     Q  f; H! f! Y. T
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);
   % u& `( U- h/ C5 b* d' X, G! a
6. 8 r\" e- }5 O# k
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';5 o* m2 }3 S* D+ O. B# g
   
2. 
   & x6 F7 Q2 S6 Y
3. z=yu97(:,[2]);

复制代码

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

1. hold on;6 Q, ^/ Z  Z\" i' c8 r1 X9 t
   
2. 5 ^( M' E2 A% c
   
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));+ ]2 g' @9 t  |: o8 D' Y
   
2. 
   % G3 u* i; A( }: T
3. zi=griddata(x,y,z,xi,yi,'cubic');
   % {( w* J$ i, f; x# w1 K% C, c* E
4. # n, x, v0 X8 P! ]' m. s$ {5 l
   

复制代码

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

1. hold on;- u8 w; A6 ~3 q4 ]! t
   
2. 
   6 ~' D0 n' z+ Q1 e1 |# B
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   ; K0 f0 e7 U, g
2. 
   ' h! b& g) l# ~, K& T
3. hold off;

复制代码

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

5 u5 \& k8 e0 j- \7 S  p

最后结果图如下：

具体附件代码如下：

& l/ d' X1 x! `1 l3 x5 |1 X
! \0 z5 @. }  b' q/ u$ _8 |


5 Q/ p) p% z$ a5 W1 K