mass.m

matlab从入门到精通第6章书后源码对初学者很有帮助

function M=mass(t,y,m1,m2,L,g)
% 创建矩阵，存储加权函数矩阵
M=zeros(6,6);
M(1,1)=1;
M(2,2)=m1+m2;
M(2,6)=-m2*L*sin(y(5));
M(3,3)=1;
M(4,4)=m1+m2;
M(4,6)=m2*L*cos(y(5));
M(5,5)=1;
M(6,2)=-L*sin(y(5));
M(6,4)=L*cos(y(5));
M(6,6)=L^2;
%M(:,1)=1;