anitagyt 发表于 2006-5-5 00:32

[求助]ABS算法 ,请指点一、二

刚编的程序，还是很不完整的，希望借助大家的力量最终能够完成，请帮我看看，谢谢！！<BR>这是原题http://photo1.hexun.com/p/2006/0430/19389/b_0F6AC971285C52FD.jpg<BR>http://photo1.hexun.com/p/2006/0430/19389/b_8AF49A6F3354BE30.jpg<BR>http://photo1.hexun.com/p/2006/0430/19389/b_663B6D96EA1BF735.jpg程序如下:clear
x=input('输入x1的值')
h=input('输入h1的值')
a=input('输入方程的系数矩阵A')
b=input('输入方程的常熟矩阵B')
m=input('输入迭代次数')
i=1
iflag=0
syms xi t w p hi
t=a'*x-b%表示不出来一个带x矩阵和实数阵相减
s=h*a
while s~=0
z=input('输入Z')
break
while z'*h*a=0
z=input('重新输入Z')
break
end
p=h'*z
i=i+1
while i~=m
xi=x-(t/(a'*p))*p%倒数不会，t也表示不出来
end
%应该有一个循环，不知道怎么弄(好像还少一个x=xi)
fprintf('X是Ax=b的一个解'),x
w=1/(h*a)
h=h-h*a*w'*h
i=i+1
%不知道如何转到第2步
end
if t==0
xi=x
hi=h
iflag=iflag+1
while i<m
i=i+1
%转到第二步，不会
end
fprintf('X是Ax=b的一个解'),x
else iflag=-i
break
fprinf('方程组Ax=b不相容')
end
h=h
iflag=iflag+1
if i<m
i=i+1
else fprintf('X是Ax=b的一个解'),x
end
else iflag=-i
break
fprinf('方程组Ax=b不相容')
end
不论是什么批评指正，都欢迎，再次谢谢大家的热心帮助！！

[ 本帖最后由 lxq 于 2007-4-26 18:12 编辑 ]

happy 发表于 2006-5-10 15:48

回复：（anitagyt）[求助]ABS算法 ,请指点一、二

不错

gczljhyn 发表于 2006-6-9 19:15

偶也是做ABS算法的!

希望可以讨论!QQ:386079617,540263301!编过简单的Huang算法,隐式LU,LX算法,希望跟你讨论!
页: [1]
查看完整版本: [求助]ABS算法 ,请指点一、二