| 查看: 774 | 回复: 17 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 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 ] |
回复此楼» 猜你喜欢
2026山东省优青
已经有4人回复
-
求助大佬sci投稿哪个好中
已经有5人回复
-
献血感触
已经有4人回复
-
评审有感
已经有16人回复
-
反应很差，大量原料没有反应
已经有3人回复
-
Sci. Bull. 悲剧经验
已经有4人回复
-
找博士生导师
已经有7人回复
-
上海大学实验技术岗位非升即走
已经有11人回复
-
26/27申博自荐-锂/钠电池方向
已经有4人回复
-
同样的基金本子，换个专家直接从C变A！
已经有3人回复
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1376.7
- 散金: 969
- 红花: 8
- 帖子: 862
- 在线: 226.4小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
wuqingming2003(金币+20,VIP+0):非常感谢，辛苦了！明天试试看，如果可以成功，50个金币全部赠送。 10-14 22:05
nono2009(金币+2,VIP+0):谢谢专家！麻烦再看一下此问题的进展。 10-15 13:37
wuqingming2003(金币+20,VIP+0):非常感谢，辛苦了！明天试试看，如果可以成功，50个金币全部赠送。 10-14 22:05
nono2009(金币+2,VIP+0):谢谢专家！麻烦再看一下此问题的进展。 10-15 13:37
|
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) ///////////////////////////////////////////////////////////////////////////// |
2楼2009-10-14 19:16:43
wuqingming2003
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1585.2
- 帖子: 74
- 在线: 34.2小时
- 虫号: 447056
- 注册: 2007-10-30
- 专业: 极端条件下使用的金属材料
3楼2009-10-14 22:14:29
wuqingming2003
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1585.2
- 帖子: 74
- 在线: 34.2小时
- 虫号: 447056
- 注册: 2007-10-30
- 专业: 极端条件下使用的金属材料
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) 这样得到的结果中旋转后的平面并没有经过这三个点，如附件中的图所示 |
4楼2009-10-15 12:12:23
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1376.7
- 散金: 969
- 红花: 8
- 帖子: 862
- 在线: 226.4小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
5楼2009-10-15 15:42:33
wuqingming2003
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1585.2
- 帖子: 74
- 在线: 34.2小时
- 虫号: 447056
- 注册: 2007-10-30
- 专业: 极端条件下使用的金属材料
6楼2009-10-15 16:22:53
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1376.7
- 散金: 969
- 红花: 8
- 帖子: 862
- 在线: 226.4小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
7楼2009-10-15 18:38:03
wuqingming2003
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1585.2
- 帖子: 74
- 在线: 34.2小时
- 虫号: 447056
- 注册: 2007-10-30
- 专业: 极端条件下使用的金属材料
8楼2009-10-15 23:16:29
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1376.7
- 散金: 969
- 红花: 8
- 帖子: 862
- 在线: 226.4小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
9楼2009-10-15 23:26:19
wuqingming2003
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1585.2
- 帖子: 74
- 在线: 34.2小时
- 虫号: 447056
- 注册: 2007-10-30
- 专业: 极端条件下使用的金属材料
10楼2009-10-16 09:37:14
5