MATLAB中施密特正交化的实现

1、施密特正交化的概念

2、MATLAB程序

 1 function b=Schmidt_orthogonalization(a)
 2 [m,n] = size(a);
 3 if(m<n)
 4     error('行小于列，无法计算，请转置后重新输入');
 5     return
 6 end
 7 b=zeros(m,n);
 8 %正交化
 9 b(:,1)=a(:,1);
10 for i=2:n
11     for j=1:i-1
12         b(:,i)=b(:,i)-dot(a(:,i),b(:,j))/dot(b(:,j),b(:,j))*b(:,j);
13     end
14     b(:,i)=b(:,i)+a(:,i);
15 end
16
17 %单位化
18 % for k=1:n
19 %     b(:,k)=b(:,k)/norm(b(:,k));
20 % end

说明1：如果需要单位化，把单位化后的注释删除即可

说明2：输入a为列向量矩阵，输出b也为列向量矩阵

转载于:https://www.cnblogs.com/Martin-Soaring/p/9719555.html

MATLAB中施密特正交化的实现相关推荐

python/sympy计算施密特正交化向量
sympy的符号计算功能很强大,学习矩阵分析,重温了线性代数中施密特正交化的方法,正好可以用sympy解决一些计算问题.施密特正交化,也称 Gram-Schmidt 正交化过程 (Gram–Schmi ...
matlab格拉姆施密特,改进的格拉姆-施密特正交化（modified Gram-Schmidt Process）
最近在重新学习线性代数,学习的教材是MIT Gilbert Strang 教授的<INTRODUCTION TO LINEAR ALGEBRA>,在第4.4章节格拉姆-施密特正交化时,书中 ...
MATLAB施密特正交化
施密特正交化的矩阵写法(MATLAB) clear clc a=input('所需要正交化向量构成的矩阵'); [lklkl,oioi]=size(a); for i=1:oioi if i= ...
matlab施密特正交化,浅谈压缩感知（十九）：MP、OMP与施密特正交化
浅谈压缩感知(十九):MP.OMP与施密特正交化关于MP.OMP的相关算法与收敛证明,这里仅简单陈述算法流程及二者的不同之处. 主要内容: MP的算法流程及其MATLAB实现OMP的算法流程以及MA ...
施密特正交 matlab,matlab-线性代数施密特正交化
matlab : R2018a 64bit OS : Windows 10 x64 typesetting : Markdown blog : my.oschina.net/zhichengjiu c ...
matlab中help所有函数功能的英文翻译
doc funname 在帮助浏览器中打开帮助文档 help funname 在命令窗口打开帮助文档 helpbrowser 直接打开帮助浏览器 lookfor funname 搜索某个关键字相关函数 ...
matlab中的矩阵
我们知道,求解线性方程组是线性代数课程中的核心内容,而矩阵又在求解线性方程组的过程中扮演着举足轻重的角色.下面我们就利用科学计算软件MATLAB来演示如何使用矩阵,同时,也使学生对线性代数的认识更加理 ...
漫步线性代数十八——正交基和格拉姆-施密特正交化(下)
格拉姆-施密特声明:以后博主会把文章的pdf版本陆续发布到的网上,免费供大家下载正交基和格拉姆-施密特正交化假设我们有是是三个无关向量a,b,ca,b,c,如果他们是正交的,那么会多问题都变得容 ...
漫步线性代数十七——正交基和格拉姆-施密特正交化(上)
对于一个正交基,每个向量和其他所有向量垂直,坐标轴就是互相正交的.我们还可以进一步改善:每个向量除以它的长度得到单位向量,这样的话正交基变成了标准正交基: 16.如果 qTiqj={01i≠j,给出正 ...
施密特正交化_格拉姆-施密特正交化--QR分解法的来源
人类的数学发展到了今天,已经如浩瀚的宇宙般深不见底,不管时代怎么变迁,数学总是以超前于人类的科技水平而存在,如果说宇宙有尽头的话,那么宇宙的尽头也是数学的尽头,数学的停止就是世界的末日.所以我后悔当初 ...

MATLAB中施密特正交化的实现

MATLAB中施密特正交化的实现相关推荐

最新文章

热门文章