求解振动常微分方程组如何加入边界条件？
%机车质量矩阵Mv1=;
Mv2=;
%机车阻尼、刚度矩阵元素
Cv_11=;
Cv_12=[-Cs1 -Cs2;Cs1*a1*S -Cs2*a2*S];
Cv_21=[-Cs1 Cs1*a1*S;-Cs2 -Cs2*a2*S];
Cv_22=;
Kv_11=;
Kv_12=[-Ks1 -Ks2;Ks1*a1*S -Ks2*a2*S];
Kv_21=[-Ks1 Ks1*a1*S;-Ks2 -Ks2*a2*S];
Kv_22=;
…
M=;
C=;
K=;
g=9.8;
Ms=[(m1+a2*mv)*g;(m2+a1*mv)*g];%机车静载
F=;
%简支梁边界条件
%theta1=0,Omega1=0;Omega(j+1)=0,theta'(j+1)=0
function xdot=B_v_system(time,x,M,C,K,F,j)
NN=2*(j+1);%自由度数
xdot(1:(NN+4),1)=x((NN+5):2*(NN+4),1);
xdot((NN+5):2*(NN+4),1)=M\(F-K*x(1:(NN+4),1)-C*x((NN+5):2*(NN+4),1));
%车-桥系统
%X=;系统坐标
NN=2*j;
%仿真程序
t0=0;
h=0.01;
tfinal=5;
tspan=;
x0=zeros(2*NN+8,1);
=ode45('B_v_system',tspan,x0,options);
H为三次样条插值形函数，Omega1;theta1…为梁单元节点转角与挠度，因为M，C，K并非方阵，因此方程组中未知数多与方程数（少两个方程），是否应该加入简支梁边界条件，怎样加？？谢谢
页:
[1]