结构力学大作业一、编号二、输出结果的过程（编写完程序后）1.输入A=N6的矩阵（单元数N不限） 2.通过函数k(A)得到结点位移列矩阵W4.通过函数li(W,A)得到各单元的局部坐标下的杆端力矩阵L。（该题的操作过程详见附录一，详细程序见附录二）三、程序设计过程根据上述的输出结果过程，可得，后期工作量较小，前期编程工作中，得到整体坐标下的单元刚度矩阵以及单元等效结点荷载都仅仅是码字工作，唯一有些难度的是定位向量的利用，即将整体坐标下的单元刚度矩阵以及单元等效结点荷载变换成该单元对应的整体刚度矩阵以及整体等效荷载。设计思路：定位向量的作用是定位，则可以试想根据定位向量得到行列变换矩阵，从而得到该单元对应的整体刚度矩阵。假设共有N个位移编号，则建立r=(N+6)(N+6)的单位矩阵。kk=(N+6)(N+6)的空矩阵，单元定位向量为rr(1,2,5,13,14,0)，整体坐标系下单元刚度矩阵为ke。其中kk左上角66置放该刚度矩阵ke；目标：通过矩阵行列变换，将kk的右下角NN变成该单元对应的整体刚度矩阵，即得到相应的行变换
