ea219.m

来自「用MATLAB程序解决高阶问题；通过20多个实例来说明线性代数在后续课和工程中的」· M 代码 · 共 12 行

% 《工程线性代数（MATLAB版）》第二章例题2.19程序ea219
% 陈怀琛，高淑萍，杨威合编，电子工业出版社，2007年6月
% 弹性梁的柔度矩阵
% 
clear
D=0.001*[5,4,3;4,4,3;3,3,6]	    % 输入柔度矩阵
f=[30;50;20]		% 给定力，求挠度
y=D*f				% 解题（1）
y1=[0;0;0.4] 		% 给定挠度，求力
K=inv(D)			% 求刚度矩阵
f1=K*y1			    % 解题（2）