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

uiabluqp

function dy=test_fun(t,x)
1 a, s2 [; q+ Y9 P' c6 ]$ ^Mtotal=6.8;
' @9 [. T1 I4 Q4 a/ ?. y" V' bMvp   =0.25;
4 j$ A- v& u  O8 VMf    =8.1;
5 v) Z6 p: \2 ?Mm    =8.6;
3 D: [: W" L' O6 ]: D6 {Mb    =8.5;
8 E$ T" [  j: I9 y* s& IMsp   =2.1;
Kfe1 =2.5*10^8;
Kfe2 =3*10^8;
Kb1  =5.11*10^8+1.3*10^7*sin(438*pi*t);
3 L9 M* o& j' t+ UKb2  =6.04*10^8+1.55*10^7*sin(1216*pi*t);
2 I! K/ Z% L# b. M5 A0 k" E2 BK1=5.2*10^8;
K2=6.3*10^8;
5 O( f7 x( U, a$ q7 s4 y; V: T4 AKvp=4.71*10^8;
Kvm=10^9;
( t0 m# P7 s* g& F3 H, [F0=250;
n=1460;
# I# N, Q9 e! x/ S3 J" f7 z3 E( fCvp=1100;
Csb=370;
) d8 i* @* T, M9 E9 M8 \5 H- Y! xCb1=1900+100*sin(438*pi*t);
; ]( `6 |' D2 W9 P6 ?$ _. l/ t8 XCb2=2500+500*sin(1216*pi*t);
dy=zeros(12,1);
- R" v+ l4 X: }4 j8 C9 Ndy(1)=x(7);
" T! [0 w5 g: V4 W4 y) `8 i% }0 f$ ?dy(2)=x(8);
; Q/ ?/ @* k" N7 T" d8 wdy(3)=x(9);
dy(4)=x(10);
( l. M+ |3 \- q' Ydy(5)=x(11);
dy(6)=x(12);
' W& J! `. L+ a5 ^  [4 W: _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;
" E! K* k7 l* X) N$ o( E5 l5 V$ Idy(8)=(Cvp*(x(7)-x(8))-Kvp*(x(2)-x(4)))/Mvp;
7 x3 ~2 d* i" t9 N' f+ l/ vdy(9)=(Cb1*(x(7)-x(9))-Kfe2*x(3)+Kb1*(x(1)-x(3))-K1*(x(3)-x(4)))/Mf;
dy(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;
dy(11)=(K2*(x(4)-x(5))-Kvm*(x(5)-x(6)))/Mb;
# I* ^* `* s# f% C: Z, gdy(12)=(Csb*(x(7)-x(12))+Kvm*(x(11)-x(12)))/Msp;
end

NNNei256

看不懂啊

CCxiaom

帮顶

ubeautqq

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