24小时热门版块排行榜

返回列表

jydtjm

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 97
帖子: 16
在线: 4.1小时
虫号: 3763602
注册: 2015-03-25
专业: 机械摩擦学与表面技术

[求助] 帮忙看一下这段程序的思路是什么

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
帮忙看一下这段程序的思路是什么

回复此楼

» 猜你喜欢

1楼 2015-04-10 15:49:51

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

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

大家帮忙看一下我的编辑回信已经有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人回复

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 jydtjm 的主题更新

返回列表