25、MATLAB 矩阵运算与特殊矩阵全解析

MATLAB 矩阵运算与特殊矩阵全解析

1. 特殊矩阵及相关代码示例

1.1 节点受力分析代码

在节点 1 处施加水平方向 1000 lbf 的载荷，以下是修改后的 MATLAB 代码：

theta1=30 % angle in degrees
theta2=60 % angle in degrees
F1x=1000 % horizontal load
F1y=0 % vertical load
A= [-cosd(theta1),cosd(theta2),0,0,0,0
    -sind(theta1),-sind(theta2),0,0,0,0
    cosd(theta1),0,1,1,0,0
    sind(theta1),0,0,0,1,0
    0,-cosd(theta2),-1,0,0,0
    0,sind(theta2),0,0,0,1];
B= [F1x,-F1y,0,0,0,0]';
x=inv(A)*B;
x=A\B;

运行结果为：

x =
 -866.03  500.00  -250.00  1000.00  433.01  -433.01

数组中的第四个值对应节点 2 在 x 方向的反作用力，为 1000，符合预期。

1.2 特殊矩阵生成函数

1.2.1 ones 和 zeros 函数

• ones 函数用于创建全为 1 的矩阵， zeros <