施密特正交化程序
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
function P = shimite(A)
[M,N] = size(A);
REF = rref(A); %%判断向量组是否线性相关
if(REF(M,:) == zeros(1,N))
fprintf('向量组线性相关,无法正交化!!!');
return;
end
%%施密特正交化
P = ones(M,1); %%正交化
A = double(A);
P(:,1) = A(:,1);
for i = 2:N
k = A(:,i);
for j = 1:i-1
k = k - (sum(A(:,i).*P(:,j))/sum(P(:,j).^2))*P(:,j);
end
P = [P,k];
end
for i = 1:N %%单位化
P(:,i) = P(:,i)/sqrt(sum(P(:,i).^2));
end
P;