黄金分割和 Fibonacci数 比较一下 goldfract(6)和 fibonacci(7 连分式:g=g+f, 黄金分割:f(k)=fk-1)+(k-2); n+1 n+1 H→>∞ MAtLaB n=40 f(2n)./f(1:n-1)-phi f- fibonacci(n) f2n)f(1n-1)
黄金分割和Fibonacci数 比较一下goldfract(6)和fibonacci(7) 连分式：g = g + f; 黄金分割：f(k) = f(k-1) +f(k-2); n n n f f = +1 = + → n n n f f 1 lim n = 40; f = fibonacci(n); f(2:n)./f(1:n-1) f(2:n)./f(1:n-1) - phi
Fibonacci的兔子 Fibonacci的兔子以黄金分割的速度增长。 p+,n=fm-I+f, 2 有两个解p和(1phi)fn=c1+c2(1- 由初始条件: 2φ-1 2φ-1 fn=2-1 (n+-(1-)”+) 是有理分式 积
Fibonacci的兔子 Fibonacci的兔子以黄金分割的速度增长。 n f n = c n = n−1 + n−2 ＋ f f f 1 2 = − n n f n c c (1 ) 有两个解phi和(1-phi) = 1 + 2 − 由初始条件： 2 1 (1 ) , 2 1 1 2 − − = − − = c c ( (1 ) ) 2 1 1 +1 +1 − − − = n n n f 是有理分式
Fibonacci的兔子 注意:没有半只兔子 MATLAB format long e n=(1:40)y; f=(phi^(n+1)-(1-phi).^n+1)/(2*phi-1) f=round(f
Fibonacci的兔子 注意：没有半只兔子☺ format long e n = (1:40)'; f = (phi.^(n+1)-(1-phi).^(n+1))/(2*phi-1) f = round(f)
魔方阵 MAtLaB 816 A-- magic(3) 357sm(A) 492 sum(a) sum(diag(a) sum(diag(fliud(a))) sum(1:9/3
魔方阵 A = magic(3) sum(A) sum(A’) sum(diag(A)) sum(diag(flipud(A))) sum(1:9)/3 8 1 6 3 5 7 4 9 2
魔方阵的八种组合 MAtLaB for k0: 3 逆时阵旋转 k*90度 rot90(A, k) rot90(A, k) e
魔方阵的八种组合 for k=0:3 rot90(A,k) rot90(A’,k) end 逆时阵旋转 k*90度