怎样画这些数据的三维曲面图？？

董雪娇

我是个新手，做有限元分析后产生了这些数据，前两列是节点的x,y坐标，最后一列是节点z向变形值（姑且代表z坐标值），现在想要将这些节点的xyz坐标绘制成三维曲面图，迫切需要高手的指导，自己也尝试了一些代码编程，但是都没有实现，谢谢大家啦！！！

0.15	0.75	0.14982E-04
0.15	0.3	0.16309E-04
0.15	0.345	0.16427E-04
0.15	0.39	0.16500E-04
0.15	0.435	0.16510E-04
0.15	0.48	0.16475E-04
0.15	0.525	0.16382E-04
0.15	0.57	0.16220E-04
0.15	0.615	0.16013E-04
0.15	0.66	0.15733E-04
0.15	0.705	0.15395E-04
0.501	0.3	0.60191E-05
0.45713	0.3	0.73848E-05
0.41325	0.3	0.87039E-05
0.36938	0.3	0.99648E-05
0.3255	0.3	0.11188E-04
0.28163	0.3	0.12428E-04
0.23775	0.3	0.13703E-04
0.19387	0.3	0.15005E-04
0.501	0.75	0.54980E-05
0.45713	0.75	0.67822E-05
0.41325	0.75	0.79920E-05
0.36938	0.75	0.91438E-05
0.3255	0.75	0.10267E-04
0.28163	0.75	0.11407E-04
0.23775	0.75	0.12586E-04
0.19387	0.75	0.13784E-04
0.501	0.345	0.59298E-05
0.501	0.39	0.58660E-05
0.501	0.435	0.58052E-05
0.501	0.48	0.57601E-05
0.501	0.525	0.57076E-05
0.501	0.57	0.56479E-05
0.501	0.615	0.55895E-05
0.501	0.66	0.55370E-05
0.501	0.705	0.54897E-05
0.21489	0.71568	0.13465E-04
0.25953	0.71208	0.12238E-04
0.30356	0.71135	0.11050E-04
0.34744	0.71135	0.98781E-05
0.39205	0.71129	0.86981E-05
0.43611	0.71568	0.74500E-05
0.21489	0.33432	0.14445E-04
0.25819	0.33963	0.13163E-04
0.3034	0.34065	0.11845E-04
0.3479	0.33994	0.10594E-04
0.39131	0.33865	0.93510E-05
0.43611	0.33432	0.80006E-05
0.46545	0.68571	0.66081E-05
0.46215	0.64016	0.67749E-05
0.46211	0.59402	0.68567E-05
0.46207	0.54849	0.69376E-05
0.46224	0.50302	0.69911E-05
0.46217	0.45714	0.70507E-05
0.46219	0.41056	0.71027E-05
0.46513	0.36507	0.70640E-05
0.18548	0.6864	0.14507E-04
0.18886	0.64011	0.14707E-04
0.18918	0.59363	0.14932E-04
0.18906	0.54797	0.15103E-04
0.18832	0.5025	0.15249E-04
0.18804	0.45698	0.15326E-04
0.18857	0.41032	0.15314E-04
0.18559	0.36503	0.15378E-04
0.23148	0.66888	0.13283E-04
0.28104	0.67288	0.11844E-04
0.32609	0.67254	0.10604E-04
0.37063	0.67243	0.94055E-05
0.4199	0.66792	0.80214E-05
0.23075	0.38238	0.14014E-04
0.28056	0.38159	0.12511E-04
0.32666	0.38273	0.11191E-04
0.37138	0.37988	0.99109E-05
0.41977	0.38219	0.84693E-05
0.42304	0.6187	0.80371E-05
0.42208	0.57362	0.81605E-05
0.42244	0.52708	0.82406E-05
0.42326	0.47961	0.82715E-05
0.42421	0.43344	0.82880E-05
0.22766	0.61746	0.13656E-04
0.23024	0.57171	0.13763E-04
0.22664	0.525	0.14004E-04
0.22524	0.47991	0.14139E-04
0.22612	0.43414	0.14155E-04
0.26315	0.63834	0.12506E-04
0.30356	0.63406	0.11376E-04
0.3503	0.63313	0.10086E-04
0.3893	0.63811	0.89778E-05
0.26061	0.41788	0.13105E-04
0.30356	0.43218	0.11832E-04
0.353	0.42328	0.10412E-04
0.39044	0.41495	0.93271E-05
0.38253	0.60024	0.92674E-05
0.37903	0.55227	0.94861E-05
0.38217	0.50242	0.94764E-05
0.3863	0.45459	0.94042E-05
0.272	0.5984	0.12416E-04
0.27874	0.54238	0.12375E-04
0.25837	0.49984	0.13085E-04
0.26209	0.45911	0.13038E-04
0.32887	0.58365	0.10828E-04
0.34332	0.47149	0.10641E-04
0.3359	0.52506	0.10761E-04
0.299	0.49143	0.11895E-04
0.17891	0.32851	0.15522E-04
0.47209	0.32851	0.68740E-05
0.17885	0.72142	0.14438E-04
0.47215	0.72142	0.63633E-05

kokyo52

楼主看一下
meshgrid
mesh
这两个命令的帮助

董雪娇

我这样编写的命令：
x=[];
y=[];
z=[];     %对应数据的xyz坐标分别设为三个向量
[X,Y]=meshgrid(min(x):max(x),min(y):max(y));
Z=griddata(x,y,z,X,Y);
mesh(X,Y,Z)
可是结果出现这样的错误：??? Error using ==> set
View [azimuth elevation] must be a two element vector.

Error in ==> mesh at 128
      set(cax,'View',y);
而且图形只有两轴坐标，没有图形。应该怎样修改程序代码呢？？

无水1324

help mesh 和meshgrid 是没有用的，你这样子做不可能做出来图，我原来试了很长时间

cyt07

可以的，有限元分析得出的数据我都是这么处理的。主要问题是你使用meshgrid时不合适。我用你的数据试了一下
>> x=indata(:,1);
>> y=indata(:,2);
>> z=indata(:,3);
>> xtemp=linspace(min(x),max(x),100);
>> ytemp=linspace(min(y),max(y),100);
>> [X,Y]=meshgrid(xtemp,ytemp);
>> Z=griddata(x,y,z,X,Y);
>> mesh(X,Y,Z)
图像的xy视图如下：

zhaoyan24

怎样生成XY、XZ或YZ视图？

messenger

用鼠标按着图旋转就行了
.

怎样生成XY、XZ或YZ视图？
zhaoyan24 发表于 2009-11-25 18:00

laihai718

太好了，正向请教这方面的问题

clystar

z点对应的是meshgrid产生的网格点,并非(x1,y1)所对应的z值,希望高手再指点下

messenger

已经讨论过了，参考此贴http://forum.simwe.com/thread-906236-1-3.html

z点对应的是meshgrid产生的网格点,并非(x1,y1)所对应的z值,希望高手再指点下
clystar 发表于 2009-11-26 19:31

xukyle

其实用Surfer软件来实现很方便的

xukyle

或利用delaunay 及trimesh函数
D=importdata('*.txt');
tri=delaunay(D(:,1),D(:,2));
trimesh(tri,D(:,1),D(:,2),D(:,3));
colormap gray


view(0,90)   % In-plane view

xukyle

plane view