谢谢!

这时一线性系统,数值计算应该没问题,按下列程序计算的特征值最大正值是0.7023,但用basic计算的结果其增长率与正的特征值为4.79时一般大,但这里的结果为-4.79.很迷惑.谢谢!

m=2;k0=0.65;k1=0.6;k3=k1-k0;k2=k1+k0;k12=k1^2+m^2;k22=k2^2+m^2;k32=k3^2+m^2;

v0=0.5;arr=0.01;arf=0.04;vc=v0;a01=(4-k12)^2+16*k1^2;a02=(4-k22)^2+16*k2^2;a03=(4-k32)^2+16*k3^2;

a1=(4-k12-4*k1*i)/a01;a2=(4-k22-4*k2*i)/a02;a3=(4-k32-4*k3*i)/a03

qq11=0.5*v0*m*(1+a1*(1+k0^2)); qq12=0.5*v0*m*(1+a2*(1+k0^2)); qq13=0.5*v0*m*(1+a3*(1+k0^2));

qq21=-0.5*v0*m*(1+a1*(1+k0^2));qq22=-0.5*v0*m*(1+a2*(1+k0^2));qq23=-0.5*v0*m*(1+a3*(1+k0^2));

qqv1=i*m*vc*(1+a1);qqv2=i*m*vc*(1+a2);qqv3=i*m*vc*(1+a3);

qd11=-0.5*v0*((1+k0^2)*(k1-2*i)*a1-k0+i);qd12=-0.5*v0*((1+k0^2)*(k2-2*i)*a2-k0+i);qd13=-0.5*v0*((1+k0^2)*(k3-2*i)*a3-k0+i);

qd21=0.5*v0*((1+k0^2)*(k1-2*i)*a1+k0+i); qd22=0.5*v0*((1+k0^2)*(k2-2*i)*a2+k0+i); qd23=0.5*v0*((1+k0^2)*(k3-2*i)*a3+k0+i);

qdv1=-i*v0*(k1-2*i)*a1+vc;qdv2=-i*v0*(k2-2*i)*a2+vc;qdv3=-i*v0*(k3-2*i)*a3+vc;

dq11=-v0*(k0-(k0+i)*(2+4*i*k1-k1^2)*a1); dq12=-v0*(k0-(k0+i)*(2+4*i*k2-k2^2)*a2); dq13=-v0*(k0-(k0+i)*(2+4*i*k3-k3^2)*a3);

dq21=-v0*(k0-(k0-i)*(2+4*i*k1-k1^2)*a1); dq22=-v0*(k0-(k0-i)*(2+4*i*k2-k2^2)*a2); dq23=-v0*(k0-(k0-i)*(2+4*i*k3-k3^2)*a3);

dd11=v0*m*(1+2*(k0+i)*(i-k1)*a1)*0.5; dd12=v0*m*(1+2*(k0+i)*(i-k2)*a2)*0.5; dd13=v0*m*(1+2*(k0+i)*(i-k3)*a3)*0.5;

dd21=-v0*m*(1-2*(k0-i)*(i-k1)*a1)*0.5; dd22=-v0*m*(1-2*(k0-i)*(i-k2)*a2)*0.5; dd23=-v0*m*(1-2*(k0-i)*(i-k3)*a3)*0.5;

dqv1=-2*vc*(2+4*i*k1-k1^2)*a1;dqv2=-2*vc*(2+4*i*k2-k2^2)*a2;dqv3=-2*vc*(2+4*i*k3-k3^2)*a3;

ddv1=i*m*vc-2*m*vc*(i-k1)*a1;ddv2=i*m*vc-2*m*vc*(i-k2)*a2;ddv3=i*m*vc-2*m*vc*(i-k3)*a3;

zz0=1

lq1=zz0*i*m*a1;lq2=zz0*i*m*a2;lq3=zz0*i*m*a3;

ld1=-zz0*(2+i*k1)*a1;ld2=-zz0*(2+i*k2)*a2;ld3=-zz0*(2+i*k3)*a3;

g1=-sqrt(k12)*(arr*sqrt(k12)-arf);g2=-sqrt(k22)*(arr*sqrt(k22)-arf);g3=-sqrt(k32)*(arr*sqrt(k32)-arf);

b11=qqv1;b12=qq22;b13=qq13;b14=qdv1;b15=qd22;b16=qd13;b17=0;b18=0;b19=0;

b21=qq11;b22=qqv2;b23=0;b24=qd11;b25=qdv2;b26=0;b27=0;b28=0;b29=0;

b31=qq21;b32=0;b33=qqv3;b34=qd21;b35=0;b36=qdv3;b37=0;b38=0;b39=0;

b41=dqv1;b42=dq22;b43=dq13;b44=ddv1;b45=dd22;b46=dd13;b47=g1;b48=0;b49=0;

b51=dq11;b52=dqv2;b53=0;b54=dd11;b55=ddv2;b56=0;b57=0;b58=g2;b59=0;

b61=dq21;b62=0;b63=dqv3;b64=dd21;b65=0;b66=ddv3;b67=0;b68=0;b69=g3;

b71=lq1;b72=0;b73=0;b74=1+ld1;b75=0;b76=0;b77=i*m*vc;b78=-0.5*v0*m*i;b79=0.5*v0*m*i;

b81=0;b82=lq2;b83=0;b84=0;b85=1+ld2;b86=0;b87=0.5*v0*m*i;b88=i*m*vc;b89=0;

