matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   ) B3 O( q' C( \+ C  v2 a$ w
2. 
   \" u) F9 y# `9 r' w/ F
3. hnb=fgetl(fphn);
   3 U! `: z6 @* }: ?
4. 6 V6 g% _' X4 ~2 q' \
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
   6 |- _3 P7 r7 u. _; R
6. 
   \" S: O+ o8 b1 ?
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';
   $ k% S\" g- w) ^8 C) s6 f% `3 m
2. 
   $ p( s; e7 W8 r) R# o
3. xa=hnmap(:,[1]);# S' j4 A+ }5 l$ y
   
4. / b# `1 M& B; Z4 Q
   
5. ya=hnmap(:,[2]);: f; a8 ?' v9 I; P
   
6. 
   # Z0 U- J& c+ i- v. M; s: n

复制代码

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

1. fp=fopen('LATLON57.txt','r');* |! v: C( R4 o\" Q( l$ e/ [
   
2. / G\" T  w0 g+ ]4 m& R$ s
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);! l& m* F2 q\" ?
   
4. 
   $ t0 X6 L- l7 q8 H( m, b
5. fclose(fp);; b2 e6 J/ t7 W3 y# R
   
6. + W% A) |8 {2 v7 V8 `0 D
   

复制代码

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

1. LL57=LL57';, J  [- P9 m+ j. t0 V
   
2. 
   2 _2 Z: D0 L: j1 p
3. x=LL57(:,[3])/10;
   ! @* |) S  z3 `  l+ W9 ?* h5 w
4. 2 E, _& f& c3 c5 Z3 O7 U! i
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');
   / F% e7 a8 y$ u0 u6 N
2. 8 L* {! j5 D$ w! y& a: X* r
   
3. ymd57=fscanf(fpy,'%d',[3,1]);
   \" R% O% B; C$ {
4. ( I% @' E$ B6 H% u; w
   
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);\" d6 H, j4 E! Q. m* l2 D
   
6. ; N% T+ c1 j2 a9 j, Q  l* Z; G: X
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   3 {+ L. @; W% C0 k' k
2. 
   $ c' p8 s! Q  }5 V8 }9 ^0 w4 W$ ^
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   \" S4 K+ v4 [) N, e; v1 G& O# L
2. 
   ; v1 k+ y1 }0 Z. r7 Y6 K
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));
   ) T4 ~$ X9 O\" g- f. D- {
2. 
   % @7 A# }& |8 i+ i; q, i2 i
3. zi=griddata(x,y,z,xi,yi,'cubic');
   & L9 i+ j5 C7 z2 |& `. p- B6 }4 D
4. - {* |' m\" G4 D7 \5 I
   

复制代码

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

1. hold on;
   1 L0 p7 ^8 t3 c
2. 
   ' d9 a/ l7 q+ ]; A  p
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   % @# C, o4 D6 b1 z
2. 
   \" Q1 i( [; |1 _: Q
3. hold off;

复制代码

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

最后结果图如下：

; ?0 m. c" ^0 J& L) \" k: {

具体附件代码如下：

+ c7 W0 E0 k, R( O" ^
' n$ q0 a8 U. A
- P5 D! R+ L* ^4 V' \7 H, c6 q& f
) e7 v: W% Z5 i- X
4 |6 ?6 [* P" B& Y& D+ n3 G