上学期学了数值分析,现在有空分享一下当时作业的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;