matlab 画等温线

2744557306

解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');2 ^: g. t! }7 x( s' q
   
2. 
   1 M( ]- c' Z; |# @- r
3. hnb=fgetl(fphn);
   # l* e7 A, b# q* Q0 p9 J3 w
4. 4 K/ l( D, x7 w) W6 Q3 ^# V
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
   0 U& G- J) I( u1 c, M2 E\" C0 X0 D% r- u
6. & X' D- w! F- z
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';
   ( y- p\" ^\" X. W: T$ b
2. 
   5 t0 e6 r5 ?9 X' j- b& D
3. xa=hnmap(:,[1]);0 t# N6 z\" ^* Y2 O
   
4. ( n: a4 E% A( o
   
5. ya=hnmap(:,[2]);. J+ u0 G: h2 j\" D$ ?% Q, H% U) G
   
6. 
   3 D\" q% G6 f( r3 Q7 B% A( y* a\" U

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   2 S4 ^# Z/ i; e
2. 
   8 K7 L, a) j* d1 b( Y- N
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   6 o& p* B, }9 b\" x6 m; b
4. ! w% {- y' L6 q9 t
   
5. fclose(fp);  `\" F) [1 f# `3 c# J
   
6. 
   ( Z  d1 D3 t, a- S\" s

复制代码

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

1. LL57=LL57';, M: K# t8 `, a7 \  k
   
2. 
   ( r  g) P, Q) v4 u$ B! h
3. x=LL57(:,[3])/10;
   $ i& V\" m3 u+ {; Y
4. 
   % {: n! n9 f* p
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');7 s; {( ~& _6 A0 l$ Y* M7 r7 z/ e
   
2. 1 o% q5 I% l+ ]% [$ n' U2 y% y
   
3. ymd57=fscanf(fpy,'%d',[3,1]);
   4 h# y' d8 {+ J. G) e0 H7 M
4. 1 p  X. T8 {9 g& w+ T* V
   
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);
   ( j3 V. j3 Z* L
6. % A) X2 x$ ^* b# i( U0 U0 M
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   ) X2 V% |% j4 D, [/ E+ U: R
2. 0 K- C- H, {5 Z( ^/ ^1 F7 R# J
   
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   \" R+ V; ~) ^3 ~+ J0 S
2. 
   \" y# Q& D) M( J\" x
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));! M\" v3 k- y5 n! ]
   
2. 
   , k0 w  ]) v/ Y4 T
3. zi=griddata(x,y,z,xi,yi,'cubic');$ F) n, {9 s3 C7 Q+ J$ P) q
   
4. 
   ( X- j' n6 G4 g0 D8 c( v\" K0 x

复制代码

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

1. hold on;
     V% \- G5 q2 w6 E
2. ' `. Y7 x, m+ T) v/ K* W& [
   
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   # a4 x; t9 n4 i
2. 
   : f% |( T0 J/ o; ~% j
3. hold off;

复制代码

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


* P: ^! |- S0 l, J7 {0 J9 g

最后结果图如下：

, l8 i3 M, X. X! o) u

! ^$ u5 g1 d5 ~2 u5 x! X6 I
8 \8 O  p1 T9 k* {8 b4 E

具体附件代码如下：

" m' r2 ?1 o$ o% r' G5 a" U8 M

+ h/ X6 \! q, R& r7 K% w