24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1269  |  回复: 5

Aemika_z

新虫 (初入文坛)

[求助] 研一菜鸟求电容层析成像敏感场的仿真代码已有3人参与

在做ECT敏感场仿真,求Matlab的代码

[ 发自手机版 http://muchong.com/3g ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sxjeff

木虫 (小有名气)

【答案】应助回帖

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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oeeker

新虫 (初入文坛)

3楼2015-11-26 03:21:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangdl123

新虫 (初入文坛)

【答案】应助回帖

jjdg: 直接pm最有效 2017-03-28 00:00:37
引用回帖:
2楼: Originally posted by sxjeff at 2014-12-22 09:56:11
function =landweber(CirS,Z,n_meas,n_elem)
%%%%%%%%初始化%%%%%%%%%
%m_fLandFactR1=0.05
m_fLandFactR1=1.5473e-009;
Trans_CirS=CirS';%Trans_CirS为灵敏度矩阵的转置
%%%%%%%%%%%迭代过程%%%%%
gk=m_fLa ...

能交流交流吗。。。我是一个初学者

发自小木虫IOS客户端
4楼2017-03-27 19:46:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Arthur_zyp

新虫 (初入文坛)

【答案】应助回帖

最近也在考虑做不做这个方向,感觉资料很少而且不太好做啊!!!
5楼2017-12-24 16:31:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leilzy

新虫 (初入文坛)

6楼2018-05-17 13:06:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Aemika_z 的主题更新
信息提示
请填处理意见