v3d_import_block.m

来自「matlab 三维显示。所有代码都有。可以显示三维网格。」· M 代码 · 共 41 行

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [x,y,z,v] = v3d_import_block(xx,yy,zz,vv)
%
% in Blockkoordinaten umrechnen - Abweichung zu Blockkordinaten = 1.0% 
% --------------------------------------------------------------------
%
% Dimension von vv gleich (yy-1 * xx-1 * zz-1)
% Dimension von v gleich (y * x * z)
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Bl鯿ke berechnen - Abweichung zu Blockkordinaten = 1.0%
for i=1:length(xx)-1
    x(i*2-1)=xx(i)+((xx(i+1)-xx(i))/100);    
    x(i*2)=xx(i+1)-((xx(i+1)-xx(i))/100);
end
for i=1:length(yy)-1
    y(i*2-1)=yy(i)+((yy(i+1)-yy(i))/100);    
    y(i*2)=yy(i+1)-((yy(i+1)-yy(i))/100);
end
for i=1:length(zz)-1
    z(i*2-1)=zz(i)+((zz(i+1)-zz(i))/100);    
    z(i*2)=zz(i+1)-((zz(i+1)-zz(i))/100);
end

% v in allen neuen Punkten setzen
for i=1:length(xx)-1
    for j=1:length(yy)-1
        for k=1:length(zz)-1
            v(j*2-1,i*2-1,k*2-1)=vv(j,i,k);
            v(j*2-1,i*2-1,k*2)=vv(j,i,k);
            v(j*2-1,i*2,k*2-1)=vv(j,i,k);
            v(j*2-1,i*2,k*2)=vv(j,i,k);
            v(j*2,i*2-1,k*2-1)=vv(j,i,k);
            v(j*2,i*2-1,k*2)=vv(j,i,k);
            v(j*2,i*2,k*2-1)=vv(j,i,k);
            v(j*2,i*2,k*2)=vv(j,i,k);
        end
    end
end