第1个回答 2010-12-31
clear;clc;
A=[3,5,7;6,1,2;-8,9,13];
B=[1,-4,6;2,9,4;7,3,-2];
C=[A,B]
D=[A;B]
C(1,:)=[];
C;
C1=C(:,1);
C2=C(:,4);
C3=C(:,5);
P=[C1,C2,C3]
syms c;
E=[
2 7 3 1
3 5 2 2
9 4 1 7];
b=[6 4 2]';
b=b(1:3)+E(1:3,4).*c;
E=E(1:3,1:3);
E=sym(E);b=sym(b);
jie=E\b %通解
第2个回答 2010-12-31
解:如下建立jzys.m文件:
A=[3 5 7;6 1 2;-8 9 13]
B=[1 -4 6;2 9 4;7 3 -2]
C=[A B] %生成3*6矩阵
D=[A;B] %生成6*3矩阵
P=[C(2,1) C(2,4) C(2,5);C(3,1) C(3,4) C(3,5)] %提取C的第2和3行,1/4/5列
AA=[2 7 3 1;3 5 2 2;9 4 1 7] %系数矩阵
b=[6;4;2]
[m,n]=size(AA)
R=rank(AA) %求系数矩阵的秩
BB=[AA b] %增广矩阵
Rr=rank(BB) %求增广矩阵的秩
format rat
if R==Rr&R==n %n为未知数的个数,判断是否有唯一解
x=AA\b %唯一解
elseif R==Rr&R<n %判断是否有无穷解
x=AA\b %求特解
C=null(AA,'r') %求AX=0的基础解系,所得C为n-R列矩阵,这n-R列即为对应的基础解系
%这种情形方程组通解xx=k(p)*C(:,P)(p=1…n-R)
else X='No solution!' %判断是否无解
end
X='K1*C(:,1)+K2*C(:,2)+x' %非齐次线性方程组的通解
-----------------------------------------------------------------------------------
运行结果:
A = 3 5 7
6 1 2
-8 9 13
B = 1 -4 6
2 9 4
7 3 -2
C = 3 5 7 1 -4 6
6 1 2 2 9 4
-8 9 13 7 3 -2
D = 3 5 7
6 1 2
-8 9 13
1 -4 6
2 9 4
7 3 -2
P = 6 2 9
-8 7 3
AA =2 7 3 1
3 5 2 2
9 4 1 7
b = 6
4
2
m = 3
n = 4
R = 2
BB = 2 7 3 1 6
3 5 2 2 4
9 4 1 7 2
Rr =2
x = -2/11
10/11
0
0
C = 1/11 -9/11
-5/11 1/11
1 0
0 1
X =K1*C(:,1)+K2*C(:,2)+x本回答被提问者采纳
第3个回答 2010-12-31
第一问中
C为3*6矩阵,D为6*3矩阵,E为提取C的第二三行第一列,F为提取D的第二三行四五列,G就是想要的矩阵。
A=[3 5 7;6 1 2;-8 9 13];
B=[1 -4 6;2 9 4;7 3 -2];
C=[A B];
D=[A;B];
E=C(2:3,1);
F=C(2:3,4:5);
G=[E F]
第二:
解方程,把x4设为a,那么编程建立函数jfcz
function jfcz(A,B)
X=inv(A)*B
不过这个程序还要输入A矩阵还有B矩阵,求出X就是通解啦。
或者不建立函数直接进行
syms a;
A=[2 7 3;3 5 2;9 4 1];
B=[6-a;4-2*a;2-7*a];
X=inv(A)*B
这个程序直接得到X极为结果本回答被网友采纳