| 查看: 861 | 回复: 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 ] |
回复此楼» 猜你喜欢
为什么网上很多人说本周 12号出结果
已经有4人回复
-
基金中了
已经有15人回复
-
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
静等基金结果
已经有16人回复
wuqingming2003
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1585.2
- 帖子: 74
- 在线: 34.2小时
- 虫号: 447056
- 注册: 2007-10-30
- 专业: 极端条件下使用的金属材料
18楼2009-10-19 23:31:46
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1356.7
- 散金: 969
- 红花: 8
- 帖子: 863
- 在线: 226.5小时
- 虫号: 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
100