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

24小时热门版块排行榜    

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

hitzhang

木虫 (正式写手)
★ ★ ★ ★ ★
wuguocheng(金币+5,VIP+0): 谢谢专家指点 10-16 16:39
引用回帖:
Originally posted by wuqingming2003 at 2009-10-16 09:37:
结果不是切面的等高线呀，请问是我哪里出错了，谢谢！

我给你的步骤4，5你没运行吧
赞一下(10人) 回复此楼
11楼2009-10-16 12:56:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuqingming2003

木虫 (小有名气)
完成步骤3后，生成了图片，然后运行了步骤4和5，但是没有得到等高线的图，不知道是哪里出问题了，请指教，谢谢！
引用回帖:
Originally posted by hitzhang at 2009-10-16 12:56:

我给你的步骤4，5你没运行吧

赞一下 回复此楼
12楼2009-10-16 13:21:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)
引用回帖:
Originally posted by wuqingming2003 at 2009-10-16 13:21:
完成步骤3后，生成了图片，然后运行了步骤4和5，但是没有得到等高线的图，不知道是哪里出问题了，请指教，谢谢！

哦，那是假象，你在运行步骤4，5之前执行一下

hold off
赞一下 回复此楼
13楼2009-10-16 14:30:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuqingming2003

木虫 (小有名气)
>> [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]);
w=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=acos(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)
>> hold off
>> xd = get(hsp,'XData');
yd = get(hsp,'YData');
zd = get(hsp,'ZData');
>> contourslice(x,y,z,w,xd,yd,zd)
hold on;
contourslice(x,y,z,w,[],1,[])
scatter3(20,20,20)
scatter3(26,26,26)
scatter3(38,10,22)

我运行到这个程序，得到的结果还是没有等高线，图形如下，请问是哪里出错了？谢谢！真是麻烦你了！
引用回帖:
Originally posted by hitzhang at 2009-10-16 14:30:

哦，那是假象，你在运行步骤4，5之前执行一下

hold off

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

hitzhang

木虫 (正式写手)
★ ★
wuguocheng(金币+2,VIP+0): 谢谢帮助 10-17 07:28
引用回帖:
Originally posted by wuqingming2003 at 2009-10-16 15:18:
>> [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 ...

把hold off换成figure试试。
赞一下(10人) 回复此楼
15楼2009-10-16 17:46:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuqingming2003

木虫 (小有名气)
换成figure后还是不行，和原来得到的图一样，就是只有三个点的一张白图。
引用回帖:
Originally posted by hitzhang at 2009-10-16 17:46:

把hold off换成figure试试。

赞一下 回复此楼
16楼2009-10-17 18:39:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
nono2009(金币+2,VIP+0):谢谢专家！ 10-18 10:16
wuqingming2003(金币+20,VIP+0):感谢几天来尽心尽力的帮助！ 10-19 23:32
引用回帖:
Originally posted by wuqingming2003 at 2009-10-17 18:39:
换成figure后还是不行，和原来得到的图一样，就是只有三个点的一张白图。



试试这段程序，如果还不行，我也没辙了。
>> [x y z]=meshgrid(1:40);
>> C=sin((x.^2+y.^2+z)/20);
>> 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=acos(dot(der2,der1)/sqrt(sumsqr(der1)*sumsqr(der2)));%两平面夹角
>> [xx zz]=meshgrid(-10:50);yy=ones(size(xx));
>> hsp = mesh(xx,yy,zz);
>> rotate(hsp,der3,the*180/pi,[1 1 1])
>> xd = get(hsp,'XData');
>> yd = get(hsp,'YData');
>> zd = get(hsp,'ZData');
>> delete(hsp)
>> contourslice(x,y,z,C,xd,yd,zd)
>> hold on;
>> scatter3(20,20,20)
>> scatter3(26,26,26)
>> scatter3(38,10,22)
>>

[ Last edited by hitzhang on 2009-10-18 at 10:12 ]
赞一下(24人) 回复此楼
17楼2009-10-18 10:11:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuqingming2003

木虫 (小有名气)
★ ★ ★ ★ ★
nono2009(金币+5,VIP+0):很好的交流。欢迎常来仿真模拟版。 10-20 06:39
感谢你这几天的帮助！虽然我还是没有得到想要的结果，估计是和我的数据的C的提取有关，我自己试着去解决吧，就当是学习了，再次感谢！
引用回帖:
Originally posted by hitzhang at 2009-10-18 10:11:


试试这段程序，如果还不行，我也没辙了。
>> [x y z]=meshgrid(1:40);
>> C=sin((x.^2+y.^2+z)/20);
>> a=[20 20 20];
>> b=[26 26 26];
>> c= ...

赞一下(10人) 回复此楼
18楼2009-10-19 23:31:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuqingming2003 的主题更新
182/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode +8 documentary 2026-08-10 8/400 2026-08-10 19:20 by loufangrui
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 帮忙看看fileCode +3 wwncly 2026-08-10 5/250 2026-08-10 19:06 by 2000zf36392
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 有没有H口的？有收到消息的吗？ +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
信息提示
关闭
请填处理意见
关闭