Matlab四位數(shù)據(jù)作圖(Slice)

問題：對如下數(shù)據(jù)，做出截面圖，并按相應(yīng)的顏色值設(shè)置數(shù)據(jù)；xyzs1222513321144202551926731sx,y,zslice函數(shù)繪出切面圖，也就x=1s=252120這三種用顏色表示的函數(shù)值。首先，看到slice命令的幫助文件中有個例子：[x,y,z]=meshgrid(-2:.2:2,-2:.25:2,-2:.16:2);v=x.*exp(-x.^2-y.^2-z.^2);xslice=[-1.2,.8,2];yslice=2;zslice=[-2,0];slice(x,y,z,v,xslice,yslice,zslice)colormaphsv這里的slice前三個參數(shù)，必須構(gòu)成網(wǎng)格數(shù)據(jù)，如圖：通過：[x,y,z]=meshgrid(-2:.2:2,-2:.25:2,-2:.16:2);figure;view(3);holdon;fori=1:size(x,1)forj=1:size(x,2)fork=1:size(x,3)plot3(x(i,j,k),y(i,j,k),z(i,j,k),'.');endendend這些點構(gòu)成了網(wǎng)格數(shù)據(jù)所以，我們需要制造網(wǎng)格數(shù)據(jù)，進行截面操作；closeall;clc;clear;A=[1222513321144202551926731];x=A(:,y=A(:,2)';z=A(:,3)';s=A(:,4)';%提取有效數(shù)據(jù)xb=min(x);xe=max(x);yb=min(y);ye=max(y);zb=min(z);ze=max(z);sb=min(s);se=max(s);x1=linspace(xb,xe,30);y1=linspace(yb,ye,30);z1=linspace(zb,ze,30);s0=linspace(sb,se,30);%對有效數(shù)據(jù)分割[x1y1z1]=meshgrid(x1,y1,z1);%生成網(wǎng)格數(shù)據(jù)點s1=zeros(size(x1,1),size(x1,2));fori=1:size(s1,1)s1(i,:)=s0;endst=zeros(size(x1,1),size(x1,2),size(x1,3));fori=1:size(st,3)st(:,:,i)=s1;ends1=st;%構(gòu)造第四維數(shù)據(jù)s1figure;holdon;view(3);slice(x1,y1,z1,s1,x,[],[]);這樣就做出來所謂的截面圖這個是通過slice(x1,y1,z1,s1,x,y,[]);來做的，但沒有顏色值設(shè)置選項；closeall;clc;clear;A=[1222513321144202551926731];x=A(:,y=A(:,2)';z=A(:,3)';s=A(:,4)';%對有效數(shù)據(jù)分割xb=min(x);xe=max(x);yb=min(y);ye=max(y);zb=min(z);ze=max(z);y1=linspace(yb,ye,30);z1=linspace(zb,ze,30);x1=ones(size(y1,1),size(y1,2))*x(1);[x1,y1]=meshgrid(x1,y1);z=zeros(length(z1));fori=1:size(z,1)z(i,:)=z1;endsurf(x1,y1,z);%這個時候產(chǎn)生了對應(yīng)于x1的網(wǎng)格截面圖，下面來定顏色[a,b]=size(z);rgb=ones(a,b);len=length(s)+1;slic=floor(linspace(1,a,len));fori=1:length(slic)-1rgb(slic(i):sli