【求助】如何做一个由三个点确定的切面的等高线？ - 仿真模拟 - MATLAB

11楼2009-10-16 12:56:47

wuqingming2003

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1585.2
帖子: 74
在线: 34.2小时
虫号: 447056
注册: 2007-10-30
专业: 极端条件下使用的金属材料

完成步骤3后，生成了图片，然后运行了步骤4和5，但是没有得到等高线的图，不知道是哪里出问题了，请指教，谢谢！

引用回帖:

Originally posted by hitzhang at 2009-10-16 12:56:

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

12楼2009-10-16 13:21:01

hitzhang

木虫 (正式写手)

仿真EPI: 1
应助: 0 (幼儿园)
贵宾: 2.15
金币: 1376.7
散金: 969
红花: 8
帖子: 862
在线: 226.4小时
虫号: 390575
注册: 2007-06-02
性别: GG
专业: 无机非金属类电介质与电解

引用回帖:

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

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

hold off

13楼2009-10-16 14:30:10

wuqingming2003

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1585.2
帖子: 74
在线: 34.2小时
虫号: 447056
注册: 2007-10-30
专业: 极端条件下使用的金属材料

>> [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

hitzhang

木虫 (正式写手)

仿真EPI: 1
应助: 0 (幼儿园)
贵宾: 2.15
金币: 1376.7
散金: 969
红花: 8
帖子: 862
在线: 226.4小时
虫号: 390575
注册: 2007-06-02
性别: GG
专业: 无机非金属类电介质与电解

★ ★
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

wuqingming2003

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1585.2
帖子: 74
在线: 34.2小时
虫号: 447056
注册: 2007-10-30
专业: 极端条件下使用的金属材料

换成figure后还是不行，和原来得到的图一样，就是只有三个点的一张白图。

引用回帖:

Originally posted by hitzhang at 2009-10-16 17:46:

把hold off换成figure试试。

16楼2009-10-17 18:39:05

hitzhang

木虫 (正式写手)

仿真EPI: 1
应助: 0 (幼儿园)
贵宾: 2.15
金币: 1376.7
散金: 969
红花: 8
帖子: 862
在线: 226.4小时
虫号: 390575
注册: 2007-06-02
性别: GG
专业: 无机非金属类电介质与电解

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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

wuqingming2003

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 1585.2
帖子: 74
在线: 34.2小时
虫号: 447056
注册: 2007-10-30
专业: 极端条件下使用的金属材料

★ ★ ★ ★ ★
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