400电话西安网站制作 彩铃制作,常德论坛最新消息,企业网站wap源码,上海建筑工程网#最近在学习matlab#xff0c;刚好和线代论文重合了
于是心血来潮用matlab建了一个模型来解决施密特正交规范化矩阵。
我们知道这个正交化矩阵挺公式化的#xff0c;一般公式化的内容我们都可以用计算机来进行操作#xff0c;节约我们人工的时间。
我们首先把矩阵导入进去…#最近在学习matlab刚好和线代论文重合了
于是心血来潮用matlab建了一个模型来解决施密特正交规范化矩阵。
我们知道这个正交化矩阵挺公式化的一般公式化的内容我们都可以用计算机来进行操作节约我们人工的时间。
我们首先把矩阵导入进去之后计算投影的长度进行缩减。进行正交化
之后单独进行规范化
代码如下
%%% 施密特正交化
aones(3,3); % 假设3行3列视题目而定
% 必须是线性无关的向量组%带入题目的向量得到
a(:,1)[1;2;-1];
a(:,2)[-1;3;1];
a(:,3)[4;-1;0];a1 a(:,1);
csize(a1); % 去判断a1多少行bzeros(c(1),3);%与题目中向量个数行列一样
b(:,1)a(:,1);
%正交化
for i2:3%从第二个开始递归for j1:i-1b(:,i)b(:,i)-dot(a(:,i),b(:,j))/dot(b(:,j),b(:,j))*b(:,j);endb(:,i)b(:,i)a(:,i);
end%单位化for k1:3b(:,k)b(:,k)/norm(b(:,k));enddisp(b);%b即为正交规范组 我们带入一道例题
设a1[1,2,-1]T,a2[-1,3,1]T,a3[4,-1,0]T 试用施密特正交化过程把这组向量规范正交化. 解取b1a1 b2a2-[a2,b1]/(||b1||)2b1[-1,3,1]T-4/6[1,2,-1]T5/3[-1,1,1]T; B3a3-[a3,b1]/(||b1||)2b1-[a3,b2]/(||b2||)2b2[4,-1,0]T-1/3[1,2,-1]T5/3[-1,1,1]T2[1,0,1]T; 再把它们单位化
取e1b1/||b1||1/√6[1,2,-1]Te2b2/||b2||1/√3[-1,1,1]Te3b3/||b3||1/√2[101]T 以上为数学人工方法求解而如果直接带入我们建立的模型我们可以直接得到 谢谢大家啦以上就是代码和展示部分
