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

最新推荐文章于 2024-07-10 22:11:32 发布

weixin_34275734

最新推荐文章于 2024-07-10 22:11:32 发布

阅读量1.9k

点赞数 1

文章标签： matlab python markdown

原文链接：https://my.oschina.net/zhichengjiu/blog/3013335

版权

2019独角兽企业重金招聘Python工程师标准>>>

    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

>>