| 查看: 1474 | 回复: 5 | |||
[求助]
研一菜鸟求电容层析成像敏感场的仿真代码 已有3人参与
|
|
在做ECT敏感场仿真，求Matlab的代码 [ 发自手机版 http://muchong.com/3g ] |
回复此楼» 收录本帖的淘帖专辑推荐
 图像重建算法相关 |
» 猜你喜欢
护理论文 晋升
已经有4人回复
-
西南大学考核制博士
已经有5人回复
-
求碳排放博导；方向是LCA、生命周期可持续发展以及碳排放
已经有8人回复
-
评审感受-评审感受-评审感受
已经有11人回复
-
26年申博自荐-计算机视觉
已经有5人回复
-
窗边初夏的小雨
已经有9人回复
-
导师各种操作恶心咋办
已经有12人回复
-
现在不知道怎么办，感觉很痛苦
已经有5人回复
-
生物学博士
已经有6人回复
-
2026年申博-电池方向
已经有11人回复
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
5楼2017-12-24 16:31:00
leilzy
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 23.2
- 帖子: 13
- 在线: 11.1小时
- 虫号: 7857770
- 注册: 2018-01-31
- 性别: MM
- 专业: 用于检测、分析、成像及治
6楼2018-05-17 13:06:29