| 查看: 634 | 回复: 2 | |||
[交流]
[交流贴]分享我的源代码（很短的代码实现矩量法计算圆柱散射） 已有1人参与
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%% % TM波： clear; N = 128; Z = zeros(N); ka = 2 * 2 * pi; kl = ka / N * 2 * pi; phi = (0 : N - 1) * 2 * pi / N; Z(1, 1) = 1 - 1i * 2 / pi * (log(kl / 2) - 1); Z(1, 2 : N) = besselh(0, 2, ka * 2 * abs(sin(phi(2 : N) / 2))); for nn = 2 : N Z(nn,  = circshift(Z(1, , [0, nn - 1]);end phiAngle = 360; phiAngle = (0 : phiAngle - 1) * 2 * pi / phiAngle; Esca = exp(1i * ka * cos(bsxfun(@minus, phiAngle', phi))); Einc = exp(1i * ka * cos(phi')); rcs = 4 / (pi * ka) * abs(Esca * (Z \ Einc)) .^ 2; hold on; plot(phiAngle * 180 / pi, 10 * log10(rcs), 'r', 'linewidth', 2); %%%%%%%%%%%%%%%%%%%%%%%%%%%% %TE波 clear; N = 128; Z = zeros(N); ka = 2 * 2 * pi; kl = ka / N * 2 * pi; phi = (0 : N - 1) * 2 * pi / N; Z(1, 1) = 0.5; Z(1, 2 : N) = abs(sin(phi(2 : N) / 2)) .* besselh(1, 2, ka * 2 * abs(sin(phi(2 : N) / 2))) * kl / 4i; for nn = 2 : N Z(nn, = circshift(Z(1, , [0, nn - 1]);end phiAngle = 360; phiAngle = (0 : phiAngle - 1) * 2 * pi / phiAngle; temp = cos(bsxfun(@plus, phiAngle', phi)); Hsca = temp .* exp(1i * ka * temp); Hinc = exp(1i * ka * cos(phi')); rcs = kl ^ 2 / (4 * pi * ka) * abs(Hsca * (Z \ Hinc)) .^ 2; hold on; plot(phiAngle * 180 / pi, 10 * log10(rcs)); %%%%%%%%%%%%%%%%%%%%%%%%%%%% 希望能对需要的人有所帮助，尤其入门搞矩量法的   |
回复此楼» 猜你喜欢
参与限项
已经有5人回复
-
推荐一本书
已经有7人回复
-
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【分享】大学爆笑文集 折腾岁月 第一部 第二部《不断更新！》【已搜索无重复】
已经有18人回复
- 分享一篇SCI文章历程
已经有160人回复
- 【分享】复合材料结构力学 英文版【已搜索无重复】
已经有22人回复
- 【分享】操作系统GUI进化史【已搜索无重复】
已经有4人回复
- 【分享】基于OpenCV的智能视频监控设计
已经有11人回复
- 【分享】态密度处理程序
已经有17人回复
- 【分享】C++程序设计总结——学C++的朋友多注意！！！
已经有10人回复
- 【分享】Python核心编程(第二版)(pdf)【已搜索无重复】
已经有5人回复
2楼2011-08-20 12:22:56
3楼2016-05-20 11:19:38