matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   , g- K5 W2 ^9 Q5 s( ^9 ?
2. 9 c, y7 H# b  \% Q2 j: ^9 [
   
3. hnb=fgetl(fphn);
   6 B, {4 W7 V/ Z3 f9 r: U
4. 
   3 k( G/ G' ^6 j0 {8 {) z, S! s& T
5. hnmap=fscanf(fphn,'%f %f',[2,59]);6 q$ E/ @; I$ N
   
6. 
   \" _- x# u0 J/ t6 d, C4 K
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';5 t; X7 Z: O. ?, m/ g1 ~' D
   
2. 
   4 R# R3 R9 a+ B4 [
3. xa=hnmap(:,[1]);1 i4 \4 K! g( F7 b1 }
   
4. ! A5 J- z( ?4 c
   
5. ya=hnmap(:,[2]);
   . W. W( t( Z* e1 a
6. : C' A9 W/ o$ ^5 s2 L. X
   

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   5 q: k: H/ p7 s0 g7 Q5 A
2. 
   - `6 y2 F; c1 q) F! g
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   0 \8 j\" t& L$ G
4. 1 [2 q/ J4 i( k' r
   
5. fclose(fp);
   8 H( D% l! K$ y% z( k+ l7 X. F7 D\" f  R
6. % A# Y, ?. M) I8 p
   

复制代码

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

1. LL57=LL57';
   ' u  p7 W2 q- Q+ L7 J4 W
2. 
   : w) p  s1 q5 u\" t0 E/ |
3. x=LL57(:,[3])/10;  }' ]' I+ V, g! F& A% c# `
   
4. 
   7 b' x, l. n8 L4 ?5 O3 Y% X
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');/ D5 s, I) a. |8 Q7 @9 q! Z
   
2. 
   0 k6 H$ k8 Q: L5 @
3. ymd57=fscanf(fpy,'%d',[3,1]);+ Y( q# m0 d9 B0 O
   
4. 
   7 J0 Z9 q- u# c' @4 l- B& @
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);, m' g' P; |- b! Q$ O8 k8 t
   
6. : \' y% r3 l! E. W
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   # O5 _! g5 V; n. y
2. 
   1 g! M4 o0 o1 l9 C0 w; O; B+ X& W
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
     K: D$ p# n) {# D
2. 
   2 L9 X$ j+ g* E# x7 w$ k5 o
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));
   6 E& D5 T6 A' p+ H8 Y
2. 9 K6 H\" G3 C( z* w8 @
   
3. zi=griddata(x,y,z,xi,yi,'cubic');5 l4 n- s4 Y# w$ R! i
   
4. ! K  T$ {) S2 P6 O
   

复制代码

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

1. hold on;0 @4 \& a% W* N# ?9 u: v# E+ z- _1 N
   
2. 7 G% O\" l  R9 T- f6 d8 \0 H
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);! p! M( ~- z5 q7 x3 W- q7 X
   
2. 0 x8 Y3 f# w' x; M; \$ S
   
3. hold off;

复制代码

这里用于在等值线上标注温度值，然后通过 hold off 结束图形的绘制。
请注意，代码中有一个拼写错误，i 应该是 ;，即 hold off;i 应该为 hold off;。
- v0 E: O8 F+ ~4 }3 c
8 \; @% c+ v4 Z: `% y# x

最后结果图如下：

. \4 A, v+ J, |5 q; i+ [

具体附件代码如下：

1 a7 b1 X( d: ^; e/ j% D
! \1 [, |, i  V" F9 o2 _6 k
0 o: H# G0 z; P+ S* F$ D8 Z
/ T' y0 `% P, A# Z# }4 U
8 t2 C) _* T% ]
' T! v' @+ k$ z% s+ u" g$ W