| 查看: 1545 | 回复: 5 | ||||
[求助]
研一菜鸟求电容层析成像敏感场的仿真代码 已有3人参与
|
|
在做ECT敏感场仿真，求Matlab的代码 [ 发自手机版 http://muchong.com/3g ] |
回复此楼» 收录本帖的淘帖专辑推荐
 图像重建算法相关 |
» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
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
20