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，用於仿真和基於模型的設計，值得學習。
該博文僅可用於測試與參考。