ke ={"语文","数学","英语","心理","管理","检修","护理","哲学","历史","地理","化学"} x1 =["语文","数学"] x2 =["管理","英语"] x3 =["检修","哲学"] x=set(x1+x2+x3) 6 print(f"有人选:{x}") ke-x print(f"没人选:{s}")

这是一段 Python 代码，其输出结果为：

有人选: {'心理', '历史', '管理', '数学', '英语', '哲学', '护理', '检修', '地理', '化学', '语文'}
没人选: {'护理', '历史', '心理', '地理', '化学'}

首先定义了一个包含所有可选课程的集合 `ke`，然后定义了三个学生选择的课程集合 `x1`、`x2` 和 `x3`，将它们合并成一个集合 `x`，表示有人选的课程。然后通过 `ke-x` 操作得到了没人选的课程集合 `s`，最后将两个集合分别输出。

相关问题

刚度矩阵Ｋｅ ＝Ｋｅ Ｌ ＋Ｋｅ Ｎ怎么分别计算出来

在有限元分析中，梁单元的刚度矩阵可以分解为两个部分：局部坐标系下的弯曲刚度矩阵（KeL）和转移矩阵（N）。

局部坐标系下的弯曲刚度矩阵（KeL）是指在梁单元的局部坐标系下，考虑弯曲变形所引起的刚度变化。它可以通过对梁单元进行离散化，采用一定的数值方法进行计算得到。KeL是一个4x4的矩阵，其中的元素可以通过梁单元的材料参数和几何参数计算得到。

转移矩阵（N）是由梁单元的几何形状和节点间的位移关系所决定的。它可以将局部坐标系下的弯曲刚度矩阵转化为全局坐标系下的刚度矩阵。N是一个6x4的矩阵，其中的元素可以通过梁单元的节点坐标计算得到。

梁单元的刚度矩阵Ke可以表示为：

Ke = KeL + KeN

其中，KeN = N'KeLN，N'表示转置矩阵，KeLN表示由KeL通过转移矩阵N变换得到的全局坐标系下的弯曲刚度矩阵。

具体的，KeL的计算可以采用之前提到的 beam_stiffness_matrix 函数实现。而KeN的计算可以通过以下代码实现：

function [KeN] = beam_transfer_matrix(nodes)
% nodes: 梁单元的节点坐标

% 计算梁单元的长度、方向余弦、正交方向余弦
L = sqrt((nodes(2,1)-nodes(1,1))^2 + (nodes(2,2)-nodes(1,2))^2 + (nodes(2,3)-nodes(1,3))^2);
cx = (nodes(2,1)-nodes(1,1))/L;
cy = (nodes(2,2)-nodes(1,2))/L;
cz = (nodes(2,3)-nodes(1,3))/L;
sx = (nodes(2,2)-nodes(1,2))/(L*cx);
sy = (nodes(1,1)-nodes(2,1))/(L*cx);
sz = 0;
qx = cy*sz - cz*sy;
qy = cz*sx - cx*sz;
qz = cx*sy - cy*sx;

% 计算转移矩阵
KeN = zeros(6, 4);
KeN(1, 1) = cx;
KeN(2, 1) = cy;
KeN(3, 1) = cz;
KeN(4, 2) = qx;
KeN(5, 2) = qy;
KeN(6, 2) = qz;
KeN(4, 3) = sx;
KeN(5, 3) = sy;
KeN(6, 3) = sz;
KeN(4, 4) = qx;
KeN(5, 4) = qy;
KeN(6, 4) = qz;
end

该函数接受梁单元的节点坐标作为输入，返回一个6x4的转移矩阵。在计算梁单元的刚度矩阵时，我们可以先计算出KeL和KeN，然后将它们相加得到梁单元的刚度矩阵Ke。

Ｋｅ ＝ＫｅＬ ＋ＫｅN，刚度矩阵要计算这两个之和

在有限元分析中，单元刚度矩阵可以分为局部坐标系下的刚度矩阵 KeL 和全局坐标系下的刚度矩阵 KeN 两部分，即：

Ke = KeL + KeN

其中 KeL 和 KeN 都是 n x n 的矩阵，n 表示该单元中自由度的数量。KeL 是局部坐标系下的刚度矩阵，表示在局部坐标系下该单元的刚度矩阵；KeN 是全局坐标系下的刚度矩阵，表示将局部坐标系下的刚度矩阵通过旋转矩阵转换到全局坐标系下的刚度矩阵。

在程序实现中，我们先计算出局部坐标系下的刚度矩阵 KeL，然后通过旋转矩阵将其转换到全局坐标系下的刚度矩阵 KeN，最后将两个矩阵相加得到单元的刚度矩阵 Ke。

下面是一个MATLAB代码示例：

% 定义局部坐标系下的单元刚度矩阵 KeL
KeL = [k11, k12;
       k21, k22];

% 定义旋转矩阵 R
cos_theta = cos(theta);
sin_theta = sin(theta);
R = [cos_theta, sin_theta, 0, 0;
     -sin_theta, cos_theta, 0, 0;
     0, 0, cos_theta, sin_theta;
     0, 0, -sin_theta, cos_theta];

% 计算全局坐标系下的单元刚度矩阵 KeN
KeN = R' * KeL * R;

% 计算单元的刚度矩阵 Ke
Ke = KeL + KeN;

上述代码中，我们首先定义了局部坐标系下的单元刚度矩阵 KeL，然后通过旋转矩阵 R 将其转换到全局坐标系下的刚度矩阵 KeN。最后，我们将 KeL 和 KeN 相加得到单元的刚度矩阵 Ke。