桁架的数据编程通常涉及以下步骤：
原始数据输入
输入节点数（Node Number, NE）、单元数（Element Number, NR）、杆件数（Number of Members, NP）等基本信息。
输入每个单元的类型（Unit Type）、截面面积（Cross-sectional Area, A）、弹性模量（Elastic Modulus, E）和材料容重（Weight Density, W）。
输入支座位移信息（Support Displacements）。
数据分类与存储
将具有相同弹性模量、截面面积和材料容重的单元归为一类，并存储在相应的数组中，如单元类型号数组（MT）、各类单元的弹性模量（EAW）、截面积（A）和材料容重（W）。
单元刚度矩阵计算
根据单元类型和几何参数，计算每个单元的刚度矩阵。这通常涉及复杂的数学公式和力学分析，可能包括单元应变能的计算和方向余弦的应用。
结构总刚度矩阵形成
将所有单元的刚度矩阵组装成结构的总刚度矩阵（K）。这通常涉及矩阵的三角化（Triangulation）和组装过程。
荷载施加与求解
输入外部荷载数据，形成右端项总外力向量（F）。
通过回代法或其他数值方法求解结构的平衡方程，得到节点的位移向量（u）和杆件的内力向量（q）。
结果输出
输出节点的位移、杆件的内力以及其他相关计算结果。
```matlab
% 初始化数据
E = 1.0e7； % 弹性模量
A = 1.5； % 截面面积
node_number = 6； % 节点数
element_number = 8； % 单元数
nc = [0,0；0,40；40,0；40,40；80,0；80,40]； % 支座位移信息
en = [1,3；1,4；2,4；3,4；3,5；4,5；4,6；5,6]； % 单元编号
ed = zeros（node_number, 2）； % 节点坐标
constraint = [1,1；1,2；2,1；2,2]； % 约束条件
% 设置自由度
dof = 0；
for i = 1:node_number
for j = 1:2
dof = dof + 1；
ed（i, j） = dof；
end
end
% 初始化刚度矩阵
ek = E * A * sparse（[1,0,-1,0；0,0,0,0；-1,0,1,0；0,0,0,0]）； % 单位矩阵
k = zeros（dof, dof）； % 结构总刚度矩阵
% 计算单元刚度矩阵并组装到总刚度矩阵中
for i = 1:element_number
% 这里需要根据具体的单元类型和几何参数计算单元刚度矩阵
% 假设每个单元的刚度矩阵为km
km = pin_jointed（i）； % 自定义函数，计算单元刚度矩阵
% 将单元刚度矩阵km组装到总刚度矩阵k中
k = k + km；
end
% 输出总刚度矩阵
disp（k）；
```
请注意，上述代码仅为示例，实际编程中需要根据具体的桁架结构和荷载条件进行详细计算和编程。