matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   ; x/ ~$ |\" s; O0 j; Z& G
2. ! p9 I# Z% \6 `& b
   
3. hnb=fgetl(fphn);
   & p8 m) P/ H( i# I\" s5 _5 w
4. 
   9 A4 ~- T- `% k2 w4 T7 Z' c3 A
5. hnmap=fscanf(fphn,'%f %f',[2,59]);1 Q& H3 \5 n' X( n% b
   
6. 
   $ I+ v6 |' }5 m! Z1 ]: f- L7 O% Z
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';$ L$ P* w6 z* @
   
2. 
   2 l2 c9 o0 W1 x! x
3. xa=hnmap(:,[1]);) a5 }  N2 t9 y3 a. r\" _5 b
   
4. 
   & b5 j* M$ l' f5 o! t7 ~. u
5. ya=hnmap(:,[2]);' Z\" S$ n( s1 j( F( P0 b% m
   
6. + ?# F& U% N# Z9 I; X1 ^
   

复制代码

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

1. fp=fopen('LATLON57.txt','r');+ U% W! \# \5 h9 k  Y
   
2. ! l1 e/ u& ~8 l
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);& F5 l\" f' C) Q. Z5 Q; Z! \$ N
   
4. 
   ! L: n0 B  S% f* f
5. fclose(fp);! F5 ~5 K, A\" v# N# e' ]$ _\" J
   
6. 
   . F* j! n0 i' _: K! A: u8 p0 A+ ~6 r, ~

复制代码

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

1. LL57=LL57';\" {4 F9 k. |0 `& E9 u\" m# j
   
2. 6 E+ Z. i7 u) \( l4 _( a- w
   
3. x=LL57(:,[3])/10;
   ( m' h% f- J+ k! ~* t5 g1 Q/ d
4. 
   * B/ c: i\" Z4 u4 K+ r
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');
   0 y+ U- E' [\" I4 o
2. 9 Y6 V& o6 K8 V\" I  |
   
3. ymd57=fscanf(fpy,'%d',[3,1]);, w+ F+ s- `8 @  T! {
   
4. \" @8 k  n9 G! L, p
   
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);
   7 A% r! f/ ~1 w3 X4 k
6. 5 A3 z( K- t2 K: Q
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   ( [+ I8 m1 G6 T/ X- r
2. 
   9 c! v6 a3 k3 t  p
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   3 n/ D. E+ [4 `6 N7 |5 c
2. , n- y0 g. W; f' q
   
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));- Y; [8 S2 h) k3 d- P+ S+ y* z
   
2. 
   1 f- F2 m6 N! Z% B2 _
3. zi=griddata(x,y,z,xi,yi,'cubic');
   6 q\" C2 Z9 D. ?# T
4. / ~% g/ E$ M5 g- c' D, b5 ]
   

复制代码

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

1. hold on;2 \1 O% S2 H4 F) y( U- F8 V3 q
   
2. 
   ; C' |% }, v! s- C8 e3 a( S# I+ a
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);% X$ }( g9 ]9 a3 q
   
2. 
   + ]( k! O/ `$ g+ ^& i1 D2 Y; h, ^
3. hold off;

复制代码

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

7 l1 e6 ^. t( m! V2 ~

最后结果图如下：

& i$ s" X1 A; @7 ]) K9 m5 Q

具体附件代码如下：

# I9 L1 j) B- f) ]$ Z  K2 \
; D2 o5 R9 i* R3 V. b3 p4 z

! c: [; O% o, n9 A