例6-9.m

来自「这是一个MATLAB7.0基础与提高例题的所有源码」· M 代码 · 共 21 行

p=4000;
E=300e9;
L=2;                                     
A1=150e-6;                            
A2=200e-6;
A3=400e-6;
a1=pi/3;
a2=pi/2;
a3=3*pi/4;
L1=L/sin(a1);
L2=L/sin(a2);
L3=L/sin(a3);                    			  %计算杆长
k1=E*A1/L1;
k2=E*A2/L2;
k3=E*A3/L3;                      		  %计算刚度系数
D=[cos(a1),cos(a2),cos(a3),0,0;sin(a1),sin(a2),sin(a3),0,0;...
    1/k1,0,0,-cos(a1),-sin(a1);0,1/k2,0,-cos(a2),-sin(a2);...
    0,0,1/k3,-cos(a3),-sin(a3)];    		      %给系数矩阵赋值
B=[p;0;0;0;0];
format long,X=D\B						 %求解线性方程组，用长格式显示结果