| 查看: 1017 | 回复: 7 | |||
[求助]
初学者的困惑，希望大神解答 已有2人参与
|
|
数值计算课的老师给布置了一道用matlab求解的题目，有一点困惑还望解答 这是我在网上找的高斯赛德尔迭代方法的源程序： function [x1 k1]=gauss_seidel(A,b,e) n=max(size(A)); for i=1:n if A(i,i)==0 'eig can not be zero' return end end k=0; kend=50; r=1; x=zeros(n,1); while k<=kend && r>e x0=x; for i=1:n s=0; for j=1:i-1 s=s+A(i,j)*x(j); end for j=i+1:n s=s+A(i,j)*x0(j); end x(i)=b(i)/A(i,i)-s/A(i,i); end r=norm(x-x0,inf); k=k+1; end if k>kend 'die dai bu shou lian ,failure' return else 'qiu jie cheng gong' end x1=x; k1=k; 这是我要求的问题： e=0.000000001; for n=4:12 A=hilb(n); b=A(:,1); [x k]=gauss_seidel(A,b,e) end 为什么在源程序中我把x=zeros（n，1）改成x=ones（n，1）就会报错，错误： Error in gauss_seidel (line 2) n=max(size(A)); Output argument "x1" (and maybe others) not assigned during call to "gauss_seidel". 非常困惑，我看源程序没有规定x必须取0，为什么不能改呢？ |
回复此楼» 猜你喜欢
316求调剂
已经有6人回复
-
346求调剂[0856]
已经有7人回复
-
一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！
已经有3人回复
-
310求调剂
已经有3人回复
-
机械专硕299求调剂至材料
已经有4人回复
-
070300化学319求调剂
已经有7人回复
-
08工科 320总分 求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有12人回复
-
307求调剂
已经有3人回复
-
265求调剂
已经有3人回复
mygt_hit
专家顾问 (职业作家)
-
专家经验: +362 - 应助: 438 (硕士)
- 贵宾: 0.019
- 金币: 19571.1
- 散金: 5130
- 红花: 135
- 沙发: 3
- 帖子: 4938
- 在线: 990.5小时
- 虫号: 1489764
- 注册: 2011-11-13
- 性别: GG
- 专业: 结构工程
- 管辖: 土木建筑
2楼2015-10-24 12:55:08
lvweizhong
铜虫 (著名写手)
- 应助: 57 (初中生)
- 金币: 3181.3
- 红花: 10
- 沙发: 1
- 帖子: 2064
- 在线: 167.6小时
- 虫号: 3683385
- 注册: 2015-02-09
- 专业: 数论
3楼2015-10-24 14:52:09
4楼2015-10-25 15:55:07
lvweizhong
铜虫 (著名写手)
- 应助: 57 (初中生)
- 金币: 3181.3
- 红花: 10
- 沙发: 1
- 帖子: 2064
- 在线: 167.6小时
- 虫号: 3683385
- 注册: 2015-02-09
- 专业: 数论
5楼2015-10-25 16:01:59
6楼2015-10-25 16:08:13
7楼2015-10-25 16:11:11
mygt_hit
专家顾问 (职业作家)
-
专家经验: +362 - 应助: 438 (硕士)
- 贵宾: 0.019
- 金币: 19571.1
- 散金: 5130
- 红花: 135
- 沙发: 3
- 帖子: 4938
- 在线: 990.5小时
- 虫号: 1489764
- 注册: 2011-11-13
- 性别: GG
- 专业: 结构工程
- 管辖: 土木建筑
8楼2015-10-25 18:38:47
