matlab线性代数例题.docVIP

《数学实验》在线习题3 Matlab程序设计部分 分析向量组，，的线性相关性，找出它们的最大无关组，并将其余向理表示成最大无关组的线性组合。 解， a1=[1 2 3]; a2=[-1 -2 0]; a3=[0 0 1]; a4=[1 -2 -1]; a5=[2 4 6]; A=[a1,a2,a3,a4,a5] ; [R,S]=rref(A) r=length(S) R = 1.0000 0 0.3333 0 2.0000 0 1.0000 0.3333 0 0 0 0 0 1.0000 0 S = 1 2 4 r = 3 线性相关 a1,a2,a3,a4,a5 最大无关组是a1,a2,a4 其余向量的线性组合是a3=1/3a1+1/3a2 a5=2a1 计算行列式 的值。其中。 解， syms x1 x2 x3 x4 y1 y2 y3 y4 D=[x1^3 x1^2*y1 x1*y1^2 y1^3;x2^3 x2^2*y2 x2*y2^2 y2^3;x3^3 x3^2*y3 x3*y3^2 y3^3;x4^3 x4^2*y4 x4*y4^2 y4^3]; d=det(D) x1=2;x2=3;x3=5;x4=7;y1=4;y2=5;y3=6;y4=7; eval(d) d = ans = 153664 已知向量，，求向量a与b的夹角的度数。 解， a=[1 -1 0];b=[-1 0 -1]; x=a.*b; x1=sum(x,2); x2=norm(a); x3=norm(b); y=x1/(x2*x3) y1=acos(y) y = -0.5000 y1 = 2.0944 已知线性方程组 ，求系数矩阵的秩和方程组的通解。 clear a=[2 -1 3 2;9 -1 14 2;3 2 5 -4;4 5 7 -10]; b=[0 1 1 2]; [R,s]=rref([a,b]); [m,n]=size(a); x0=zeros(n,1); r=length(s); x0(s,:)=R(1:r,end); disp(非齐次线性方程组的特解为：) x0 disp(对应的线性方程组的基础解系为：) x=null(a,r) 非齐次线性方程组的特解为： x0 = 0.1429 0.2857 0 0 对应的齐次线性方程组的基础解系： x = -1.5714 0 -0.1429 2.0000 1.0000 0 0 1.0000 则方程组的通解为： 求齐次方程组的通解。 clear a=[-1 1 0 1;0 2 1 0;2 3 -1 -1]; b=[1 2 0]; [R,s]=rref([a,b]); [m,n]=size(a); x0=zeros(n,1); r=length(s); x0(s,:)=R(1:r,end); disp(非齐次线性方程组的特解为：) x0 disp(对应的线性方程组的基础解系为：) x=null(a,r) 非齐次线性方程组的特解为： x0 = -0.4286 0.5714 0.8571 0 对应的齐次线性方程组的基础解系： x = 0.8571 -0.1429 0.2857 1.0000 ，求正交矩阵P及对角形矩阵B，使。 clear a=[2 3 -2;3 6 11;-2 11 5]; [v,d]=eig(a) v = -0.3684 0.9280 0.0562 0.6512 0.2144 0.7280 -0.6635 -0.3047 0.6833 d = -6.9057 0 0 0 3.3500 0 0 0 16.5556 求下列向量的秩和最大无关组，并将其余向量用该最大无关组线性表出： a1=[1 2 1 3]; a2=[4 -1 -5 -6]; a3=[1 -3 -4 -7]; A=[a1,a2,a3] ; [R,S]=rref(A) r=length(S) R = 1.0000 0 -1.