应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】如何做一个由三个点确定的切面的等高线？
北京石油化工学院2026年研究生招生接收调剂公告
181/2返回列表12下一页
查看: 740  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 wuqingming2003 的 50 个金币

wuqingming2003

木虫 (小有名气)

[交流] 【求助】如何做一个由三个点确定的切面的等高线？

本人是matlab新手，现得到一组数据（在附件中），是一个40*40*40的空间内的各个点的强度值，现在想做由(26,26,26)、(20,20,20)和(38,10,22)三个点确定的平面的等高线分布图，请问如何做？
请将程序写的详细一点，万分感谢！

[ Last edited by wuqingming2003 on 2009-10-14 at 09:41 ]
回复此楼

» 猜你喜欢
1楼 2009-10-14 09:25:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
wuqingming2003(金币+20,VIP+0):非常感谢，辛苦了！明天试试看，如果可以成功，50个金币全部赠送。 10-14 22:05
nono2009(金币+2,VIP+0):谢谢专家！麻烦再看一下此问题的进展。 10-15 13:37
引用回帖:
Originally posted by wuqingming2003 at 2009-10-14 09:25:
本人是matlab新手，现得到一组数据（在附件中），是一个40*40*40的空间内的各个点的强度值，现在想做由(26,26,26)、(20,20,20)和(38,10,22)三个点确定的平面的等高线分布图，请问如何做？
请将程序写的详细一点， ...

40*40*40的数组比较大，估计你的附件白贴了。

下面是一个例子，可供参考。


\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\

生成一个数据
>> [x y z]=meshgrid(1:40);
>> C=sin((x.^2+y.^2+z)/20);

定义制定的斜面


1，三个点的坐标；
>> a=[20 20 20];
>> b=[26 26 26];
>> c=[38 10 22];

2，计算x-z平面与预制平面的夹角

>> der1=[0 1 0]; %x-z平面法相方向
>> der2=cross(b-a,c-a);%预制平面法向方向
>> der3=cross(der1,der2);%两平面交线方向
>> the=dot(der2,der1)/sqrt(sumsqr(der1)*sumsqr(der2));%两平面夹角


3，画出这两个平面，看看是不是落在那三个点上？

>> [xx zz]=meshgrid(-10:50);yy=ones(size(xx));
>> mesh(xx,yy,zz),hold on;hsp = mesh(xx,yy,zz);
>> rotate(hsp,der3,the*180/pi,[1 1 1])%旋转平面，如附件一
>> scatter3(20,20,20)
>> scatter3(26,26,26)
>> scatter3(38,10,22)


4，得到旋转后平面的数据
>> xd = get(hsp,'XData');
>> yd = get(hsp,'YData');
>> zd = get(hsp,'ZData');


5，画一下且过那三个点的等值线，如附件2

>> contourslice(x,y,z,C,xd,yd,zd)
>> hold on;
>> contourslice(x,y,z,C,[],1,[])
>> scatter3(20,20,20)
scatter3(26,26,26)
scatter3(38,10,22)

/////////////////////////////////////////////////////////////////////////////
赞一下(16人) 回复此楼
2楼2009-10-14 19:16:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuqingming2003

木虫 (小有名气)
非常感谢你的回答，写的很详细，我明天试试看，再次感谢！
引用回帖:
Originally posted by hitzhang at 2009-10-14 19:16:



40*40*40的数组比较大，估计你的附件白贴了。

下面是一个例子，可供参考。


\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\

