MATLAB编程求解水仙花数.docxVIP

【原创】MATLAB编程--求解水仙花数 【例】何为水仙花数？ 水仙花数是指一个n位数（n 3,）它的每个位上的数字的n次幕之和等于它本 身。（例如：1A3 + 5A3 +3A3 = 153 ） 三位的水仙花数共有 4 个: 153, 370, 371, 407; 四位的水仙花数共有 3 个: 1634 , 8208 , 9474 ; 五位的水仙花数共有 3 个: 54748 , 92727 , 93084 ; 六位的水仙花数只有 1 个: 548834 ; 七位的水仙花数共有 4 个: 1741725 , 4210818 , 9800817 , 9926315 八位的水仙花数共有 3 个:,,程序如下： for m=100:999 m1=fix(m/100); m2=rem(fix(m/10),10); m3=rem(m,10); if m==m1A3+m2A3+m3A3 disp(m) end end 153 370 371 407 while循环结构 【例】Fibo nacci数组的元素满足Fib on acci规则：ak 2 ak ak 1 , (k 1,2,)； 且a1 a2 1。现要求该数组中第一个大于10000的元素。 a(1)=1;a (2)=1;i=2; while a(i)=10000 a(i+1)=a(i-1)+a(i); i=i+1; end; i,a(i), 21 ans = 10946 【例】用for循环指令来寻求Fib on acc数组中第一个大于10000的元素 n=100;a=ones(1,n); for i=3:n a(i)=a(i-1)+a(i-2); if a(i)=10000 a(i), break; end; end,i ans = 10946 i = 21 switch-case 结构 【例】学生的成 绩管理，用来演示switch结构的应用 clear; for i=1:10;a{i}=89+i;b{i}=79+i;c{i}=69+i;d{i}=59+i;end;c=[d,c]; Name={ Jack,Marry,Peter, Rose, Tom}; Mark={72,83,56,94,100};Rank=cell(1,5); S=struct(Name,Name,Marks,Mark,Rank,Rank); for i=1:5 满分;优秀; 满分; 优秀; 良好; 及格;不及格 ; 得分 等级]);disp() S(i).Rank= case a S(i).Rank= case b S(i).Rank= case c S(i).Rank= otherwise S(i).Rank= end end disp([ 学生姓名 for i=1:5; disp([S(i).Name,blanks(6),num2str(S(i).Marks),blanks(6),S(i) .Rank]); end; 学生姓名 得分 等级 Jack 72 及格 Marry 83 良好 不及格Peter 56 不及格 Rose 94 优秀 Tom 100 满分 【例】 matlab 魔方矩 阵的程序 设计 function M = magic(n) %MAGIC Magic square. % MAGIC(N) is an N-by-N matrix constructed from the integers % 1 through NA2 with equal row, colum n, and diago nal sums. % Produces valid magic squares for all N 0 except N = 2. % Copyright 1984-2002 The MathWorks, Inc. % $Revision: 5.15 $ $Date: 2002/04/15 03:44:23 $ % Historically, MATLABs magic was a built-in function. % This M-file uses a new algorithm to generate the same matrices. n = floor(real(double(n(1)))); % Odd order. if mod(n,2) == 1 [J,I] = meshgrid(1:n); A = mod(I+J-(n+3)/2,n); B = mod(I+2*J-2,n); M = n*A + B + 1; % Doubly even order. elseif mod(n,4) == 0 [J,I] = me