在电影《功夫足球》里，周星驰有句台词，叫做：“做人如果没梦想，跟咸鱼有什么分别?”
国外喜欢把挤在人海里讨生活的人比做沙丁鱼，而国内特别倾向于咸鱼。
有一回，我见到了咸鱼的制作过程。先是宰杀，宰杀之后涂盐，然后挂在墙头上风干，一直到鱼眼变得泛白而空洞。
那时候，我忽然觉得，原来用咸鱼来形容失去梦想的人生，是一个很生动的比喻。
我们中的大多数人，或许就算保留梦想，为梦想奋斗了，可能这辈子都是这样普通，不会有什么大的成绩。
可梦想还是要有的，万一实现了呢？
哪怕梦想终不能实现，最后的结果依然是一条咸鱼，我也要做最不甘心的那条，饱含往事被风干了的那条，最有滋有味的那条。
是的，即使是做咸鱼，也要做最咸的那一条。
美国数学月刊
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 endend[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 endend[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 endend 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 endend[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 endend[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 endend 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 endend[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 endend[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 endend 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 endend[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 endend[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 endend 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 endend