Question

MATLAB数学实验。高分悬赏答案!!!!!!!!!!!!!!!!!!!!!!!

A=[357;6 1 2;-8 9 13]
B=[1 -4 6;2 9 4;7 3 -2]建立jzys.m文件求解下列各题
2） 将A B 扩展为 3*6的矩阵 与 的矩阵6*3的距阵D
3） 提取C 得第2，3行，1，4，5列元素构成子阵 P
4)建立jfcz.m文件，求方程组
2X1+7X2+3X3+X4=6
3X1+5X2+2X3+2X4=4
9X1+4X2+X3+7X4=2的通解
[问题解答]
1．解：如下建立jzys.m文件：
（matlab代码）
运行结果为：（把在命令窗口中得到的结果复制粘贴过来，注意去掉多余的空格）
不好意思
是
） 将A B 扩展为 3*6的矩阵 C与 的矩阵6*3的距阵D
3） 提取C 得第2，3行，1，4，5列元素构成子阵 P

其余的是对的 麻烦高人给解答以下啊

麻烦帮我解一下第2问吧 第2大体 通解的~~!!!

Answer 1

1.先用下面这段代码找到方程的一个特解
clc
clear all
format long

A=[2 7 3 1;3 5 2 2;9 4 1 7];
B=[6;4;2];
format rat;
X=A\B
解为
X =

-2/11
10/11
0
0
因为矩阵中有两个非零元素，所以矩阵A的秩为2，基本解有n-rank(A)个
2．求方程组的基本解
Z=null(A,'r')
解为
Z =

1/11 -9/11
-5/11 1/11
1 0
0 1
所以该方程组的通解为
X通解=X+Z1*q1+Z2*q2
q1,q2为任意向量
完毕

Answer 2

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 %通解

Answer 3

解：如下建立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

Answer 4

第一问中
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极为结果