matlab求零空间,有关线性代数的Matlab代码笔记（2）行空间、零空间

今天继续，尝试加入一些范例依然是简单的内容： %%%%%%%%%%%%%%%%%%%%% 说明 %%%%%%%%%%%%%%%%%%% 行空间的基：按行的角度来看待矩阵，更多介绍在代码说明里，简单的利用了昨天的代码。 %%%%%%%%%%%%%%%%%%%%% CODE %%%%%%%%%%%%%%%%%% function B = rowba

今天继续，尝试加入一些范例

依然是简单的内容：

%%%%%%%%%%%%%%%%%%%%%说明 %%%%%%%%%%%%%%%%%%%

行空间的基：按行的角度来看待矩阵，更多介绍在代码说明里，简单的利用了昨天的代码。

%%%%%%%%%%%%%%%%%%%%%CODE %%%%%%%%%%%%%%%%%%

function B = rowbasis(A)

% rowbasis 行空间的基

% B = rowbasis(A) 返回A的行空间的一组基

% 以B的列表示

% A的行空间相当于A的转置的列空间

% rowbasis 找出A'中线性无关的

% 前r列，r代表A的秩

B = colbasis(A');

%%%%%%%%%%%%%%%%%%%%%范例%%%%%%%%%%%%%%%%%%

>> A=[1,1,1;1,3,2;2,3,4]

A =

1 1 1

1 3 2

2 3 4

>> rowbasis(A)

ans =

1 1 2

1 3 3

1 2 4

%%%%%%%%%%%%%%%%%%%%%说明 %%%%%%%%%%%%%%%%%%%

零空间的基：这个零空间的代码值得好好看看，体现了线性代数里关于列意义的基本思想。

直接看去你可能会不知道这些代码在干什么。因为一个简单的语句可能包含了很多操作。

我自己不太明白的时候是喜欢举个例子，把代码一句一句执行来看看的。这里要注意的是

零空间里的列指示了矩阵里构成零向量的列的线性组合。而代码里用到了很多构成矩阵的办法。

%%%%%%%%%%%%%%%%%%%%%CODE %%%%%%%%%%%%%%%%%%

function N = nulbasis(A)

% nulbasis 零空间的基.

% N = nulbasis(A) 在N的列中返回A的零空间的一组基

% 包含了Ax=0的n-r个特殊解

% freecol是自由列.

% Example:

% >> A = [1 2 0 3;

% [0 0 1 4];

% >> N = nulbasis(A)

% N = [-2 -3]

% [ 1 0]

% [ 0 -4]

% [ 0 1]

[R, pivcol] = rref(A, sqrt(eps));%精度sqrt(eps)=1.4901e-008

[m, n] = size(A);

r = length(pivcol);

freecol = 1:n;

freecol(pivcol) = [];%主元所在列

N = zeros(n, n-r);

N(freecol, : ) = eye(n-r);%由自由列序号给零空间的行按单位矩阵赋值

N(pivcol, : ) = -R(1:r, freecol);%主元所在列赋行简化阶梯矩阵对应值相反数

%这样组合能得到零向量

%这个方法接近我们手动求零空间的算法

本文原创发布php中文网，转载请注明出处，感谢您的尊重！

matlab求零空间,有关线性代数的Matlab代码笔记（2）行空间、零空间相关推荐

matlab fisher检验,FISHER线性判别MATLAB实现.doc
FISHER线性判别MATLAB实现 Fisher线性判别上机实验报告班级: 学号: 姓名: 算法描述 Fisher线性判别分析的基本思想:选择一个投影方向(线性变换,线性组合),将高维问题降低到一 ...
matlab求x对应y值,matlab 不知道函数表达式，已知y值求x
^根据给定的数来据,我们可以假自定函数表达式为baiy=b1+b2*x+b3*x^du2+b4*x^3+b5*x^4;所以上述函zhi数可以用matlab的regress()多元线性dao回归分析函数 ...
matlab求两方程交点坐标,matlab求交点坐标
b xi ? 0(i ? 1, 2, , n) 三.线性规划问题的求解方法二元线性规划问题的图解法线性规划问题的理论解法线性规划问题的MATLAB软件解法线性规划问题的图解...... MATL ...
如何用matlab求线性方程的通解,利用Matlab求线性方程组的通解
教育科学从摆箭利用求线性方程组的通解简绍勇杜玲陈勇新余高等专科学校数学与信息科学系江西新余 [摘要 ]讨论线性方程组解的个数及求线性方程组的通解问题是线性代数中的常见问 ...
matlab求被21整除,用matlab求[100,999]之间能被21整除的数的个数,还有建立一个字符串向量,删除大写字母._...
怕散!|用户:的问题用 matlab求[100,999]之间能被21整除的数的个数,还有建立一个字符串向量,删除大写字母. Promise丶陌影|用户:问问回答: clc;clear; m=0; f ...
matlab求球的体积,【matlab计算不规则物体体积资讯】matlab计算不规则物体体积足球知识与常识 - 足球百科 - 599比分...
matlab,求不规则三维图形体积? 我咋看也是个光滑曲面呀,哪来的体积?当然假如xy面到曲面之间全被填充为实体,应该能用循环或积分计算每个xy单位面积上的曲顶柱体的体积,然后求和,我用循环尝试了下: ...
如何用matlab求线性方程的通解,用matlab求线性方程的解
用矩阵除法,求线性方程组的特解利用矩阵除法求线性方程组的形式为:AX=b,其中A为系数矩阵,b为方程组右侧的一列常数. 由AX=b,得X=A\b,在这种条件下,首先要保证A为满秩矩阵,否则无法求解. ...
matlab求系统全响应,利用matlab求LTI连续系统的响应
<利用matlab求LTI连续系统的响应>由会员分享,可在线阅读,更多相关<利用matlab求LTI连续系统的响应(3页珍藏版)>请在人人文库网上搜索. 1.实验3 利用mat ...
3Blue1Brown:“线性代数的本质”完整笔记
我最早系统地学习线性代数是在大二时候,当时特意选修了学校物理系开设的4学分的线代,大概也就是比我们自己专业的线代多了一章向量空间的内容,其实最后上完发现,整个课程内容还是偏向于计算,对线性代数的几何直 ...

matlab求零空间,有关线性代数的Matlab代码笔记（2）行空间、零空间

matlab求零空间,有关线性代数的Matlab代码笔记（2）行空间、零空间相关推荐

最新文章

热门文章