matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   ) t: m: W1 K2 Y2 Y# Y6 t
2. 9 j$ P, V  X+ o0 `0 t. I5 \9 i
   
3. hnb=fgetl(fphn);6 ]% X7 s6 J& _/ S* g/ _% H+ U
   
4. 
   0 j4 f( O8 S: Z/ s
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
   + l\" u- @! Q3 \6 o
6. 
   \" G, e# ?  b2 h; W  S
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';' c7 M3 O' w4 E
   
2. ( d/ j5 d2 S  U7 O
   
3. xa=hnmap(:,[1]);+ t9 a! M# Q2 e5 S# b! h% P# C$ q
   
4. 
   ' r9 _( p- T2 S/ I4 I2 \
5. ya=hnmap(:,[2]);
   : w2 m- R1 j- y, W7 l& E$ B
6. 
   - C1 X& W  R, n. ^

复制代码

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

1. fp=fopen('LATLON57.txt','r');/ {7 g  U3 I& J2 L5 c: _
   
2. 
   * `$ ^' H8 f0 r
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   ( y; N; ^& w) X* d0 J7 c
4. \" |- v: |0 i! v4 y
   
5. fclose(fp);
   # h6 a/ L2 U- S. G
6. 0 N# G8 c8 ^9 R, u) ~0 b
   

复制代码

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

1. LL57=LL57';
   7 K+ ]0 |' }& h& q  y5 X- I1 V
2. 
   % p. v- q+ L# {' H% Y% y# [, q
3. x=LL57(:,[3])/10;
   - x5 d; M3 W5 b\" `7 {
4. 
   : R0 y  U$ j: O+ M/ N$ T4 ]
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');
   / }  r5 _0 T% X; N0 k
2. 
   # G+ }7 d1 ?! [5 q\" d' u, s
3. ymd57=fscanf(fpy,'%d',[3,1]);
   ) r; k$ k. j8 b( O* u; w: L
4. ) ?1 h& y! e' s8 S1 u) j* V* r# \
   
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);( f* R- r/ k; m% ^! `. M
   
6. 
   0 z2 d  r9 \1 E4 L& L. d4 s6 j
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   . ~- o4 n\" r. k  z2 c
2. 
   \" o' a# n& u# v0 u1 Z( x  x
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   4 J& Q' i* Q' J\" z\" `; x
2. , W6 _( I$ f4 J
   
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));) \% g& R6 j7 H, N& B\" l3 g; I: Z% i
   
2. 9 g% _% I& f* X7 M
   
3. zi=griddata(x,y,z,xi,yi,'cubic');2 e* S# `( n0 [$ N  @: g
   
4. 
   3 ?$ P7 @% O$ p; z) p

复制代码

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

1. hold on;* J7 i1 K# H8 Q0 D, o- F0 j/ N
   
2. 
   $ d% O0 {1 h7 {$ E2 y( u: I  ^
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);2 K9 L6 \( Q& h3 O8 }% A* l
   
2. : Q& s- v+ P, ~2 w7 @3 m1 [
   
3. hold off;

复制代码

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


: k$ ~$ E: f" l

最后结果图如下：

/ w3 H5 @5 s# D0 P* _7 T, J

具体附件代码如下：

1 C- L2 s# I/ G

  D# K- e8 l3 G: R7 {( P' z1 O& X% x2 f
3 @0 V3 M5 m. [0 K0 N# X
* J, `/ f6 t+ Z( l( j9 I% \* x. u