| 查看: 1547 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
研一菜鸟求电容层析成像敏感场的仿真代码 已有3人参与
|
|||
|
在做ECT敏感场仿真，求Matlab的代码 [ 发自手机版 http://muchong.com/3g ] |
回复此楼» 收录本帖的淘帖专辑推荐
 图像重建算法相关 |
» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有7人回复
-
这样的filecode谁见过
已经有18人回复
-
关于代码变化问题，想知道的进来
已经有17人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
leilzy
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 23.2
- 帖子: 13
- 在线: 11.1小时
- 虫号: 7857770
- 注册: 2018-01-31
- 性别: MM
- 专业: 用于检测、分析、成像及治
6楼2018-05-17 13:06:29
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
20