第六讲 程序设计.docVIP

第六讲 程序设计(续) 如果想看一些内部函数,可以输入edit 函数名, 例如,edit magic. 书上内容解释 例1 应用switch，建立一个学生成绩管理数据库。 在editor中编写文件student.m S=struct(Name,Name,Score,Score,Rank,Rank); for i=1:5 switch S(i).Score case 100 S(i).Rank=满分; case a S(i).Rank=优秀; case b S(i).Rank=良好; case c S(i).Rank=中等; case d S(i).Rank=及格; otherwise S(i).Rank=不及格; end end disp([学生姓名 ,得分 ,等级]) for i=1:5 disp([S(i).Name,blanks(6),num2str(S(i).Score),blanks(6),S(i).Rank]); end 在命令窗口输入 for i=1:10 a{i}=89+i; b{i}=79+i; c{i}=69+i; d{i}=59+i; end Name={zhang,wang,li,zhao,liu}; Score={82,65,94,55,100}; Student 例2建立34矩阵A，使. a=ones(3,4); %该语句可要可不要. m=3;n=4; for i=1:m for j=1:n a(i,j)=1/(i+j-1); end end format rat a 例3 使用while语句,取任意整数，如是偶数，则用2除；否则乘3加1，重复此过程，直到整数变为1。 function c=collatz(n) c=n; while n1 if rem(n,2)==0 n=n/2; else n=3*n+1; end c=[c,n]; end 例4 生成非波那契（Leonardo Pisano Fibonacci was born in 1170）数。 如果每对兔子每个月生出一对兔子，且新生的兔子从第二个月开始就能生育，那么一年后由最初的那对兔子一共产生多少对兔子？ 用表示个月后兔子的对数，则有(第个月新生兔子对,成年兔子对，所以第个月有兔子对),这是一个递归函数，初始条件为。 （1） function f=fibonacci(n) % FIBONACCI(n) generates the first n Fibonacci numbers. f=zeros(n,1); f(1)=1; f(2)=2; for k=3:n f(k)=f(k-1)+f(k-2); end （2） function f=fibnum(n) %FIBNUM(n) generates the nth Fibonacci number. if n=1 f=1; else f=fibnum(n-1)+fibnum(n-2); %递归调用。 end 注：（1）与（2）的结果有什么区别？ 例5 function c=collatz(n) c=n; while n1 if rem(n,2)==0 n=n/2; else n=3*n+1; end c=[c,n]; end plot(c) hold on xlabel(X label) ylabel(Y label) title(collatz(n)) legend(图形) 例6 求使n！的n。 n=1; while prod(1:n)1.e100,n=n+1;end n 例7：已知以及对应的函数值，利用拉格朗日插值公式，计算插值函数 在编辑器中编写如下函数： function y=lagr1(X,Y,n,x) y=0.0; for k=1:n p=1.0; for j=1:n if j~=k p=p*(x-X(j))/(X(k)-X(j)); end end y=p*Y(k)+y; end 在命令窗口中输入如下命令： X=[0.0,0.1,0.195,0.3,0.401,0.5]; Y=[0.39894,0.39695,0.39142,0.38138,0.36812,0.35206]; n=6; y=lagr1(X,Y,n,0.15) y = 0.3945 如果在命令窗口中输入： syms x; lagr1(X,Y,n,x) 则输出的表达式