matlab广州等温线地形图,MATLAB 画等温线.doc
测量到不同坐标点的高度值,如何用matlab画三维图
附上部分数据:
A=[-210.6627 -33391.1192 5.0273
-221.3052 -33387.7415 4.5969
-210.9391 -33393.0068 5.5647
-221.8901 -33390.7396 5.0077
-211.384 -33394.7093 5.6505
-222.6117 -33392.778 5.0554
-212.7074 -33397.5459 5.7381
-225.8973 -33397.5869 5.5587];
解:代码在matlab2009a版以上均可运行。
A=[-210.6627 -33391.1192 5.0273
-221.3052 -33387.7415 4.5969
-210.9391 -33393.0068 5.5647
-221.8901 -33390.7396 5.0077
-211.384 -33394.7093 5.6505
-222.6117 -33392.778 5.0554
-212.7074 -33397.5459 5.7381
-225.8973 -33397.5869 5.5587];
xData = A(:,1);
yData = A(:,2);
zData = A(:,3);
fitresult = fit( [xData, yData], zData, linearinterp);
figure( Name, 三维图 );
plot( fitresult, [xData, yData], zData );
xlabel( x );
ylabel( y );
zlabel( z );
grid on
view( -53, 50 );
试一试:
A=[-210.6627 -33391.1192 5.0273
-221.3052 -33387.7415 4.5969
-210.9391 -33393.0068 5.5647
-221.8901 -33390.7396 5.0077
-211.384 -33394.7093 5.6505
-222.6117 -33392.778 5.0554
-212.7074 -33397.5459 5.7381
-225.8973 -33397.5869 5.5587];
xData = A(:,1);
yData = A(:,2);
zData = A(:,3);
minx = min(xData);
maxx = max(xData);
miny = min(yData);
maxy = max(yData);
tx = linspace(minx,maxx);
ty = linspace(miny,maxy);
[XI,YI] = meshgrid(tx,ty);
ZI = griddata(xData,yData,zData,XI,YI);
mesh(XI,YI,ZI), hold
plot3(xData,yData,zData,o), hold off
xlabel( x );
ylabel( y );
zlabel( z );
grid on
view( -53, 50 );
.................
x=[7392 7392 7392 7392 5608 5608 5608 5608 8000 5000 6500 6500 6500 6500 6500 6500 6500 6500 6500 6500];
y=[476 404 476 404 476 404 476 404 440 440 440 440 500 380 440 440 440 440 440 440];
z=[17.06 16.93 17.8 17.4 17.17 16.89 17.28 17.26 17.04 16.94 16.92 17.3 17.37 17.17 17.34 17.11 17.22 17.24 17.11 17.2];
minx = min(x);
maxx = max(x);
miny = min(y);
maxy = max(y);
[X,Y,Z]=griddata(x,y,z,linspace(minx,maxx),linspace(miny,maxy),v4); %插值
figure,surf(X,Y,Z) %三维曲面
...
%% Initialization.
% Convert all inputs to column vectors.
x = x(:);
y = y(:);
z = z(:);
%% Fit: fit 1.
ft = linearinterp;
opts = fitoptions( ft );
opts.Weights = zeros(1,0);
opts.Normalize = on;
[fitresult, gof] = fit( [x, y], z, ft, opts );
% Plot fit with data.
figure( Name, untitled fit 1 );
h = plot( fitresult, [x, y], z );
% legend( h, fit 1, z vs. x, y, Location, NorthEast );
% Label axes
xlabel( x );
ylabel( y );
zlabel( z );
grid on
.......................
x=[7392 7392 7392 7392 5608 5608 5608 5608 8000 5000 6500 6500 6500 6500 6500 6500 6500 6500 6500 6500];
y=[476 404 476 404 476 404 476 404 440 440 440 440 500 380 440 440 440 440 440 440];
z=[17.06 16.93 17.8 17.4 17.17 16.89 17.28 17.26 17.04 16.94 16.92 17.3 17.37 17.17 17.34 17.11 17.22 17.24 17.11 17.2];
% Convert all inputs to column vectors.
x = x(:);
y = y(:);
z = z(:);
%% Fit: fit 1.
ft = fittype( poly22 );
opts = fitoptions( ft );
opts.Lower = [-Inf -Inf -Inf -Inf -Inf -Inf];
opts.Upper = [Inf Inf Inf Inf Inf Inf];
opts.Weights = zeros(1,0);
[fitresult, gof] = fit( [x, y], z, ft, opts );
fitresult
% Plot fit with data.
figure( Name, z );
h = plot( fitresult, [x, y], z );
% legend( h, z fit 1, z vs. x, y, Location, NorthEast );
% Label axes
xlabel( x );
ylabel( y );
zlabel( z );
grid on
view( -57.5, 42 );
.................
A=[
18467 17001 308
20177 17642 276
26852 16114 225
23785 17643 194
27823 14737 189
18906 16346 173
25981 18051 173
22046 17634 171
23238 6502 169
27696 11621 169
27700 11609 165
16607 17365 155
23146 15382 153
14405 18032 152
27232 14482 150
27380 18202 136
26591 13715 126
14074 16516 124
15801 3966 115
21684 13101 114
15255 5110 110
24580 13319 107
23325 16701 105
24065 7353 104
25021 16290 104
15952 18397 103
27346 13331 100
24685 14278 98
17414 15476 97
12778 5799 93
17044 10691 93
17981 18449 93
20983 15862 93
22605 14301 93
5451 2757 92
19041 15769 90
24003 15286 90
7100 2449 89
18413 11721 88
21475 8540 85
19007 11488 84
17008 4775 82
17079 5894 81
13855 3345 79
13920 5354 79
22193 12185 79
16569 6055 78
18993 12371 78
22965 13535 78
15517 17034 77
24631 9422 76
4153 2299 73
11702 4480 71
24153 12450 71
15007 5535 70
19569 7348 70
15087 3512 69
22846 9149 69
25461 9834 68
16823 4207 67
21766 12348 67
12625 16259 66
];
x=A(:,1);
y=A(:,2);
z=A(:,3);
minx = min(x);
maxx = max(x);
miny = min(y);
maxy = max(y);
[X,Y,Z]=griddata(x,y,z,linspace(minx,maxx),linspace(miny,maxy),v4); %插值
figure,surf(X,Y,Z) %三维曲面
...
%% Initialization.
% Convert all inputs to column vectors.
x = x(:);
y = y(:);
z = z(:);
%% Fit: untitled fit 1.
ft = linearinterp;
opts = fitoptions( ft );
opts.Weights = zeros(1,0);
opts.Normalize = on;
[fitresult, gof] = fit( [x, y], z, ft, opts );
% Plot fit with data.
figure( Name, fit 1 );
h = plot( fitresult, [x, y], z );
% legend( h, fit 1, z vs. x, y, Location, NorthEast );
% Label axes
xlabel( x );
ylabel( y );
zlabel( z );
grid on
展开阅读全文
matlab广州等温线地形图,MATLAB 画等温线.doc相关推荐
- matlab脑电地形图怎么画,脑电地形图绘制.ppt
脑电地形图绘制 脑电地形图绘制 报告内容 数据可视化 脑电地形图的意义 脑电地形图绘制的关键 离线地形图 实时地形图 Ⅰ 数据可视化 概念:借助图形化手段,将数据以图形图像形式表示,并利用数据分析和开 ...
- matlab画三维图电磁场,利用matlab实现矩形波导电磁场分布图的绘制.doc
利用matlab实现矩形波导电磁场分布图的绘制.doc 利用MATLAB实现矩形波导电磁场分布图的绘制(附源程序)通过MATLAB计算并绘出任意时刻金属矩形波导的主模TE10模的电磁场分布图.波导尺寸 ...
- matlab 一条水平线,matlab画一条水平线
凸轮廓线的MATLAB画法_工学_高等教育_教育专区.基于matlab的凸轮画法实例 凸轮廓线的 MATLAB 画法 1 凸轮轮廓方程 X ? OE ? EF ? E * Cos( J ) ? ( . ...
- 频域串联滞后校正matlab,自动控制原理实验七基于MATLAB控制系统频域法串联校正设计.doc...
自动控制原理实验七基于MATLAB控制系统频域法串联校正设计.doc 实验七 基于MATLAB控制系统频域法串联校正设计 一.实验目的 (1)对给定系统设计满足频域或时域指标的串联校正装置: (2)掌 ...
- 摄像机投影成像 matlab,使用matlab仿真三维物点的透视投影成像.doc
使用matlab仿真三维物点的透视投影成像.doc 使用matlab仿真三维物点的透视投影成像 ⒈仿真的目的和要求: 理解摄像机透视投影模型中的每一个参数的意思,然后在matlab中仿真三维物点根据摄 ...
- 热力学分布用matlab,matlab在热物理学中的应用.doc
matlab在热物理学中的应用.doc 1MATLAB在热物理学中的应用摘要本文阐述了基于MATLAB的数值计算.可视化图形处理.开放式以及可扩充体系结构的特点,并介绍了高性能语言MATLAB在大学物 ...
- matlab进行fm调制与解调,基于matlab的fm信号调制与解调.doc
基于matlab的fm信号调制与解调.doc 调制就是将基带信号的频谱搬移到信道通带中或者其中的某个频段上的过程,而解调是将信道中来的频带信号恢复为基带信号的反过程.调制的目的是把要传输的模拟信号或数 ...
- matlab 汽车 流场,matlab画流场图
基于 Matlab 分布式工具箱的流场计算及其可视化 蔡群;周美莲;段杰峰;李青... 基于 MATLAB 和 CFD 数据库的流场可视化的实现 [J], 晏畅 5.基于 VB 与 MATLAB 混合 ...
- matlab水印剪切攻击程序,可以运行的水印matlab程序(嵌入,提取,攻击测试等).doc
可以运行的水印matlab程序(嵌入,提取,攻击测试等).doc IMREAD YUANSHITUXIANGBMP SUBPLOT2,2,1IMSHOWTITLE 原始图像 WATERIMREAD W ...
最新文章
- pinpoint zipink skywalking
- 360首席安全官谭晓生宣布离职
- LeetCode 20 有效的括号
- java中获取绝对值的方法_Java判断绝对值的方法总结
- android 调整cpu频率,[IMX6DL][Android4.4] CPU频率调节模式以及降频方法
- C++远征之封装篇——对象数组,对象成员
- textedit实时显示位置_加什么地形就看什么等高线!等高线实时预览就是这么爽...
- 运用ENVI bandmath对多波段进行计算
- JAVA 反射(JAVA面试题)
- Word XP 中目录的编制方法,Word XP 功能键使用大全
- 科学计算机求年金,科普如何科学使用普通计算器(CPA考试小技巧)
- 河北省教育厅计算机制度,河北省中小学计算机室规章制度.doc.docx
- python+nodejs+vue自驾游旅游网站系统
- Meta研发触觉手套,有人想在元宇宙撸猫,有人产生大胆的想法
- 毁人不倦 - 许嵩 - 苏格拉没有底
- 遇到无法打出的字怎么办?
- iPhone开发Blog
- ReFusion: 3D Reconstruction in Dynamic Environments for RGB-D Cameras Exploiting Residuals 论文阅读
- 教育系统APP(三)
- 网络安全渗透测试的常用工具