反幂法

上学期学了数值分析,现在有空分享一下当时作业的matlab代码,防止以后遗忘。这只是具体题目中的一个例子,可以参考

clc;
format long
A=[-3 1 0;1 -3 -3;0 -3 4];
[x,y]=fanmi(A,5,0.0001)
子函数:
function [w,y]=fanmi(A,p,es,maxit)
if nargin<2,error('at last 2 input arguments required'),end
if nargin<4|isempty(maxit),t=50;end
if nargin<3|isempty(es),es=0.00001;end
[m,n]=size(A);
if m~=n,error('Matrix A must be square');end
u=zeros(n,1);
y=zeros(n,1);
for c=1:n
    a(c)=p;
end
B=diag(a);
u(1)=1;%设定u的初值
t=0;   %记录迭代次数
while(1)
  for j=1:n
    v=sqrt(u'*u);
    y(j)=u(j)/v;
  end
    t=t+1;
    if(t<2)
        w=y'*u;
    end
    u=(A-B)\y;  %求解Au=y的方程
    w0=w;
    w=y'*u;
    ea=abs(w^(-1)-w0^(-1))/abs(w^(-1));
    if(ea<=es)
        m
        break
    end
end
w=1/w+p;

 

  • 2
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值