| 查看: 1579 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
研一菜鸟求电容层析成像敏感场的仿真代码 已有3人参与
|
|||
|
在做ECT敏感场仿真，求Matlab的代码 [ 发自手机版 http://muchong.com/3g ] |
回复此楼» 收录本帖的淘帖专辑推荐
 图像重建算法相关 |
» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有12人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有16人回复
-
初秋的晨风
已经有3人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有11人回复
-
Ei源刊怎么投
已经有4人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有6人回复
-
两块石头
已经有5人回复
-
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生
已经有3人回复
-
找导师
已经有8人回复
5楼2017-12-24 16:31:00
sxjeff
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 2630.8
- 散金: 100
- 帖子: 172
- 在线: 316.3小时
- 虫号: 2167039
- 注册: 2012-12-05
- 性别: GG
- 专业: 机械测试理论与技术
【答案】应助回帖
|
function [gk,PkR]=landweber(CirS,Z,n_meas,n_elem) %%%%%%%%初始化%%%%%%%%% %m_fLandFactR1=0.05 m_fLandFactR1=1.5473e-009; Trans_CirS=CirS';%Trans_CirS为灵敏度矩阵的转置 %%%%%%%%%%%迭代过程%%%%% gk=m_fLandFactR1*Trans_CirS*Z;%zeros(812,1); for k=1:500 s1=gk-m_fLandFactR1*Trans_CirS*(CirS*gk-Z); gk=s1; % for i=1:812 % if gk(i)<0 % gk(i)=0; % end % if gk(i)>1 % gk(i)=1; % end % end end %%%%%%%%预迭代Pre-Iteration%%%%%%% A0R=Trans_CirS*m_fLandFactR1; M_IR=eye(n_meas); PkR=A0R; %%%%%%%%二阶迭代%%%%% % for i=0:1:50 % PkR=PkR*(M_IR+M_IR-CirS*PkR); % end %%%%%%%%一阶迭代%%%%% for i=1:500 PkR=PkR+A0R*(M_IR-CirS*PkR); end %%%%%%%%%%%%%%%%%%%% 我也在摸索当中，相互学习  |
2楼2014-12-22 09:56:11
3楼2015-11-26 03:21:03
4楼2017-03-27 19:46:34