生成一个数据
>> [x y z]=meshgrid(1:40);
>> C=sin(( ...

赞一下 回复此楼
3楼2009-10-14 22:14:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuqingming2003

木虫 (小有名气)
nono2009(金币+0,VIP+0):建议PM或引用回复该帖，以便别人及时知道你的message！ 10-15 13:36
我用你说的方法试了一下，出现了点问题，我把程序写在下面，麻烦帮我看看哪里出错了，谢谢！
[path,fn]=uigetfile('*.txt','Open');
fp=fopen([fn,path],'r');
head=fscanf(fp,'%s',4);
data=fscanf(fp,'%f',[4,40*40*40]);
fclose(fp);
x=reshape(data(1,,[40,40,40]);
y=reshape(data(2,,[40,40,40]);
z=reshape(data(3,,[40,40,40]);
C=reshape(data(4,,[40,40,40]);
>> a=[20 20 20];
>> b=[26 26 26];
>> c=[38 10 22];
>> der1=[0 1 0]; %x-z平面法相方向
>> der2=cross(b-a,c-a);%预制平面法向方向
>> der3=cross(der1,der2);%两平面交线方向
>> the=dot(der2,der1)/sqrt(sumsqr(der1)*sumsqr(der2));%两平面夹角
>> [xx zz]=meshgrid(-10:50);yy=ones(size(xx));
>> mesh(xx,yy,zz),hold on;hsp = mesh(xx,yy,zz);
>> rotate(hsp,der3,the*180/pi,[1 1 1])%旋转平面
>> scatter3(20,20,20)
>> scatter3(26,26,26)
>> scatter3(38,10,22)
这样得到的结果中旋转后的平面并没有经过这三个点，如附件中的图所示
引用回帖:
Originally posted by hitzhang at 2009-10-14 19:16:



40*40*40的数组比较大，估计你的附件白贴了。

下面是一个例子，可供参考。


\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\

生成一个数据
>> [x y z]=meshgrid(1:40);
>> C=sin(( ...

赞一下(6人) 回复此楼
4楼2009-10-15 12:12:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kuhailangyu(金币+1,VIP+0):谢谢参与 10-15 16:01
wuqingming2003(金币+10,VIP+0):改成这个代码后，画的平面就经过三个点了 10-15 16:20
引用回帖:
Originally posted by wuqingming2003 at 2009-10-15 12:12:
我用你说的方法试了一下，出现了点问题，我把程序写在下面，麻烦帮我看看哪里出错了，谢谢！
[path,fn]=uigetfile('*.txt','Open');
fp=fopen([fn,path],'r');
head=fscanf(fp,'%s',4);
data=fscanf(fp,'%f', ...

呵呵，不好意思，代码贴错了。角度忘记求反余弦了，这样改

>> the=acos(dot(der2,der1)/sqrt(sumsqr(der1)*sumsqr(der2)));%两平面夹角
赞一下(16人) 回复此楼
5楼2009-10-15 15:42:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuqingming2003

木虫 (小有名气)
nono2009(金币+0,VIP+0):再提醒一次：建议PM或引用回复该帖，以便别人及时知道你的message！ 10-15 16:26
非常感谢！按照你的方法，点共面的问题已经解决了，可我将后续的程序输入后并没有得到等高线的图，得到的图是下面的样子

请问是怎么回事？不好意识，我是新手，又要麻烦你了！

我想得到类似第二张图那样的结果，请问该怎么弄？无需是彩色，只要二维面及其等高线即可，万分感谢！
引用回帖:
Originally posted by hitzhang at 2009-10-15 15:42:



呵呵，不好意思，代码贴错了。角度忘记求反余弦了，这样改

>> the=acos(dot(der2,der1)/sqrt(sumsqr(der1)*sumsqr(der2)));%两平面夹角

[ Last edited by wuqingming2003 on 2009-10-15 at 16:45 ]
赞一下(6人) 回复此楼
6楼2009-10-15 16:22:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)
引用回帖:
Originally posted by wuqingming2003 at 2009-10-15 16:22:
非常感谢！按照你的方法，点共面的问题已经解决了，可我将后续的程序输入后并没有得到等高线的图，得到的图是下面的样子

请问是怎么回事？不好意识，我是新手，又要麻烦你了！

我想得到类似第二张图那样的结 ...

没明白啥意思，旋转一下看看。
赞一下 回复此楼
7楼2009-10-15 18:38:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuqingming2003

木虫 (小有名气)
我把输出的结果旋转了一下，麻烦你看看，谢谢！
引用回帖:
Originally posted by hitzhang at 2009-10-15 18:38:

没明白啥意思，旋转一下看看。

赞一下 回复此楼
8楼2009-10-15 23:16:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)
引用回帖:
Originally posted by wuqingming2003 at 2009-10-15 23:16:
我把输出的结果旋转了一下，麻烦你看看，谢谢！

看啥？
回复此楼
9楼2009-10-15 23:26:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuqingming2003

木虫 (小有名气)
结果不是切面的等高线呀，请问是我哪里出错了，谢谢！
引用回帖:
Originally posted by hitzhang at 2009-10-15 23:26:

看啥？

赞一下 回复此楼
10楼2009-10-16 09:37:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuqingming2003 的主题更新
181/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学调剂求助 +6 LULONG1 2026-04-03 6/300 2026-04-03 23:13 by qzxyhcsy
[考研] 311求调剂 +20 zchqwer 2026-04-01 22/1100 2026-04-03 22:09 by lglzsd
[考研] 335求调剂 +7 沈清璃 2026-04-03 7/350 2026-04-03 18:55 by lijunpoly
[考研] 一志愿生物与医药，296分，求调剂 +8 66鹿 2026-04-03 9/450 2026-04-03 14:22 by 化学化工硕士招�
[考研] 求调剂 +3 usbdndj 2026-04-03 3/150 2026-04-03 14:10 by dxiaoxin
[考研] 0705理学294求调剂 +3 成果成果cg5 2026-04-03 3/150 2026-04-03 14:04 by simons1972
[考研] 303求调剂 +3 一色清羽 2026-04-02 4/200 2026-04-03 10:22 by 蓝云思雨
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 348求调剂 +6 吴彦祖24k 2026-04-02 6/300 2026-04-02 14:07 by 给你你注意休息
[考研] 085410 一志愿211 22408分数359求调剂 +3 123456789qw 2026-03-31 4/200 2026-04-02 00:06 by 义文wang
[考研] 379求调剂 +3 ?苦瓜不苦 2026-04-01 3/150 2026-04-01 20:09 by 暮云清寒
[考研] 301求调剂 +8 axibli 2026-04-01 8/400 2026-04-01 09:51 by 我的船我的海
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +5 双马尾痞老板2 2026-03-31 5/250 2026-04-01 09:04 by oooqiao
[考研] 材料科学与工程求调剂 +13 深V宿舍吧 2026-03-29 13/650 2026-03-31 19:50 by Dyhoer
[考研] 375求调剂 +7 雨夏整夜 2026-03-29 7/350 2026-03-31 18:52 by xhai2011
[考研] 江苏苏北高校诚邀调剂同学 +3 zzll406 2026-03-31 3/150 2026-03-31 16:54 by 及时行乐fan
[考研] 哈尔滨工业大学材料与化工专硕378求调剂 +3 塔比乌斯 2026-03-30 3/150 2026-03-30 22:55 by 无际的草原
[考研] 085701求调剂初试286分 +5 secret0328 2026-03-28 5/250 2026-03-30 12:54 by fangnagu
[考研] 2026年华南师范大学欢迎化学，化工，生物，生医工等专业优秀学子加入！ +3 llss0711 2026-03-28 6/300 2026-03-29 10:26 by llss0711
[考研] 药学105500求调剂 +3 Ssun。。 2026-03-28 3/150 2026-03-28 11:24 by lxf170613
信息提示
关闭
请填处理意见
关闭