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;