精细时程积分法matlab程序,MATLAB内燃机轴系振动计算程序

用于计算内燃机轴系振动的MTALAB程序,用基于高斯勒让德积分的精细时程积分法进行方程组求解,精度非常高。需要输入内燃机载荷。

djwy910=max(xyz1(10,0.0*nn+1:nn-0));

mixdjwy12=min(xyz1(2,0.0*nn+1:nn-0));

mixdjwy23=min(xyz1(3,0.0*nn+1:nn-0));

mixdjwy34=min(xyz1(4,0.0*nn+1:nn-0));

mixdjwy45=min(xyz1(5,0.0*nn+1:nn-0));

mixdjwy56=min(xyz1(6,0.0*nn+1:nn-0));

mixdjwy67=min(xyz1(7,0.0*nn+1:nn-0));

mixdjwy78=min(xyz1(8,0.0*nn+1:nn-0));

mixdjwy89=min(xyz1(9,0.0*nn+1:nn-0));

mixdjwy910=min(xyz1(10,0.0*nn+1:nn-0));

wn12=pi*(0.0001)^3/16;

wn23=pi*(0.21)^3/16;

wn34=pi*(0.21)^3/16;

wn45=pi*(0.21)^3/16;

wn56=pi*(0.21)^3/16;

wn67=pi*(0.21)^3/16;

wn78=pi*(0.21)^3/16;

wn89=pi*(0.3)^3/16;

wn910=pi*(0.22)^3/16;

xyz12=Txyz(2,0.0*nn+1:nn-