matlab-线性代数 施密特正交化

     matlab : R2018a 64bit
      OS : Windows 10 x64
typesetting : Markdown
       blog : my.oschina.net/zhichengjiu
      gitee : gitee.com/zhichengjiu

 

code

clear
clc

% 施密特正交化
% 这些代码是本身写的，不具有权威性，只能作个简单的参考。
% 写这份代码的时候，大约是2016年。如今2019了，代码的细节都忘记了。
% 嗨...逆水行舟，不进则退呀。

a=ones(4,2);   % 四行，两列
% 必须是线性无关的向量组
a(:,1)=[1;1;2;3];
a(:,2)=[-1;1;4;-1];
 
a  % 看看a
a1 = a(:,1)
c=size(a1);  % 去判断a1多少行
 
b=zeros(c(1),2);    
 
b(:,1)=a1;   % 根据施密特，b1=a1
 
for i=2:1:2   % 须要改的是最后一个2，要把它改为共有几个向量的数目，好比说有a3，那么就改为3   
   for j=1:1:i-1 
   b(:,i)=b(:,i)+(-1*dot(b(:,j),a(:,2))/dot(b(:,j),b(:,j)))*b(:,j);
   end
end
 
n=input('想看第几个正交向量组的向量呢：')
b(:,n)

 

result

a =

     1    -1
     1     1
     2     4
     3    -1


a1 =

     1
     1
     2
     3

想看第几个正交向量组的向量呢：2

n =

     2


ans =

   -0.3333
   -0.3333
   -0.6667
   -1.0000

>>

 

resource

• [文档] ww2.mathworks.cn/help/matlab
• [文档] ww2.mathworks.cn/help/simulink
• [平台] www.oschina.net
• [平台] gitee.com

 

---

感谢帮助 志成就 的人们。
matlab优秀，值得学习。基础知识 + 专业知识 + matlab = ?
Simulink，用于仿真和基于模型的设计，值得学习。
该博文仅可用于测试与参考。