即使是做咸鱼，也要做最咸的那一条

在电影《功夫足球》里，周星驰有句台词，叫做：“做人如果没梦想，跟咸鱼有什么分别?”

国外喜欢把挤在人海里讨生活的人比做沙丁鱼，而国内特别倾向于咸鱼。

有一回，我见到了咸鱼的制作过程。先是宰杀，宰杀之后涂盐，然后挂在墙头上风干，一直到鱼眼变得泛白而空洞。

那时候，我忽然觉得，原来用咸鱼来形容失去梦想的人生，是一个很生动的比喻。

我们中的大多数人，或许就算保留梦想，为梦想奋斗了，可能这辈子都是这样普通，不会有什么大的成绩。

可梦想还是要有的，万一实现了呢？

哪怕梦想终不能实现，最后的结果依然是一条咸鱼，我也要做最不甘心的那条，饱含往事被风干了的那条，最有滋有味的那条。

是的，即使是做咸鱼，也要做最咸的那一条。

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