资源描述:
《MATLAB求线性方程组.ppt》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、第三章线性方程组§3.1求线性方程组的唯一解或特解一.用克拉默法则例3.1.1.求方程组的解.第三章线性方程组§3.1求线性方程组的唯一解或特解>>a_1=[5;1;0;0;0];a_2=[6;5;1;0;0];>>a_3=[0;6;5;1;0];a_4=[0;0;6;5;1];>>a_5=[0;0;0;6;5];b=[1;0;0;0;1];>>a_1=[5;1;0;0;0];a_2=[6;5;1;0;0];>>a_3=[0;6;5;1;0];a_4=[0;0;6;5;1];>>a_5=[0;0;0;6;5];b=[1;0;0;0;1
2、];>>D=det([a_1,a_2,a_3,a_4,a_5]);>>D_1=det([b,a_2,a_3,a_4,a_5]);>>D_2=det([a_1,b,a_3,a_4,a_5]);>>D_3=det([a_1,a_2,b,a_4,a_5]);>>D_4=det([a_1,a_2,a_3,b,a_5]);>>D_5=det([a_1,a_2,a_3,a_4,b]);>>a_1=[5;1;0;0;0];a_2=[6;5;1;0;0];>>a_3=[0;6;5;1;0];a_4=[0;0;6;5;1];>>a_5=[0;0;0;6;5
3、];b=[1;0;0;0;1];>>D=det([a_1,a_2,a_3,a_4,a_5]);>>D_1=det([b,a_2,a_3,a_4,a_5]);>>D_2=det([a_1,b,a_3,a_4,a_5]);>>D_3=det([a_1,a_2,b,a_4,a_5]);>>D_4=det([a_1,a_2,a_3,b,a_5]);>>D_5=det([a_1,a_2,a_3,a_4,b]);>>x_1=D_1/D;x_2=D_2/D;x_3=D_3/D;x_4=D_4/D;>>x_5=D_5/D;>>a_1=[5;1;0;0;0
4、];a_2=[6;5;1;0;0];>>a_3=[0;6;5;1;0];a_4=[0;0;6;5;1];>>a_5=[0;0;0;6;5];b=[1;0;0;0;1];>>D=det([a_1,a_2,a_3,a_4,a_5]);>>D_1=det([b,a_2,a_3,a_4,a_5]);>>D_2=det([a_1,b,a_3,a_4,a_5]);>>D_3=det([a_1,a_2,b,a_4,a_5]);>>D_4=det([a_1,a_2,a_3,b,a_5]);>>D_5=det([a_1,a_2,a_3,a_4,b]);>>
5、x_1=D_1/D;x_2=D_2/D;x_3=D_3/D;x_4=D_4/D;>>x_5=D_5/D;>>formatrat,X=[x_1,x_2,x_3,x_4,x_5]>>a_1=[5;1;0;0;0];a_2=[6;5;1;0;0];>>a_3=[0;6;5;1;0];a_4=[0;0;6;5;1];>>a_5=[0;0;0;6;5];b=[1;0;0;0;1];>>D=det([a_1,a_2,a_3,a_4,a_5]);>>D_1=det([b,a_2,a_3,a_4,a_5]);>>D_2=det([a_1,b,a_3,a_
6、4,a_5]);>>D_3=det([a_1,a_2,b,a_4,a_5]);>>D_4=det([a_1,a_2,a_3,b,a_5]);>>D_5=det([a_1,a_2,a_3,a_4,b]);>>x_1=D_1/D;x_2=D_2/D;x_3=D_3/D;x_4=D_4/D;>>x_5=D_5/D;>>formatrat,X=[x_1;x_2;x_3;x_4;x_5]X=1507/665-229/13337/35-79/133212/665第三章线性方程组§3.1求线性方程组的唯一解或特解二.用矩阵除法>>%把该方程组记为AX
7、=b,则X=Ab>>A=[5,6,0,0,0;1,5,6,0,0;0,1,5,6,0;0,0,1,5,6;0,0,0,1,5];>>b=[1;0;0;0;1];formatrat,X=Ab>>%把该方程组记为AX=b,则X=Ab>>A=[5,6,0,0,0;1,5,6,0,0;0,1,5,6,0;0,0,1,5,6;0,0,0,1,5];>>b=[1;0;0;0;1];formatrat,X=AbX=1507/665-229/13337/35-79/133212/665第三章线性方程组§3.1求线性方程组的唯一解或特解三.用矩阵
8、的初等变换>>A=[5,6,0,0,0;1,5,6,0,0;0,1,5,6,0;0,0,1,5,6;0,0,0,1,5];>>b=[1;0;0;0;1];>>B=[A,b];%增广矩阵>>A=