[求助]帮忙看下程序

gui951 · 发表于 2006-6-2 13:55

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

function x=gauss_eli(A,b) A=[6 2 1 -1 2 4 1 0 1 1 4 -1 -1 0 -1 3]; b=[6 1 5 -5]'; n=length(b);p=1:n;LU=A; for k=1:n [c,i]=max(abs(LU(k:n,k))); ik=i+k-1; if ik~=k m=p(k);p(k)=p(ik);p(ik)=m; ck=LU(k,:);LU(k,:)=LU(ik,:);LU(ik,:)=ck; end if k==n break;end LU(k+1:n,k)=LU(k+1:n,k)/LU(k,k); LU(k+1:n,k+1:n)=LU(k+1:n,k+1:n)-LU(k+1:n,k)*LU(k,k+1:n); end L=diag(ones(n,1))+tril(LU,-1); U=triu(LU); p; y(1)=b(p(1)); for i=2:n y(i)=b(p(i))-L(i,1:i-1)*y(1:i,-1)'; end x(n)=y(n)/U(n,n); for i=n-1:-1:1 x(i)=(y(i)-U(i,i+1:n)*x(i+1:n)')/U(i,i); end

[此贴子已经被cdwxg于2006-6-2 23:51:40编辑过]

gui951 · 发表于 2006-6-2 14:00

??? Index exceeds matrix dimensions.
Error in ==> D:\MATLAB6p5\work\file1.m On line 24 ==> y(i)=b(p(i))-L(i,1:i-1)*y(1:i,-1)';

gui951 · 发表于 2006-6-2 15:23

dddddddddddddddddddddddddd

lxq · 发表于 2006-6-2 18:10

y(i)=b(p(i))-L(i,1:i-1)*y(1:i,-1)'; 最后这个分号没什么作用啊

nochance · 发表于 2006-6-2 21:00

你程序里面的y数组，应该是个一维的吧。那你调用y（1：i，-1），我觉得可能有问题，而且为什么是-1列呢？见红色的地方。

assist · 发表于 2006-6-2 23:45

矩阵维数不对

nochance · 发表于 2006-6-3 08:36

y(i)=b(p(i))-L(i,1:i-1)*y(1:i,-1)'; 当时敲错了吧，你的意思可能是y(1:i-1)'

gui951 · 发表于 2006-6-3 13:06

nochance 同志说对了问题解决了我自己可是检查了N遍都查不出来那是相当的郁闷

[此贴子已经被aspen于2006-6-3 13:55:20编辑过]

gui951 · 发表于 2006-6-3 13:08

nochance 同志说对了问题解决了我自己可是检查了N遍都查不出来那是相当的郁闷

账号		自动登录	找回密码
密码			我要加入

[求助]帮忙看下程序

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

错误信息

ddddddddddddd

谢谢大家非常感谢

谢谢大家非常感谢

[求助]帮忙看下程序

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

错误信息

ddddddddddddd

谢谢大家 非常感谢

谢谢大家 非常感谢

谢谢大家非常感谢

谢谢大家非常感谢