分块矩阵
1 iNode=[1,0;3,0;2,1]; 2 B=[iNode(2,:)-iNode(1,:);iNode(3,:)-iNode(1,:)]'; 3 H=[ B(1,1)^2,2* B(1,1)*B(2,1),B(2,1)^2;... 4 B(1,2)*B(1,1),B(1,2)*B(2,1)+B(1,1)*B(2,2),B(2,1)*B(2,2);... 5 B(1,2)^2,2* B(2,2)*B(1,2),B(2,2)^2]; 6 ve=[iNode(2,:)-iNode(1,:);iNode(3,:)-iNode(1,:);iNode(3,:)-iNode(2,:)]; 7 veLength=[norm(ve(1,:)),norm(ve(2,:)),norm(ve(3,:))]; 8 Q=?; 9 I3=eye(3); 10 D=blkdiag(I3,B',B',B',H,H,H,Q);
