| 查看: 1638 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
研一菜鸟求电容层析成像敏感场的仿真代码 已有3人参与
|
|||
|
在做ECT敏感场仿真，求Matlab的代码 [ 发自手机版 http://muchong.com/3g ] |
回复此楼» 收录本帖的淘帖专辑推荐
 图像重建算法相关 |
» 猜你喜欢
极限拉扯，如何抉择
已经有15人回复
-
申请科研项目纯靠实力能成功吗？
已经有8人回复
-
中科院海洋所招收海洋材料相关领域博士研究生及研究助理
已经有4人回复
-
40+中年人失业后，还能找到好工作吗？
已经有3人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
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
200