matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');- F# e+ N# v  T/ [+ e6 S
   
2. 
   \" P4 I\" n\" W% z8 Z8 N5 f
3. hnb=fgetl(fphn);2 D/ K: y: q, C) }( h
   
4. ( G6 D# |\" q/ e- }
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);\" k' N6 Z1 Z5 i
   
6. 
   ' q. A2 N% J; ?
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';
   # f+ z# Y5 N% |$ Z* Y' Z
2. 9 ]' Z1 I) k. |\" ~\" d# B
   
3. xa=hnmap(:,[1]);+ n9 m6 i2 B& d- Q4 }6 J
   
4. * p$ ~2 S8 ^) w8 f& g; K' ~
   
5. ya=hnmap(:,[2]);: _/ L: K* d1 a# v0 d
   
6. 8 m) u, K; I% Z, w9 X
   

复制代码

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

1. fp=fopen('LATLON57.txt','r');; `8 t$ `3 W3 @\" N\" r\" c) o& T' q
   
2. 
   0 \  W6 z, C+ ]\" \
3. LL57=fscanf(fp,'%d %f %f',[3,97]);* L# m0 [4 c' d
   
4. 
   ) H. e9 T7 @( k# U
5. fclose(fp);; a$ k/ s: y( c: P* C4 [
   
6. 0 W3 j, M  X7 L
   

复制代码

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

1. LL57=LL57';
   . O& v1 l& G: M  }- F
2. 
     y1 B2 Z, l5 i8 ^
3. x=LL57(:,[3])/10;
   . r4 A- _4 y; ~\" Z
4. 
   6 ]3 Y/ X6 f- V' E
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');  V3 L1 V! R\" u$ k
   
2. 
   $ K# a' i& l4 }
3. ymd57=fscanf(fpy,'%d',[3,1]);& r+ U0 E! g\" B  k2 c
   
4. 
   \" p( a& E. u2 o$ W6 g
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);& N! q2 ^) w( B  l# V
   
6. 6 e! o& @\" P' O
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   , m0 b) m3 @$ `+ R( G6 S# x( z
2. ) b& m& z9 m$ y  x1 a
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   % h5 k: u& u  P' `
2. 
   # a' f2 W8 |* o9 I* g: 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));
   ! Q7 \; n5 l, i6 z
2. % X! q; T+ O, V+ W7 u
   
3. zi=griddata(x,y,z,xi,yi,'cubic');
   ' a9 a. h0 ~7 N/ W* H) T& L
4. 
   # w3 X7 e: f9 L8 i* s6 c  ^! m

复制代码

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

1. hold on;2 G' K# ~& [, K. d: i9 P
   
2. 
   . C3 D7 e3 o) b8 z+ k5 L* A- R
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);; K1 n! M# u- Q, f
   
2. 0 @4 I  D* G- y
   
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。
' r8 l; n6 x2 u" L8 s- u2 W& X: f! p6 U

最后结果图如下：

) p+ e) m- U0 }
. ~" P0 `- t& j

具体附件代码如下：

. }7 p/ @; C; g
9 w0 a0 ?& J/ ]$ ]+ a, D


; L/ X. {: ^3 [5 _9 P+ `
- w) l: l$ f, u" M