| 查看: 284 | 回复: 0 | ||
[求助]
帮忙看一下这段程序的思路是什么
|
|
while resid_newton>1e-6 [d,FLAG,RELRES,ITER] = pcg(@(u)Coef(beta,Hmatrix,u,x0,Mx,Nx,Mzcol,p),f0,tol,maxit,[],[],[]); % [d,FLAG,RELRES,ITER] = pcg(@(u)Coef(beta,Hmatrix,u,x0,Mx,Nx,Mzcol,p),f0,tol,maxit,@(u)Preinverse(beta,Hmatrix,u,x0,Mx,Nx,Mzcol,p),[],[]); RELRES x1=x0+d; f1=-fxklu(beta,Hmatrix,x1,Mx,Nx,Mzcol,g,p); F1=norm(f1,'fro'); functionvalue=[functionvalue,F1]; %%%the norm of equations value at u1 resid_newton=F1/F0; resid_newton_value=[resid_newton_value,resid_newton]; ITER_set=[ITER_set;ITER]; recovered_x=reshape(real(x1(1:Mxcol)),Mx,Nx); % in1=(recovered_x<0); recovered_x(in1)=0.001; % in2=(recovered_x>1); recovered_x(in2)=1; Image_Re_Err_recovered=norm(recovered_x-X,'fro')/norm(X,'fro'); Image_Re_Err=[Image_Re_Err,Image_Re_Err_recovered]; fprintf('%3.0fth iteration: ITER1=%2.0f,ReErr=%6.4f, \n', iter,ITER,Image_Re_Err_recovered); % figure; % imagesc(recovered_x);axis image; axis off;colormap gray; x0=x1; f0=f1; %index=(x1(1:Mxcol) <0); x1(index)=0; iter=iter+1; %min(min(x1(N^2+1:end))) num=[num,iter]; if iter>5, break, end end 帮忙看一下这段程序的思路是什么 |
回复此楼» 猜你喜欢
各位基础老师们，看完中央《从严治党》的决定后，对学术界更加有信心了吧
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
找到一些相关的精华帖子，希望有用哦~
大家帮忙看一下我的编辑回信
已经有7人回复
- 请有经验的大神来看下这个水解反应
已经有5人回复
- 大家帮忙看一下投稿三周后系统出现Achiving completed on 24-Aug-2014是什么意思
已经有5人回复
- 麻烦帮忙看一下ICCSAI2014这个会议是否靠谱
已经有8人回复
- 大家帮忙看一下这句话是什么意思？是要做什么？
已经有5人回复
- 帮忙看一下这个扩散方程，谢谢！
已经有8人回复
- 请大家帮忙看一下这个SEM照片，是属于脆性断裂吗
已经有5人回复
- 这能收敛吗？高手帮忙看一下
已经有12人回复
- 求大神来帮助我,,,,关于C语言的.懂得快进来,,,
已经有7人回复
- 请帮忙看一下lammps中是否有这个势函数？
已经有3人回复
- 程序不会写怎么办呢？
已经有4人回复
- fortran编译中$符号的使用
已经有4人回复
- 请帮忙看一下这句英文怎么写比较好
已经有3人回复
- 投稿Scientific Reports，技术审查未通过，其中这一条请大家帮忙看一下是什么意思？
已经有3人回复
- 大家帮忙看一下这个电荷密度
已经有19人回复
- 帮忙看一下下面这个谱图是不是我要的？为什么11ppm左右有一个峰？
已经有8人回复
- 各位大侠，帮忙看一下这个氧化铝的吸附解吸曲线滞后环是什么型的啊，孔型咋样啊~~谢谢
已经有12人回复
- 各位帮忙看一下这样的两篇文章是否具有相似性
已经有3人回复
- 求助啊，SCI修稿后一月有余，询问编辑，给了一封回信，大家帮忙看一下是什么意思啊
已经有24人回复
- 哪位大神能帮忙总结下常见天然产物的预检识方法啊，膜拜~~~~
已经有3人回复
- 紧急求助，大家能不能帮忙提供个计算的思路
已经有19人回复
- 帮忙看一下编辑是什么意思（Materials Research Bulletin）
已经有12人回复
- 【求助】哪位能帮忙写一下DS中小分子和蛋白对接的方法啊？计算一下能量什么的？非常感
已经有12人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
