24小时热门版块排行榜    

查看: 1466  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dollapym

铁虫 (正式写手)

[求助] 急急急!!!求助用BICG或者restarted GMRES 解线性方程组A x=b 的MATLAB code 已有1人参与

求大神帮忙啊 , 很急如果需要的话可以再加金币的!!
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

http://www.netlib.org
Netlib is a collection of mathematical software, papers, and databases.
9楼2016-09-08 14:52:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
dollapym: 金币+10, 有帮助 2016-09-07 10:31:01
dollapym: 金币+10, 有帮助 2016-09-07 10:37:42
https://www.mathworks.com/help/matlab/ref/bicg.html
CODE:
load west0479;
A = west0479;
b = full(sum(A,2));
tol = 1e-12;
maxit = 20;
[x0,fl0,rr0,it0,rv0] = bicg(A,b,tol,maxit);
semilogy(0:maxit,rv0/norm(b),'-o');
xlabel('Iteration number');
ylabel('Relative residual');

2楼2016-09-06 20:26:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dollapym

铁虫 (正式写手)

引用回帖:
2楼: Originally posted by FMStation at 2016-09-06 20:26:01
https://www.mathworks.com/help/matlab/ref/bicg.html

load west0479;
A = west0479;
b = full(sum(A,2));
tol = 1e-12;
maxit = 20;
= bicg(A,b,tol,maxit);
semilogy(0:maxit,rv0/norm(b),'-o');
xla ...

不是直接用函数bicg,是要写这个方法的代码

发自小木虫IOS客户端
3楼2016-09-06 20:55:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

【答案】应助回帖

http://www.aip.de/groups/soe/local/numres/bookcpdf/c2-7.pdf
void linbcg(unsigned long n, double b[], double x[], int itol, double tol,
int itmax, int *iter, double *err)
4楼2016-09-06 21:46:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见