%《MATLAB及其在理工課程中的應(yīng)用指南》第七章例7－1－1程序
% 求平衡桿系的支撐反力
% 西安電子科技大學(xué)出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
G1=200; G2=100; L1= 2; L2 = sqrt (2) ;	%給原始參數(shù)賦值
theta1 =30*pi/180; theta2 = 45*pi/180;
% 設(shè)X=[ Nax, Nay, Nbx, Nby, Ncx, Ncy ],則按此次序，系數(shù)矩陣A,B可寫成下式
A=[1,0,0,0,1,0;0,1,0,0,0,1;0,0,0,0,-sin(theta1), cos(theta1);... 
0,0,1,0,-1,0;0,0,0,1,0,-1;0,0,0,0, sin(theta2), cos(theta2)]
B=[0;G1;G1/2*cos(theta1);0;G2;-G2/2*cos(theta2)]
X = A\B;					%用左除求解線性方程組
disp('Nax, Nay, Nbx, Nby, Ncx, Ncy')		% 顯示結(jié)果
disp(X')