应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 大神们，这程序什么意思啊
51/1返回列表
查看: 493  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

279663962

金虫 (小有名气)

[交流] 大神们，这程序什么意思啊 已有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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-05-23 21:49:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

279663962

金虫 (小有名气)
引用回帖:
2楼: Originally posted by zhangzhiguang at 2015-05-25 11:56:30
这应该是计算或构造CirS数组或矩阵。将获得的elex,eley这样的二维数组变成三维数组的一部分，然后通过下面的运算得到CirS数组。

牛
3楼2015-06-30 16:18:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

zhangzhiguang

至尊木虫 (文坛精英)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
这应该是计算或构造CirS数组或矩阵。将获得的elex,eley这样的二维数组变成三维数组的一部分，然后通过下面的运算得到CirS数组。
赞一下(1人) 回复此楼
好好学习,天天向上!
2楼2015-05-25 11:56:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 06:20 by xzcv0Tehd8lQ
[基金申请] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 10/500 2026-11-02 06:16 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:13 by xzcv0Tehd8lQ
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 5/250 2026-11-02 06:07 by xzcv0Tehd8lQ
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:07 by xzcv0Tehd8lQ
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-02 05:41 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 21/1050 2026-11-02 05:41 by xzcv0Tehd8lQ
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-02 05:37 by xzcv0Tehd8lQ
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +5 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-02 05:34 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 03:26 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:47 by xzcv0Tehd8lQ
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:01 by xzcv0Tehd8lQ
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 21:45 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 21:45 by xzcv0Tehd8lQ
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
关闭
请填处理意见
关闭