metlab的1.3例题及demo分析

1.3.2输入一个行矩阵






1.3.3分行输入一个行矩阵






1.3.4MATLAB提示出错






1.3.5用函数zeros生成全零阵






1.3.6用函数eye生成全零阵






1.3.7矩阵的加减运算






1.3.8两个矩阵的乘法运算






1.3.9矩阵的数乘运算






1.3.10向量的点积






1.3.11向量的叉乘






1.3.12向量的混合积






1.3.13左除和右除






1.3.14矩阵的乘方






1.3.15矩阵的转置






1.3.16对奇异矩阵求逆是MATLAB给出的警告信息






1.3.17用初等变换的方法求逆矩阵






1.3.18以有理格式输出结果






demo1 男孩的表面

这个例子的例子展示了如何创建男孩的表面。1901，希尔伯特给Werner Boy的任务，再次证明这是不可能的嵌入射影平面到空间无奇点。相反，男孩带着这个表面回来了。它有三个点，但没有奇点。更多信息，见男孩的表面的文章在维基百科，http://www.wikipedia.org。


n = 128;

u = linspace(0,pi,n);

v = linspace(0,pi,n);

u = repmat(u,n,1);

v = repmat(v',1,n);

x = cos(v).*sin(u);

y = sin(v).*sin(u);

z = cos(u);

f = 1/2*((2*x.^2-y.^2-z.^2) + 2*y.*z.*(y.^2-z.^2) +  ...

    z.*x.*(x.^2-z.^2) + x.*y.*(y.^2-x.^2));

g = sqrt(3)/2 * ((y.^2-z.^2) + z.*x.*(z.^2-x.^2) + x.*y.*(y.^2-x.^2));

h = (x+y+z).*((x+y+z).^3 + 4*(y-x).*(z-y).*(x-z));

clf

s = surf(f,g,h/10,u, ...

    'LineStyle','none', ...

    'FaceLighting','gouraud', ...

    'FaceColor','interp');

colormap jet;

axis off;

daspect([1 1 1]);

l1 = light;

l2 = light;

lightangle(l1,70,-40);

lightangle(l2,-30,80);

view(-40,32);

camzoom(1.5);




demo2四联系

这个例子示例说明如何生成四联解开环面旋转四偏心圆。感谢C. Henry Edwards，数学系，乔治亚大学。

ab = [0 2*pi];

rtr = [6 1 1];

pq = [10 50];

box = [-6.6 6.6 -6.6 6.6 -3 3];

vue = [200 70];

clf

tube('xylink1a',ab,rtr,pq,box,vue)

colormap(jet);

hold on

tube('xylink1b',ab,rtr,pq,box,vue)

tube('xylink1c',ab,rtr,pq,box,vue)

tube('xylink1d',ab,rtr,pq,box,vue)

hold off;