QR分解解方程的MATLAB程序.pdf资源-CSDN文库

版权申诉

5 浏览量 2021-10-30 04:55:48 上传评论收藏 4KB PDF 举报

资源推荐

资源详情

资源评论

function [Q,R,x]=QR(A,b)

n=length(A);Q=eye(n);R=A;

for i=1:n-1

x=R(i:n,i);

if all(x(2:n-i+1)~=0)

x=x/norm(x,inf);delta=norm(x,2)*sign(x(1));

v=[x(1)+delta;x(2:n-i+1)];beta=2/norm(v,2)^2;

else

beta=0;v=0;

end

h=eye(n-i+1)-beta*v*v';

H=blkdiag(eye(i-1),h);

Q=Q*H;

R=H*R;

end

R=triu(R);y=Q'*b;

for i=n:-1:1

x(i)=(y(i)-R(i,i+1:n)*x(i+1:n))/R(i,i);

end

%注解： Ax=b，A=QR。

本内容试读结束，登录后可阅读更多

下载后可阅读完整内容，剩余0页未读，立即下载

内容反馈

版权申诉

jishuyh

粉丝: 0
资源: 7万+

最新资源

资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！点击此处反馈

feedback-tip