应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 大神看一下这个程序有什么问题
321/4返回列表1234下一页
查看: 1717  |  回复: 31
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

0000_8888

金虫 (著名写手)

[求助] 大神看一下这个程序有什么问题 已有3人参与

format long
A=[10 -1 -2;-1 10 -2;-1 -1 5];
b=[7.2;8.3;4.2];
omega=1.15;X0=[0;0;0;0];eps=1e-5;N=300;
x=SOR(A,b,omega,X0,eps,N);


function x=SOR(A,b,omega,x0,eps,N)
n=length(x0);x=ones(n,1);k=0;
while k<=N
    for i=1:n
        x(i)=omega*(- A(i,1:i-1)*x(1:i-1)...
            - A(i,i+1:n)*x0(i+1:n)+b(i))/A(i,i)+(1-omega)*x0(i);
    end
    k=k+1;
    if norm(x-x0,inf)<eps
        break;
    end
    x0=x;
end
if k>N
    warning('算法超出最大迭代次数！');
else
    disp(['迭代次数=  ',num2str(k)])
    x
end


运行结果是：
Index exceeds matrix dimensions.

Error in SOR (line 5)
        x(i)=omega*(- A(i,1:i-1)*x(1:i-1)...

不知道出现什么问题请指教！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
没有做不到，只有想不到
1楼 2014-11-05 13:58:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

linccnu

木虫 (小有名气)

【答案】应助回帖

第一次循环就通不过！
当i=1时，-A(i,i-1)*...----->-A(1,0)*....矩阵的子标是不允许为零和负数的！
赞一下(1人) 回复此楼
youcandoit!
25楼2014-11-09 15:17:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

readytogo
★
0000_8888(金币+1): 谢谢参与
2楼2014-11-05 14:22:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dmbb
★
0000_8888(金币+1): 谢谢参与
3楼2014-11-05 14:28:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qqaoqq
★
0000_8888(金币+1): 谢谢参与
4楼2014-11-05 14:38:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kiss20085537
★
0000_8888(金币+1): 谢谢参与
5楼2014-11-05 14:43:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

假大空
★
0000_8888(金币+1): 谢谢参与
6楼2014-11-05 14:44:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qqllqq
★
0000_8888(金币+1): 谢谢参与
7楼2014-11-05 14:50:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuxiuling

木虫 (知名作家)
★
0000_8888(金币+1): 谢谢参与
加油
回复此楼
jiayou
8楼2014-11-05 14:53:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dpabzh
★
0000_8888(金币+1): 谢谢参与
9楼2014-11-05 15:00:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

54smile晓
10楼2014-11-05 18:22:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 0000_8888 的主题更新
321/4返回列表1234下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今年的国基金是打分制吗？ 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 本人考085602 化学工程 专硕 +12 不知道叫什么！ 2026-03-15 14/700 2026-03-16 16:45 by 我的船我的海
[考研] 0703化学调剂，求各位老师收留 +8 秋有木北 2026-03-14 8/400 2026-03-16 15:21 by 哦哦123
[考研] 308求调剂 +3 是Lupa啊 2026-03-16 3/150 2026-03-16 10:07 by 求调剂zz
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 328求调剂 +3 5201314Lsy！ 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +3 Losir 2026-03-12 3/150 2026-03-14 12:11 by 热情沙漠
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 304求调剂（085602一志愿985） +12 化工人999 2026-03-09 12/600 2026-03-13 12:02 by JourneyLucky
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考研] 一志愿：武汉理工，材料工程，英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
关闭
请填处理意见
关闭