当前位置:
文档之家› 三种迭代法matlab程序 数值分析
三种迭代法matlab程序 数值分析
• %A是n维非奇异阵。%b是n维向量。%P是初值。
• %delta是误差界。%max1是给定的迭代最高次数。%X为所求的方程组AX=b的近似解。
• N=length(b);
• for k=1:max1
• for j=1:N
•
if j==1
•
X(1)=(b(1)-A(1,2:N)*P(2:N))/A(1,1);
• %A是n维非奇异阵。%b是n维向量。%P是初值。
• %delta是误差界。%max1是给定的迭代最高次数。%X为所求的方程组AX=b的近似解。
ቤተ መጻሕፍቲ ባይዱ• N=length(b);
• for k=1:max1
• for j=1:N
•
if j==1
•
X(1)=(b(1)-A(1,2:N)*P(2:N))/A(1,1);
•
elseif j==N
•
X(N)=(b(N)-A(N,1:N-1)*(X(1:N-1))')/A(N,N);
•
else
•
X(j)=(b(j)-A(j,1:j-1)*X(1:j-1)-A(j,j+1:N)*P(j+1:N))/A(j,j);
•
end
• end
• err=abs(norm(X'-P));
X(j)=(b(j)-A(j,[1:j-1,j+1:N])*P([1:j-1,j+1:N]))/A(j,j); end err=abs(norm(X'-P)); P=X'; if(err<delta)
break end end X=X';k,err;
塞德尔迭代法的迭代程序:
• function X=gseid(A,b,P,delta,max1)
雅可比迭代法的Matlab程序
给 定 初 始 值 X P0 , 用 雅 克 比 迭 代 法 求 解 线 性 方 程 组
AX b,并生成序列Pk ,求不超过误差界的近似解。
function X=jacobi(A,b,P,delta,max1) %A是n维非奇异阵。%b是n维向量。%P是初值。%delta是误差界。 %max1是给定的迭代最高次数。%X为所求的方程组AX=b的近似解。 N=length(b); for k=1:max1 for j=1:N
• P=X';
• if(err<delta)
•
break
• end
• end
• X=X';
• err,k
雅可比迭代法的Matlab程序
给 定 初 始 值 X P0 , 用 雅 克 比 迭 代 法 求 解 线 性 方 程 组
AX b,并生成序列Pk ,求不超过误差界的近似解。
function X=jacobi(A,b,P,delta,max1) %A是n维非奇异阵。%b是n维向量。%P是初值。%delta是误差界。 %max1是给定的迭代最高次数。%X为所求的方程组AX=b的近似解。 N=length(b); for k=1:max1 for j=1:N
X(j)=(b(j)-A(j,[1:j-1,j+1:N])*P([1:j-1,j+1:N]))/A(j,j); end err=abs(norm(X'-P)); P=X'; if(err<delta)
break end end X=X';k,err;
塞德尔迭代法的迭代程序:
• function X=gseid(A,b,P,delta,max1)
X(j)=(b(j)-A(j,[1:j-1,j+1:N])*P([1:j-1,j+1:N]))/A(j,j); end err=abs(norm(X'-P)); P=X'; if(err<delta)
break end end X=X';k,err;
塞德尔迭代法的迭代程序:
• function X=gseid(A,b,P,delta,max1)
•
elseif j==N
•
X(N)=(b(N)-A(N,1:N-1)*(X(1:N-1))')/A(N,N);
•
else
•
X(j)=(b(j)-A(j,1:j-1)*X(1:j-1)-A(j,j+1:N)*P(j+1:N))/A(j,j);
•
end
• end
• err=abs(norm(X'-P));
• P=X';
• if(err<delta)
•
break
• end
• end
• X=X';
• err,k
雅可比迭代法的Matlab程序
给 定 初 始 值 X P0 , 用 雅 克 比 迭 代 法 求 解 线 性 方 程 组
AX b,并生成序列Pk ,求不超过误差界的近似解。
function X=jacobi(A,b,P,delta,max1) %A是n维非奇异阵。%b是n维向量。%P是初值。%delta是误差界。 %max1是给定的迭代最高次数。%X为所求的方程组AX=b的近似解。 N=length(b); for k=1:max1 for j=1:N
•
elseif j==N
•
X(N)=(b(N)-A(N,1:N-1)*(X(1:N-1))')/A(N,N);
•
else
•
X(j)=(b(j)-A(j,1:j-1)*X(1:j-1)-A(j,j+1:N)*P(j+1:N))/A(j,j);
•
end
• end
• err=abs(norm(X'-P));
• %A是n维非奇异阵。%b是n维向量。%P是初值。
• %delta是误差界。%max1是给定的迭代最高次数。%X为所求的方程组AX=b的近似解。
• N=length(b);
• for k=1:max1
• for j=1:N
•
if j==1
•
X(1)=(b(1)-A(1,2:N)*P(2:N))/A(1,1);
• P=X';
• if(err<delta)
•
break
• end
• end
• X=X';
• err,k