在電影《功夫足球》里，周星馳有句台詞，叫做：“做人如果沒夢想，跟咸魚有什么分別?”
國外喜歡把擠在人海里討生活的人比做沙丁魚，而國內特別傾向於咸魚。
有一回，我見到了咸魚的制作過程。先是宰殺，宰殺之后塗鹽，然后掛在牆頭上風干，一直到魚眼變得泛白而空洞。
那時候，我忽然覺得，原來用咸魚來形容失去夢想的人生，是一個很生動的比喻。
我們中的大多數人，或許就算保留夢想，為夢想奮斗了，可能這輩子都是這樣普通，不會有什么大的成績。
可夢想還是要有的，萬一實現了呢？
哪怕夢想終不能實現，最后的結果依然是一條咸魚，我也要做最不甘心的那條，飽含往事被風干了的那條，最有滋有味的那條。
是的，即使是做咸魚，也要做最咸的那一條。
http://gen.lib.rus.ec/scimag/journaltable.php?journalid=1154
美國數學月刊
A=[-4 -1;2 0],B=[-2 1;3 -2],C=[-3 1;5 -2]; D=[1 0;0 1],E=[1 0;0 1],F=[1 0;0 1]; G=[2 0;0 2],H=[2 0;0 2],I=[2 0;0 2]; J=[1 0;0 1],K=[1 0;0 1],L=[1 0;0 1]; for k=1:4 M=-4*A*A'*D-2*A*D*A-2*A'*D*A'-2*A*G-2*G*A'-2*J+K+L-3*A+2*B+2*C, N=-4*B*B'*E-2*B*E*B-2*B'*E*B'-2*B*H-2*H*B'-2*K+J+L-3*B+2*A+2*C, O=-4*C*C'*F-2*C*F*C-2*C'*F*C'-2*C*I-2*I*C'-2*L+J+K-3*C+2*A+2*B; R=-G-2*D*A-2*A'*D, S=-H-2*E*B-2*B'*E, T=-I-2*F*C-2*C'*F; J=J-2*D+E+E, K=K+D-2*E+F, L=L+D+E-2*F; [V1,J1] = jordan((M+M')/2), [V2,J2] = jordan((N+N')/2), [V3,J3] = jordan((O+O')/2), D=V1*max(J1,[0 0;0 0])*V1', E=V2*max(J2,[0 0;0 0])*V2', F=V3*max(J2,[0 0;0 0])*V3', [W1,O1] = jordan((R+R')/2), [W2,O2] = jordan((S+S')/2), [W3,O3] = jordan((T+T')/2), [m1,n1]=size(O1); for i=1:m1 for j=1:n1 if O1(i,j)>=1 B1(i,j)=O1(i,j); else B1(i,j)=0; end end end [m2,n2]=size(O2); for i=1:m2 for j=1:n2 if O2(i,j)>=1 B2(i,j)=O2(i,j); else B2(i,j)=0; end end end [m3,n3]=size(O3); for i=1:m3 for j=1:n3 if O3(i,j)>=1 B3(i,j)=O3(i,j); else B3(i,j)=0; end end end G=W1*B1*W1', H=W2*B2*W2', I=W3*B3*W3', %k=k+1; end
A1=[-4 -1;2 0],A2=[-2 1;3 -2],A3=[-3 1;5 -2]; P1=[1 0;0 1],P2=[1 0;0 1],P3=[1 0;0 1]; Q1=[2 0;0 2],Q2=[2 0;0 2],Q3=[2 0;0 2]; M1=[1 0;0 1],M2=[1 0;0 1],M3=[1 0;0 1]; for k=1:5 M=-2*A1*A1'*P1-2*P1*A1*A1'-2*A1*P1*A1-2*A1'*P1*A1'-2*A1*Q1-2*Q1*A1'-2*M1+M2+M3-3*A1+2*A2+2*A3, N=-2*A2*A2'*P2-2*P2*A2*A2'-2*A2*P2*A2-2*A2'*P2*A2'-2*A2*Q2-2*Q2*A2'-2*M2+M1+M3-3*A2+2*A1+2*A3, O=-2*A3*A3'*P3-2*P3*A3*A3'-2*A3*P3*A3-2*A3'*P3*A3'-2*A3*Q3-2*Q3*A3'-2*M3+M1+M2-3*A3+2*A1+2*A2; R=-Q1-2*P1*A1-2*A1'*P1, S=-Q2-2*P2*A2-2*A2'*P2, T=-Q3-2*P3*A3-2*A3'*P3; M1=M1-2*P1+P2+P2, M2=M2+P1-2*P2+P3, M3=M3+P1+P2-2*P3; [V1,J1] = jordan((M+M')/2), [V2,J2] = jordan((N+N')/2), [V3,J3] = jordan((O+O')/2), P1=V1*max(J1,[0 0;0 0])*V1', P2=V2*max(J2,[0 0;0 0])*V2', P3=V3*max(J2,[0 0;0 0])*V3', [W1,O1] = jordan((R+R')/2), [W2,O2] = jordan((S+S')/2), [W3,O3] = jordan((T+T')/2), [m1,n1]=size(O1); for i=1:m1 for j=1:n1 if O1(i,j)>=1 B1(i,j)=O1(i,j); else B1(i,j)=0; end end end [m2,n2]=size(O2); for i=1:m2 for j=1:n2 if O2(i,j)>=1 B2(i,j)=O2(i,j); else B2(i,j)=0; end end end [m3,n3]=size(O3); for i=1:m3 for j=1:n3 if O3(i,j)>=1 B3(i,j)=O3(i,j); else B3(i,j)=0; end end end Q1=W1*B1*W1', Q2=W2*B2*W2', Q3=W3*B3*W3', %k=k+1; end
A1=[-4 -1;2 0],A2=[-2 1;3 -2],A3=[-3 1;5 -2]; P1=[1 0;0 1],P2=[1 0;0 1],P3=[1 0;0 1]; Q1=[2 0;0 2],Q2=[2 0;0 2],Q3=[2 0;0 2]; M1=[1 0;0 1],M2=[1 0;0 1],M3=[1 0;0 1]; for k=1:3 M=-2*A1*A1'*P1-2*P1*A1*A1'-2*A1*P1*A1-2*A1'*P1*A1'-2*A1*Q1-2*Q1*A1'-2*M1+M2+M3-3*A1+2*A2+2*A3, N=-2*A2*A2'*P2-2*P2*A2*A2'-2*A2*P2*A2-2*A2'*P2*A2'-2*A2*Q2-2*Q2*A2'-2*M2+M1+M3-3*A2+2*A1+2*A3, O=-2*A3*A3'*P3-2*P3*A3*A3'-2*A3*P3*A3-2*A3'*P3*A3'-2*A3*Q3-2*Q3*A3'-2*M3+M1+M2-3*A3+2*A1+2*A2; R=-Q1-2*P1*A1-2*A1'*P1, S=-Q2-2*P2*A2-2*A2'*P2, T=-Q3-2*P3*A3-2*A3'*P3; M1=M1+2*P1-P2-P3, M2=M2-P1+2*P2-P3, M3=M3-P1-P2+2*P3; [U1,S1,V1] = svd((M+M')/2), [U2,S2,V2] = svd((N+N')/2), [U3,S3,V3] = svd((O+O')/2), P1=U1*max(S1,[0 0;0 0])*V1, P2=U2*max(S2,[0 0;0 0])*V2, P3=U3*max(S2,[0 0;0 0])*V3, [U4,O1,V4] = svd((R+R')/2), [U5,O2,V5] = svd((S+S')/2), [U6,O3,V6] = svd((T+T')/2), [m1,n1]=size(O1); for i=1:m1 for j=1:n1 if O1(i,j)>=1 B1(i,j)=O1(i,j); else B1(i,j)=0; end end end [m2,n2]=size(O2); for i=1:m2 for j=1:n2 if O2(i,j)>=1 B2(i,j)=O2(i,j); else B2(i,j)=0; end end end [m3,n3]=size(O3); for i=1:m3 for j=1:n3 if O3(i,j)>=1 B3(i,j)=O3(i,j); else B3(i,j)=0; end end end Q1=U4*B1*V4, Q2=U5*B2*V5, Q3=U6*B3*V6, %k=k+1; end
A1=[-4 -1;2 0],A2=[-2 1;3 -2],A3=[-3 1;5 -2]; P1=[1 0;0 1],P2=[1 0;0 1],P3=[1 0;0 1]; Q1=[2 0;0 2],Q2=[2 0;0 2],Q3=[2 0;0 2]; M1=[1 0;0 1],M2=[1 0;0 1],M3=[1 0;0 1]; for k=1:4 M=-2*A1*A1'*P1-2*P1*A1*A1'-2*A1*P1*A1-2*A1'*P1*A1'-2*A1*Q1-2*Q1*A1'-2*M1+M2+M3-3*A1+2*A2+2*A3, N=-2*A2*A2'*P2-2*P2*A2*A2'-2*A2*P2*A2-2*A2'*P2*A2'-2*A2*Q2-2*Q2*A2'-2*M2+M1+M3-3*A2+2*A1+2*A3, O=-2*A3*A3'*P3-2*P3*A3*A3'-2*A3*P3*A3-2*A3'*P3*A3'-2*A3*Q3-2*Q3*A3'-2*M3+M1+M2-3*A3+2*A1+2*A2; R=-Q1-2*P1*A1-2*A1'*P1, S=-Q2-2*P2*A2-2*A2'*P2, T=-Q3-2*P3*A3-2*A3'*P3; M1=M1+2*P1-P2-P3, M2=M2-P1+2*P2-P3, M3=M3-P1-P2+2*P3; [U1,S1] = eig((M+M')/2), [U2,S2] = eig((N+N')/2), [U3,S3] = eig((O+O')/2), P1=U1*max(S1,[0 0;0 0])*U1', P2=U2*max(S2,[0 0;0 0])*U2', P3=U3*max(S2,[0 0;0 0])*U3', [U4,O1] = eig((R+R')/2), [U5,O2] = eig((S+S')/2), [U6,O3] = eig((T+T')/2), [m1,n1]=size(O1); for i=1:m1 for j=1:n1 if O1(i,j)>=1 B1(i,j)=O1(i,j); else B1(i,j)=0; end end end [m2,n2]=size(O2); for i=1:m2 for j=1:n2 if O2(i,j)>=1 B2(i,j)=O2(i,j); else B2(i,j)=0; end end end [m3,n3]=size(O3); for i=1:m3 for j=1:n3 if O3(i,j)>=1 B3(i,j)=O3(i,j); else B3(i,j)=0; end end end Q1=U4*B1*U4', Q2=U5*B2*U5', Q3=U6*B3*U6', %k=k+1; end
[u,v]=eig(A)
[V,D]=eig(A) %% 如果A是實對稱陣，那么eig A求出來的特征矩陣是單位正交陣。
format long g %位數更多
%%算1-正定投影
function B=change(A)
[m,n]=size(A);
for i=1:m
for j=1:n
if A(i,j)>128
B(i,j)=1;
else B(i,j)=0;
end
end
end