24小时热门版块排行榜

>论坛更新日志 (4201)
>文献求助 (447)
>虫友互识 (361)
>导师招生 (357)
>考博 (195)
>博后之家 (176)
>硕博家园 (160)
>休闲灌水 (148)
>论文投稿 (129)
>基金申请 (92)
>考研 (69)
>招聘信息布告栏 (60)
>教师之家 (56)
>公派出国 (45)
>绿色求助(高悬赏) (43)
>找工作 (38)

返回列表

279663962

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 581.1
散金: 27
红花: 2
帖子: 123
在线: 74.7小时
虫号: 3414258
注册: 2014-09-13
专业: 自动化

[交流] 大神们，这程序什么意思啊已有1人参与

function CirS = JacobianECT(elex,eley);

% obtain the distributions
elecfieldx(:,:,1) = elex(1,:)';
elecfieldx(:,:,2) = elex(2,:)';
elecfieldx(:,:,3) = elex(3,:)';
elecfieldx(:,:,4) = elex(4,:)';
elecfieldx(:,:,5) = elex(5,:)';
elecfieldx(:,:,6) = elex(6,:)';
elecfieldx(:,:,7) = elex(7,:)';
elecfieldx(:,:,8) = elex(8,:)';
elecfieldx(:,:,9) = elex(9,:)';
elecfieldx(:,:,10) = elex(10,:)';
elecfieldx(:,:,11) = elex(11,:)';
elecfieldx(:,:,12) = elex(12,:)';
elecfieldx(:,:,13) = elex(13,:)';
elecfieldx(:,:,14) = elex(14,:)';
elecfieldx(:,:,15) = elex(15,:)';
elecfieldx(:,:,16) = elex(16,:)';

elecfieldy(:,:,1) = eley(1,:)';
elecfieldy(:,:,2) = eley(2,:)';
elecfieldy(:,:,3) = eley(3,:)';
elecfieldy(:,:,4) = eley(4,:)';
elecfieldy(:,:,5) = eley(5,:)';
elecfieldy(:,:,6) = eley(6,:)';
elecfieldy(:,:,7) = eley(7,:)';
elecfieldy(:,:,8) = eley(8,:)';
elecfieldy(:,:,9) = eley(9,:)';
elecfieldy(:,:,10) = eley(10,:)';
elecfieldy(:,:,11) = eley(11,:)';
elecfieldy(:,:,12) = eley(12,:)';
elecfieldy(:,:,13) = eley(13,:)';
elecfieldy(:,:,14) = eley(14,:)';
elecfieldy(:,:,15) = eley(15,:)';
elecfieldy(:,:,16) = eley(16,:)';

% calculate the Sensitivity
k = 1;
% for i = 1:16
%    for j = 1:16
%       CirS(k,:) = -(elecfieldx(:,:,i).*elecfieldx(:,:,j) + elecfieldy(:,:,i).*elecfieldy(:,:,j));
%       k=k+1;
%    end
% end

for i = 1:16
for j = i+2:i+14
      if j<17
         CirS(k,:) = -(elecfieldx(:,:,i).*elecfieldx(:,:,j) + elecfieldy(:,:,i).*elecfieldy(:,:,j));
         k=k+1;
      else
         CirS(k,:) = -(elecfieldx(:,:,i).*elecfieldx(:,:,j-16) + elecfieldy(:,:,i).*elecfieldy(:,:,j-16));
         k=k+1;
      end
end
end

回复此楼