new09编程2大学物理量子物理.pdf

流程控制 流程控制 循环结构 分支结构 if…else…end for…end （有单分支，双分支 （循环次数已定） 和多分支结构） switch…case… while…end otherwise…end (循环次数由条件定) （适用于多分支结构） Try结构 （尝试结构） for循环结构 Cantor set康托尔点集 长度为1的线段三等分后，去掉中间一段，再将剩下的两段再分别三等分，各去掉中间一段， 剩下更短的四段继续三等分后去中间一段，……，继续这种操作直至无穷。在极限的情况下， 线段长度趋于0，线段数目趋于无穷，相当于得到一个离散的点集，称为康托尔点集。 3 2.5 2 1.5 1 0.5 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 u=[0+3i,1+3i]; plot(u) %画直线 axis([0 1 0 3.1]), hold on for k=2:6 u=[u/3,NaN,u/3+2/3]+i*2/k; %上一段直线三等分， plot(u) % 去掉中间，画在纵坐标2/k处 end 矢量化编程 有时可以用矢量或矩阵来替代for循环，称之为矢量化编程。 在下例中，计算本来要对x,n 分别作for循环计算，即是二重循环， 实际是用矩阵运算来替代。 四条正弦曲线 1 x=[0:0.001:1]; % 变量x 是矢量 0.8 0.6 y=sin(pi*[1:4]*x); % 函数值y 是矩阵 0.4 plot(x,y) %用y 的每行对x作图