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 ...
最新文章
- Unity下的ECS框架 Entitas简介
- 实战SSM_O2O商铺_11【商铺注册】Controller层的实现
- [MATLAB调试笔记]Possion方程求网格点电场强度
- 谷歌gcp 远程计算机_引导性GCP:带有Google Cloud Pub / Sub的Spring Cloud Stream
- python jieba词频统计英文文本_python实战,中文自然语言处理,应用jieba库来统计文本词频...
- 计算机等级考试java题型_计算机等级考试报考建议
- 步骤一:入门linux基础/01Linux简介和安装/003Linux系统的多面性
- hdu-1088 Write a simple HTML Browser
- SQL数据库实战之查找当前薪水详情以及部门编号dept_no
- mac memcached_如何在Mac OS上安装Memcached Server
- 【Android】【移动应用开发】APP案列
- foobar 更换皮肤
- Java判断手机号格式是否正确
- 四、Mosquitto 高级应用之用户配置
- 第二人生的源码分析(二十六)底层网络协议
- 如何查区块链项目的服务器地址,怎么查区块链服务器地址
- 启动(程序还没执行生成dump文件的代码)就崩溃的处理流程
- 如何将镜像文件上传到服务器,通过把docker镜像保存为文件载入到别的服务器
- 警告:关于电磁辐射对孕妇的危害。
- 《孙子兵法战略运筹之始计篇》
热门文章
- Java 1.3 控制流程(条件语句、循环结构)
- algorithm头文件下的next_permutation()
- 《南溪的目标检测学习笔记》——夏侯南溪的CNN调参笔记,加油
- 面向对象语言编程的心得
- http和websocket共用同一端口
- masaic 数据增强代码
- LeetCode 14. 最长公共前缀 (单指针水平扫描)
- 无烟会议室:CMMI vs. Scrum vs. XP(QCon 2010 感受)
- ActionScript 与后台交互 RemoteObject
- java 传址或传值