Matlab期末作业 李家兴20131336013.doc

Matlab基础及应用教程期末作业 学号 20131336013 姓名 李家兴 一、M文件及注释 用Matlab绘制管状几何体 利用MATLAB绘制管状几何体，所谓管状几何体就是中空的像水管一样的几何体。MATLAB函数TubeLike，可以根据管道中心线坐标及各点处管道半径绘制管状几何体。 程序M文件 function [X,Y,Z] = TubeLike(x,y,z,r) % 绘制三维管道型立体 % TubeLike(x,y,z) 绘制三维管道型几何体，输入参数x,y,z分别为管道中心线各点处 % 的坐标，x,y和z应为等长的向量，此时管道具有统一的半径1. % % TubeLike(x,y,z,r) 输入参数r用来指定管道半径。r可以是标量，也可以是与x,y,z等 % 长的向量。当r是标量时，管道具有统一半径；当r是向量时，管道 % 各截面处可以有不同的半径。 % % TubeLike(x) 输入参数x为3行或3列的矩阵，用来指定管道中心线各点处的坐标， % 此时管道具有统一的半径1. % % TubeLike(x,r) 输入参数x为3行或3列的矩阵，用来指定管道中心线各点处的坐标， % 输入参数r（标量或向量）用来指定管道半径。 % % [X,Y,Z] = TubeLike(...) 输出三维管道型几何体的网格数据X，Y和Z. % % Example: % t = linspace(0,2*pi,50); % x = sin(t); % y = cos(t); % z = cos(t/2); % r = sin(t); % TubeLike(x,y,z,r) if nargin=1 nargin=2 % 检查数据维数是否正确 [m,n] = size(x); p = min(m,n); % 维数 if p ~= 3 error(应输入三维样本数据,并且样本容量应大于3); end % 把样本观测值矩阵转置，使得行对应变量，列对应观测 if m = n x = x; end yd = x(2,:); zd = x(3,:); xd = x(1,:); if nargin == 1 r = ones(size(xd)); else if isvector(y) if numel(y) == 1 r = y*ones(size(xd)); elseif numel(y)1 numel(y) == length(xd) r = y(:); else error(半径应为标量或与x等长的向量); end else error(半径应为标量或与x等长的向量); end end elseif nargin=3 nargin=4 if isvector(x) isvector(y) isvector(z) numxyz = [numel(x),numel(y),numel(z)]; if any(numxyz-min(numxyz)) error(管道中心坐标x,y,z应为等长的向量); else xd = x(:); yd = y(:); zd = z(:); end else error(管道中心坐标x,y,z应为等长的向量); end if nargin == 3 r = ones(size(xd)); else if isvector(r) if numel(r) == 1 r = r*ones(size(xd)); elseif numel(r)1 nume