matlab中最大无关组,matlab如何求极大线性无关组

匿名用户

3级

2013-04-05 回答

调出实验中的矩阵A、B

1.作出A的行向量组：a1、a2、a3、a4、a5、a6；

2.作出B的列向量组：b1、b2、b3、b4、b5、b6；

3.由A的一、三、五行，二、三、四列交叉点上的元素作出子矩阵A3；

4.做一个10阶矩阵A4，其分块形式为A4 ；

5.由索引向量L产生取A的第2.、4、5行所成的子矩阵A5；

6.将A 对应的行向量组正交规范为正交向量组A6，并验证所得的结果；

7.求a1与a2的内积A7；

8.完成以下初等变换：将A 的一、四行互换，再将其第三列乘以，6再将其第一行的10倍加至第五行；

9.求B的列向量组的一个极大线性无关向量组A9，并将其余向量用极大线性无关向量组线性表示

给一个例子：

clc; clear all;

%% 矩阵预处理

A = randint(4, 5, [10 100]);

r = rank(A);

s1 = size(A);

c = 1 : min(size(A));

combos = combntns(c,r);

s2 = size(combos);

%% 求出极大线性无关组

if s1(1) <= s1(2)

for i = 1 : s2(1)

B(:, :, i) = A(combos(i, 1:r), 1:s1(2));

if rank(B(:, :, i)) == r

C = B(:, :, i);

end

else

for i = 1 : s2(1)

B(:, :, i) = A(1 : s1(1), combos(i, 1:r));

if rank(B(:, :, i)) == r

C = B(:, :, i);

end

%% 正规化

for i = 1 : size(C, 1)

a1 = max(C(i, :));

for j = 1 : size(C, 2)

a1 = gcd(a1, C(i, j));

end

C(i, :) = C(i, :) / a1;

end

%% 打印

fprintf('\n矩阵: \n')

disp(A);

fprintf('基向量为: \n\n');

for i = 1 : size(C, 1)

if i > 1 & i <= size(C, 1)

fprintf(' + \n');

end

str = sprintf('k%d*[', i);

fprintf(str);

for j = 1 : size(C, 2)-1

fprintf('%.3f, ', C(i, j));

end

fprintf('%.3f ]', C(i, j+1));

end

fprintf('\n其中，ki为任意实数 \n');

matlab中最大无关组,matlab如何求极大线性无关组相关推荐

Matlab 矩阵分析（一）矩阵的秩，极大线性无关组和基
Matlab 矩阵分析(一) 矩阵的秩和极大线性无关组 r = rank(A) 计算矩阵的秩 A 是矩阵 [R,j]=rref(A) 计算行最简 R是矩阵的行最简矩阵,j是矩阵的极大线性无关组的列数 ...
matlab中trapz的用,matlab使用trapz函数对积分的运算
Matlab中find函数使用输出是对的.因为他找到的是B(2:n,2)中的第一个元素.你可以好好看看find函数的作用.你如果写下面代码应该就对了:A=B(2:n,2);A(q=find(B(2: ...
matlab中dist的命令,matlab dist函数
dist--欧式距离加权函数(Euclidean distance weight function) 语法: Z = dist(W,P) df = dist('deriv') D = dist(pos ...
matlab中错误使用fmincon,MATLAB中fmincon 函数问题
MATLAB中fmincon 函数问题 Matlab的fmincon优化问题请问: 各位高手帮忙看看我的程序又什么问题?显示错误 Error in ==> Fun at 33 [w,fval] ...
向量组秩及其极大线性无关组求解浅析
向量组秩和极大线性无关组求解问题来源阐述线性代数课程中,在学习了向量组的线性相关性和向量组的秩后,一类常见的计算问题是给出向量组,求解其秩和极大线性无关组.课程中一般给出的方法都是以向量为列组成矩 ...
matlab中常微分方法,MATLAB解常微分方程组的解法（好东西要共享）
1:问题常微分方程的初值问题的标准数学表述为:y'=f(t,y),a<=t<=b,y(a)=y(0) :我们要求解的任何高阶常微分方程都可以用替换法化为上式所示的一阶形式,其中y为向量, ...
matlab中ode指令,[转载]MATLAB中ODE的使用
ode23 解非刚性微分方程,低精度,使用Runge-Kutta法的二三阶算法.ode45 解非刚性微分方程,中等精度,使用Runge-Kutta法的四五阶算法.ode113 解非刚性微分方程,变精度 ...
matlab中 trapz函数,trapz(matlab中trapz函数用法)
采用梯形法求积分下面两种用法等效 h=pi/200; t=0:h:pi/2; y=sin(t)+cos(2*t); z1=trapz(t,y) z2=trapz(y)*h z1 = 1.0000 z ...
matlab中有趣的图像,MATLAB中有趣的演示实例
MATLAB中有趣的演示实例 MATLAB中也同样有许多有意思的实例,为提高读者对MATLAB和Simulink的兴趣,特举部分以供参考.具体如下,运行的时候只要将":"前面的代码 ...

matlab中最大无关组,matlab如何求极大线性无关组

matlab中最大无关组,matlab如何求极大线性无关组相关推荐

最新文章

热门文章