b91=0;b92=0;b93=lq3;b94=0;b95=0;b96=1+ld3;b97=-0.5*v0*m*i;b98=0;b99=i*m*vc;

A1=[real(b11),-imag(b11),real(b12),-imag(b12),real(b13),-imag(b13),real(b14),-imag(b14),real(b15),-imag(b15),real(b16),-imag(b16),real(b17),-imag(b17),real(b18),-imag(b18),real(b19),-imag(b19)];

A2=[imag(b11),real(b11),imag(b12),real(b12),imag(b13),real(b13),imag(b14),real(b14),imag(b15),real(b15),imag(b16),real(b16),imag(b17),real(b17),imag(b18),real(b18),imag(b19),real(b19)];

A3=[real(b21),-imag(b21),real(b22),-imag(b22),real(b23),-imag(b23),real(b24),-imag(b24),real(b25),-imag(b25),real(b26),-imag(b26),real(b27),-imag(b27),real(b28),-imag(b28),real(b29),-imag(b29)];

A4=[imag(b21),real(b21),imag(b22),real(b22),imag(b23),real(b23),imag(b24),real(b24),imag(b25),real(b25),imag(b26),real(b26),imag(b27),real(b27),imag(b28),real(b28),imag(b29),real(b29)];

A5=[real(b31),-imag(b31),real(b32),-imag(b32),real(b33),-imag(b33),real(b34),-imag(b34),real(b35),-imag(b35),real(b36),-imag(b36),real(b37),-imag(b37),real(b38),-imag(b38),real(b39),-imag(b39)];

A6=[imag(b31),real(b31),imag(b32),real(b32),imag(b33),real(b33),imag(b34),real(b34),imag(b35),real(b35),imag(b36),real(b36),imag(b37),real(b37),imag(b38),real(b38),imag(b39),real(b39)];

A7=[real(b41),-imag(b41),real(b42),-imag(b42),real(b43),-imag(b43),real(b44),-imag(b44),real(b45),-imag(b45),real(b46),-imag(b46),real(b47),-imag(b47),real(b48),-imag(b48),real(b49),-imag(b49)];

A8=[imag(b41),real(b41),imag(b42),real(b42),imag(b43),real(b43),imag(b44),real(b44),imag(b45),real(b45),imag(b46),real(b46),imag(b47),real(b47),imag(b48),real(b48),imag(b49),real(b49)];

A9=[real(b51),-imag(b51),real(b52),-imag(b52),real(b53),-imag(b53),real(b54),-imag(b54),real(b55),-imag(b55),real(b56),-imag(b56),real(b57),-imag(b57),real(b58),-imag(b58),real(b59),-imag(b59)];

A10=[imag(b51),real(b51),imag(b52),real(b52),imag(b53),real(b53),imag(b54),real(b54),imag(b55),real(b55),imag(b56),real(b56),imag(b57),real(b57),imag(b58),real(b58),imag(b59),real(b59)];

A11=[real(b61),-imag(b61),real(b62),-imag(b62),real(b63),-imag(b63),real(b64),-imag(b64),real(b65),-imag(b65),real(b66),-imag(b66),real(b67),-imag(b67),real(b68),-imag(b68),real(b69),-imag(b69)];

A12=[imag(b61),real(b61),imag(b62),real(b62),imag(b63),real(b63),imag(b64),real(b64),imag(b65),real(b65),imag(b66),real(b66),imag(b67),real(b67),imag(b68),real(b68),imag(b69),real(b69)];

A13=[real(b71),-imag(b71),real(b72),-imag(b72),real(b73),-imag(b73),real(b74),-imag(b74),real(b75),-imag(b75),real(b76),-imag(b76),real(b77),-imag(b77),real(b78),-imag(b78),real(b79),-imag(b79)];

A14=[imag(b71),real(b71),imag(b72),real(b72),imag(b73),real(b73),imag(b74),real(b74),imag(b75),real(b75),imag(b76),real(b76),imag(b77),real(b77),imag(b78),real(b78),imag(b79),real(b79)];

A15=[real(b81),-imag(b81),real(b82),-imag(b82),real(b83),-imag(b83),real(b84),-imag(b84),real(b85),-imag(b85),real(b86),-imag(b86),real(b87),-imag(b87),real(b88),-imag(b88),real(b89),-imag(b89)];

A16=[imag(b81),real(b81),imag(b82),real(b82),imag(b83),real(b83),imag(b84),real(b84),imag(b85),real(b85),imag(b86),real(b86),imag(b87),real(b87),imag(b88),real(b88),imag(b89),real(b89)];

A17=[real(b91),-imag(b91),real(b92),-imag(b92),real(b93),-imag(b93),real(b94),-imag(b94),real(b95),-imag(b95),real(b96),-imag(b96),real(b97),-imag(b97),real(b98),-imag(b98),real(b99),-imag(b99)];

A18=[imag(b91),real(b91),imag(b92),real(b92),imag(b93),real(b93),imag(b94),real(b94),imag(b95),real(b95),imag(b96),real(b96),imag(b97),real(b97),imag(b98),real(b98),imag(b99),real(b99)];

A=-[A1;A2;A3;A4;A5;A6;A7;A8;A9;A10;A11;A12;A13;A14;A15;A16;A17;A18]

d=eig(A,'nobalance')