matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   5 ?' A1 m* `' ^3 ^
2. , b. g  R6 c  H\" r, B
   
3. hnb=fgetl(fphn);- X2 _) z1 _5 h: K5 M\" s
   
4. : l  ~# Z2 P: h
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);2 v! b+ @3 I2 X; d1 x
   
6. ' m  k0 I$ D5 h4 ^( S. T
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';0 F: h. g4 Z# V0 r& j+ [
   
2. , w5 z* P6 G5 w$ f
   
3. xa=hnmap(:,[1]);2 t% N) n6 `6 `. ^* E& q
   
4. 
   ! c) n6 T\" ~9 O$ K' M$ S) I
5. ya=hnmap(:,[2]);$ b9 p4 W4 }; D# B
   
6. \" Q6 r) ^3 l0 r7 m4 I
   

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   0 P8 e; c: L8 t+ ?% \4 N4 Z
2. * d, K+ c+ j5 T0 w( Y1 B
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   3 [0 p/ E* S\" j2 \
4. . E+ {# I( o5 H6 b& c) _/ C( |
   
5. fclose(fp);
   # c$ J- P\" b8 e1 P' M
6. ; V5 P3 [, S( N1 [
   

复制代码

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

1. LL57=LL57';
   \" L! |- y/ O7 [  @7 ^- G- G) T
2. 
   ( x- Y; {) K- H\" r& \& ]) G, q, |  S
3. x=LL57(:,[3])/10;
   ; n) j8 l7 {2 T1 M; m, k
4. ; b6 A: t8 g# H9 _% x
   
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');7 L0 q! Q! `9 j, n
   
2. 7 G7 E6 L% G( N
   
3. ymd57=fscanf(fpy,'%d',[3,1]);
     Y\" b# Y0 @2 c
4. 
   9 l/ G$ O. |, V1 @9 h/ }5 n' e+ ?1 C
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);
   3 S6 S* ?: i+ |( [8 B
6. 
   . J6 h% q- H# s; \\" J0 t
7. fclose(fpy);

复制代码

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

1. yu97=yu97';. ^; C7 J7 @1 ]
   
2. ; M0 R1 R  z/ c( n9 d# v3 S7 e
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;5 R+ t/ W5 ~9 v4 `9 e4 S
   
2. \" t8 {0 p0 I+ M; E7 f
   
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));- l; A  _  `: a/ S$ o; K% [
   
2. 
   \" ]! k! V, U: u( s, L9 L; A2 c2 @2 m- s
3. zi=griddata(x,y,z,xi,yi,'cubic');
   ) K8 J# B( z* X$ O2 ^' t/ y
4. : O' F. I% w$ B' o0 Z# Q
   

复制代码

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

1. hold on;1 n# m2 J' \5 B$ ~. H8 i
   
2. 
   7 {+ X% J0 e- C: m
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   3 q4 x1 O4 q: V8 w\" P+ N' X
2. 
   ; L9 P4 n& q. V* b/ c% c
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。
6 I6 j) a; {4 @
6 ?% l9 W) B" Y4 Z4 x3 R, K1 k

最后结果图如下：

  B4 t! d/ h+ B% j2 E1 t

: t" A2 u$ q4 ?# Z0 f0 ~9 a' M

具体附件代码如下：

6 C8 R+ R7 }+ U% w2 ^. j. [
* J2 j+ M" {( f/ a. I7 }4 L1 Z( X


( |5 _  p# `7 t3 N8 b