数学建模社区-数学中国

标题: matlab 画等温线 [打印本页]

---

作者: 2744557306    时间: 2023-11-12 11:13
标题: matlab 画等温线
解析代码：

1. clear; close all;

复制代码

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

1. fphn=fopen('hunan.txt','r');
   
2. 
   & T& I$ R5 \( P: X) r" z$ o
3. hnb=fgetl(fphn);
   
4. 
   
5. hnmap=fscanf(fphn,'%f %f',[2,59]);
   ' Q+ ?6 J  s& ^2 [- a
6. 
   
7. fclose(fphn);

复制代码

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

1. hnmap=hnmap';
   2 `1 S) P$ i6 F. L. {4 e( {
2. 
   " z! }$ O0 x: t2 y& f
3. xa=hnmap(:,[1]);
   4 |. R% k: p  a& i9 X! V
4. 
   ( I0 R+ V6 M6 K
5. ya=hnmap(:,[2]);
   3 q/ `7 n! q( C
6. 
   

复制代码

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

1. fp=fopen('LATLON57.txt','r');
   ( J" {4 w  \0 V6 J# @' I& @' W
2. 
   * R; I  I2 k/ C
3. LL57=fscanf(fp,'%d %f %f',[3,97]);
   
4. 
   + d- r' [" X8 X+ _
5. fclose(fp);
   
6. 
   0 j, ^/ D, y7 h0 q

复制代码

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

1. LL57=LL57';
   
2. 
   
3. x=LL57(:,[3])/10;
   : K, T# ?+ |( m+ E9 s
4. 
   / {* j2 J" a  k( I; v
5. y=LL57(:,[2])/10;

复制代码

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

1. fpy=fopen('etw00100.txt','r');
   
2. 
   
3. ymd57=fscanf(fpy,'%d',[3,1]);
   # @  D; s1 `- |) t. p2 F, q
4. 
   & X% b# [! T# J5 q7 Y# z. q
5. yu97=fscanf(fpy,'%d %f %f',[3,97]);
   
6. 
   
7. fclose(fpy);

复制代码

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

1. yu97=yu97';
   ; b) G; u: [' `
2. 
     Q# R, ?- M! B# J( Y4 J- y: P0 b
3. z=yu97(:,[2]);

复制代码

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

1. hold on;
   % I8 i! ^6 q3 N3 V, t) x8 U# Q
2. 
   5 K6 x9 H% n4 }/ U
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));
   " u- U1 V# a0 O
2. 
   
3. zi=griddata(x,y,z,xi,yi,'cubic');
   % w$ ~4 A  \0 z
4. 
   . L, D+ {5 z, Y

复制代码

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

1. hold on;
   3 C7 s0 ?, ?, x% U4 S1 y
2. 
   9 ]6 X* b) v2 [; K) g$ C5 ~
3. [c,h]=contour(xi,yi,zi,'b-');

复制代码

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

1. clabel(c,h);
   
2. 
   7 Q3 l7 c3 g" \0 a
3. hold off;

复制代码

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

% P( V# P- e9 v! G

最后结果图如下：

3 E" r+ O2 B2 F

# F$ s" r& _3 V% Z: H
) F; O) Q- F8 G; V! }2 g

具体附件代码如下：

3 @8 N& k) i5 r& h$ W1 T1 E! C3 J


7 C3 W8 K2 m/ n, A
6 y( w# v5 R, e' |9 j9 y8 X
0 ~3 r, t5 D# l0 p% L2 [

画等温线.rar

2023-11-12 11:13 上传

点击文件名下载附件

下载积分: 体力 -2 点

4.09 KB, 下载次数: 0, 下载积分: 体力 -2 点

售价: 2 点体力  [记录]  [购买]

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5