matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');- Y2 v* `* V  a0 F# r+ w* ]
   
2. 0 P; e. N1 U0 I6 y( f; y  L
   
3. hnb=fgetl(fphn);* _9 i5 {' ?8 A7 u5 S; h; ]# P
   
4. 9 f( ?* {7 U+ |$ P! G. S6 T
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);/ W- \$ E/ _& ~2 y6 W
   
6. + J  S/ B6 J4 j
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';' |; z, u9 o9 o# n' T& G( Z
   
2. 
   # s! O2 M0 \4 [6 [8 q) \
3. xa=hnmap(:,[1]);
   9 G' e  ]: X& ?
4. 3 G2 o( i, d8 G/ W/ B
   
5. ya=hnmap(:,[2]);
   2 r& v) a# l2 Z, J' \/ i7 r4 w! r
6. 
   ( B- w\" b4 ]0 O, O! h

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   6 K\" r: j5 o! G! X\" ?* u
2. 
   7 b2 b: [  {# G6 a; o+ I
3. LL57=fscanf(fp,'%d %f %f',[3,97]);0 V+ ?+ a  k/ g/ ^$ k' i( e1 e\" F! y
   
4. 
   , O, E) g5 S5 `4 k- K* p& X
5. fclose(fp);- u3 N) |+ O( l; ?) `9 X
   
6. ( y) |+ B; g( U* \% Q( q0 ]
   

复制代码

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

1. LL57=LL57';
   7 {( }! n9 J5 m' [; Y( L9 s
2. 7 c4 C$ Y) `4 ?* @5 p
   
3. x=LL57(:,[3])/10;
   % a6 O3 Y( a6 C$ i
4. 8 S2 D9 P3 K$ }
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');
   $ Y. b1 M\" Y7 l* w) I
2. * G) ?; Y8 _9 `
   
3. ymd57=fscanf(fpy,'%d',[3,1]);+ v4 V3 }' Y  _' ~2 E
   
4. 0 c/ @: d) G2 T5 U\" w( M! y
   
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);3 F' a) ~7 y7 D
   
6. 
   ) P6 m$ s$ c1 n6 e8 [+ d
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   $ j$ a. D( k; u( i
2. 
   + g$ g& q3 o. p4 Y\" I$ |
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   ; h8 `4 j) u, f1 b3 [+ O& f( R
2. 
   ' p: Y% C! i' Y* d9 D* G( E
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 [( R: D' C& o& s( v. i6 H
   
2. % D! @! {7 o; i\" ]# ?. b# U\" l' K
   
3. zi=griddata(x,y,z,xi,yi,'cubic');* g* f  S8 q, J- ^
   
4. 
   9 u6 ~8 p) ?+ z& F0 c/ u

复制代码

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

1. hold on;/ r$ j7 p) B  E! j3 J$ M
   
2. 
   ' x+ D6 p3 A. m7 C' b; x
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);% e  B% k$ ~: N8 g* ^! d
   
2.   ~4 q( I2 ?) e& H0 ]; x; W5 [' W$ \0 I
   
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
' e- v6 w6 g8 M& L1 u请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。
9 J: `; V  C) q  W8 q# R0 P" y! M
# c4 }5 c# \; |4 ?0 T

最后结果图如下：

  P5 l# M7 s3 K5 v
' Y6 Z' e: W# x: ^: o

具体附件代码如下：

! B5 X$ _& q. u* G) U. L



/ P3 b+ O' i; ~/ l7 P" W
, ?* S/ y) K4 }) T5 j# X