本帖最后由 valiant 于 2010-9-27 15:11 编辑
- x=[10 20 30 40 10 20 30 40 30 40 28 40];
- z=[10 10 10 10 20 21 19 20 30 30 40 40];
- t=[100 100 100 100 500 500 500 500 500 100 500 100 ];
- %假设存在0°90°180°270°四个方向,且数据一样。
- x1=x*cosd(0);y1=x*sind(0);
- x2=x*cosd(90);y2=x*sind(90);
- x3=x*cosd(180);y3=x*sind(180);
- x4=x*cosd(270);y4=x*sind(270);
- %end
- x=[x1;x2;x3;x4];
- y=[y1;y2;y3;y4];
- z=[z;z;z;z];
- t=[t;t;t;t];
- p1=-max(x1):1:max(x1); %三维坐标下x方向
- p2=-max(x1):1:max(x1); %三维坐标下y方向
- p3=0:1:max(z); %三维坐标下z方向
- [px,py] = meshgrid(p1,p2);
- F = TriScatteredInterp(x,y,t);
- pt=F(px,py);
- pt((pt>500))=NaN;pt((pt<10))=NaN;
- mesh(px,py,pt);shading interp; colorbar;
复制代码 还是不行,请各位高手明示啊。。 |