matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   1 l4 u# _* P/ `; k8 ?9 K6 d
2. 
   3 J. K. R* y0 R7 T# y1 f1 U
3. hnb=fgetl(fphn);& g2 U! F3 B4 S6 s8 I0 z
   
4. 
   - F. y1 H3 t* p2 u1 {: c7 Y
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
   9 o# [! D! N, @; G( \
6. : K! z/ H: K6 |
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';8 [1 b; Q6 f1 B+ [
   
2. 
   \" E0 y3 L5 Y% @) ^/ Z0 u5 }5 o
3. xa=hnmap(:,[1]);
   , \9 }\" S# [0 y# q
4. ! p+ `: A! U) ]5 W' g
   
5. ya=hnmap(:,[2]);
   ) p- s) \2 ?3 r$ N, z) n' [
6. + `2 X2 h* U. S9 O\" L: ]5 t
   

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   9 H% G, d  W. e+ f
2. 1 k1 V4 X/ L* |: A\" s( D
   
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   : a5 Y7 \4 ~! }7 P( X/ W
4. ' o/ b/ g. I: m
   
5. fclose(fp);
   ) j4 m) t2 j2 D* c' e' O: g
6. ( M, ^8 n$ V\" G: [! i- q6 m; h
   

复制代码

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

1. LL57=LL57';
   2 X\" d  R- l\" o0 a$ i
2. \" S/ I7 E5 ?3 \$ T' Q9 h2 k) a\" P
   
3. x=LL57(:,[3])/10;5 v\" `9 f$ N( Q( Q* e
   
4. 
   - v( P' [; h- d\" ?, a
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');0 {( D5 ]8 B\" w# W5 O
   
2. 
   4 L7 V' K2 e5 w9 _( V, [5 _
3. ymd57=fscanf(fpy,'%d',[3,1]);% ]$ O\" ~# P' }
   
4. 
   & s  s2 S( n& N$ x' J
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);4 k+ p; h+ B9 }: Z! z7 K
   
6. 0 g% o\" A0 w5 o8 W3 H9 m
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   . R* d8 H: q2 o; p+ q. y1 \
2. 
   2 f% M7 ^, D/ f2 h3 ^$ W
3. z=yu97(:,[2]);

复制代码

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

1. hold on;' ~0 A. a& v# \/ ^3 d
   
2. 6 `* {+ K. d0 D* _- U! 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));
   , k6 P9 `3 U' E$ Z
2. 
   $ m  d0 j% c, u
3. zi=griddata(x,y,z,xi,yi,'cubic');+ b4 q5 c  q& g9 F* o. B6 _# ^
   
4. 6 A! [) b( v& u8 b- Z\" L) Q
   

复制代码

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

1. hold on;
     Q5 h+ q\" u' x: }1 @/ e4 g
2. 
   . B* a\" w1 F4 R
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   2 ?0 w' ]; W5 ?) J' r# X- L2 _
2. % O- S/ [- L0 _6 {! r$ K
   
3. hold off;

复制代码

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

% Y) Y" A4 w* x8 T

最后结果图如下：

. D9 D0 ~5 B* ]6 k8 g. H

9 p7 R4 j/ \. J

具体附件代码如下：

; A# q7 d  b$ F+ U6 l


  J' _+ x$ w+ O1 U* A
" J; ]) m6 t# X/ `6 D0 k3 o6 `! N" K