MATLAB程序设计概要.ppt

n=1000; a=1; for m=1:n; a=a*m; if a=1e100 break end end a, m 例：求解n!为100位数的最小n值。 一个数列，它的后一元素是前两个元素的积，即x(n+2)=x(n+1)*x(n)，已知x(1)=1，x(2)=2。编写求此数列的程序，n的最大值由键盘输入，程序最后输出该数列。 n=input(n=); x=1:n; for m=3:n x(m)=x(m-1)*x(m-2); end x 例：有一级数， ，编写程序，计算当此级数收敛到10－5时，级数的项数N，及SN的值。 n=1;s=0; while n^-2=1e-5 s=s+n^-2; n=n+1; end n-1 s-n^-2 n=1;s=0; for n=1:1000 if n^-21e-5 break end s=s+n^-2; end n,s 例：编写一程序，求魔方矩阵magic(20)对角线上的元素，并以一行向量表示。 a=magic(20); s=[]; for m=1:20 for n=1:20 if m==n s=[s a(m,n)]; end end end s 有如下递推关系，编程分别计算精确到前50项和前100项的值。 y=0; p=1; for n=1:50 p=p/n; y=y+p*x^n; end y y=0; p=1; for n=1:100 p=p*x/n; y=y+p; end disp(y); 存在如下递推关系： 编程计算精确到前100项的值,x由键盘输入。 y=0; p=1; for n=1:100 p=p*(-x)/n; y=y+p; end y for循环语句注意事项 (1) for循环不能用循环内重新给循环变量赋值来终止。 (2) for循环的循环变量=[表达式1:表达式2:表达式3]，其实为一行向量，它还可以是数组，其更一般的形式为 for 循环控制变量=数组表达式 语句 end 用 for 循环求行向量[-2,5,3,6,-2]各元素之和。 a=[-2,5,3,6,-2]; s=0; k=0; for n=a n k=k+1; s=s+n; end k,s sam6_5 n = -2 n = 5 n = 3 n = 6 n = -2 k = 5 s = 10 a=[-2,5,3,6,-2]; s=0; for n=a s=s+n; end s a=[-2,5,3,6,-2]; s=[]; r=0; for n=a r=r+n; s=[s r]; end s cumsum(a) sam6_6 data = 3 9 45 6 7 16 -1 5 n = 3 7 x = -4 …… k = 4 data=[3 9 45 6; 7 16 -1 5] k=0; for n=data n k=k+1; x=n(1)-n(2) end k data=[3 9 45 6; 7 16 -1 5] d=[]; for n=data d=[d,n(1)-n(2)]; end d (3) for 循环可嵌套使用。 以 for 循环求1！+2！+…+10!的值。 s=0; for m=1:10 p=1; for n=1:m p=p*n; end s=s+p; end s (4) 当有一个等效的数组方法来解给定的问题时，应避免用for 循环。 for n=1:10 x(n)=sin(n*pi/10); end x n=1:10;x=sin(n*pi/10) (5) 为了得到更快的速度，在for 循环(while 循环)被执行之前，应预先分配数组。否则在for 循环内每执行一次命令，变量x的大小增加1，迫使MATLAB 每进行一次循环都要花费时间对x分配更多的内存。为了省去这个步骤，可以在程序的首行加入： x=zeros(1,10) 2. while 循环 for循环的循环次数是固定的，而whi