求解一个振动方程的时候出现的问题，想请问下哪里出了问题呢

uiabluqp

function dy=test_fun(t,x)
Mtotal=6.8;
" Z- M" I* A. d9 q* k8 qMvp   =0.25;
Mf    =8.1;
Mm    =8.6;
# S( W  `' L0 @* b' RMb    =8.5;
Msp   =2.1;
Kfe1 =2.5*10^8;
7 A! u8 F& [# o# ^Kfe2 =3*10^8;
' G9 J& D) `* G; n& h, Z9 |Kb1  =5.11*10^8+1.3*10^7*sin(438*pi*t);
Kb2  =6.04*10^8+1.55*10^7*sin(1216*pi*t);
8 d. G) H1 O7 XK1=5.2*10^8;
* y1 G5 e9 W* }K2=6.3*10^8;
6 W% P. k, Q/ i8 Y  N4 VKvp=4.71*10^8;
Kvm=10^9;
8 D+ ^: F$ Z: K9 dF0=250;
* x( f, w4 O, ^  |3 p) In=1460;
9 p# J  t8 K- J# ~7 uCvp=1100;
Csb=370;
Cb1=1900+100*sin(438*pi*t);
/ l( G$ f4 z$ Q# ]3 n$ OCb2=2500+500*sin(1216*pi*t);
dy=zeros(12,1);
0 A- J4 `1 W7 h- vdy(1)=x(7);
6 X8 X) Q* j, C8 Z* o3 ldy(2)=x(8);
dy(3)=x(9);
dy(4)=x(10);
- n) V: r) N: ndy(5)=x(11);
dy(6)=x(12);
6 P+ {* s$ A( U1 ^" ?. [dy(7)=(-Cb1*(x(7)-x(9))-Cb2*(x(7)-x(4))-Cvp*(x(7)-x(8))-Csb*(x(7)-x(12))-Kfe1*x(1)-Kb1*(x(1)-x(3))-Kb2(x(1)-x(4))+F0*sin((2*pi*n/60)*t))/Mtotal;
) s) I6 u. k0 r+ i; I  h& w0 Q* c9 Sdy(8)=(Cvp*(x(7)-x(8))-Kvp*(x(2)-x(4)))/Mvp;
9 w) \* Q1 {7 c" Ady(9)=(Cb1*(x(7)-x(9))-Kfe2*x(3)+Kb1*(x(1)-x(3))-K1*(x(3)-x(4)))/Mf;
1 k' X0 B2 L  A7 T! B7 p3 K/ xdy(10)=(Cb2*(x(7)-x(10))+Kb2*(x(1)-x(4))+Kvp*(x(2)-x(4))+K1*(x(3)-x(4))-K2*(x(4)-x(5)))/Mm;
8 i6 J( d+ M, F3 \( fdy(11)=(K2*(x(4)-x(5))-Kvm*(x(5)-x(6)))/Mb;
dy(12)=(Csb*(x(7)-x(12))+Kvm*(x(11)-x(12)))/Msp;
6 W  a* D! X! p. Z% `end
) W1 T) Y" r* S3 P- v; H4 ]

NNNei256

看不懂啊

CCxiaom

帮顶

ubeautqq

dy(7) = ...   -Kb2(x(1)-x(4)) 这个位置应该是少了一个*号吧