clc;clear?all;close?all;
f?=?1e9;
c0?=?3e8;
lambda?=?c0/f;
step?=?lambda/10;
factor?=?1.5;
%%?sphere?創建目標三角剖分網格?
R?=?factor?*?lambda;
Jl?=?pi*R;
N_theta?=?round（Jl/step）;
theta?=?0?:?pi/N_theta?:?pi;
X?=?[];?Y?=?[];?Z?=?[];
for?n?=?1?:?length（theta）
????th?=?theta（n）;
????if?abs（th-0）<1e-4
????????X（end+1）?=?0;?Y（end+1）?=?0;?Z（end+1）?=?R;
????elseif?abs（th-pi）<1e-4
????????X（end+1）?=?0;?Y（end+1）?=?0;?Z（end+1）?=?-R;
????else
????????r?=?R*sin（th）;
????????Wl?=?2*pi*r;
????????if?Wl?<=?2*step
????????????X（end+1）?=?r*cos（0）;?Y（end+1）?=?r*sin（0）;?Z（end+1）?=?R*cos（th）;
????????????X（end+1）?=?r*cos（pi）;?Y（end+1）?=?r*sin（pi）;?Z（end+1）?=?R*cos（th）;
????????else
????????????N_phi?=?round（Wl/step）;
????????????phi?=?0?:?2*pi/N_phi?:?2*pi;
????????????for?m?=?1?:?length（phi）
????????????????X（end+1）?=?r*cos（phi（m））;?Y（end+1）?=?r*sin（phi（m））;?Z（end+1）?=?R*cos（th）;
????????????end
????????end
????end????
end
X?=?X（:）;?Y?=?Y（:）;?Z?=?Z（:）;
DT?=?delaunayTriangulation（XYZ）;
[tdatav]?=?convexHull（DT）;
pdata?=?[X?Y?Z];
trisurf（tdata?X?Y?Z‘FaceAlpha‘1.0）;axis?equal;box?on;
xlabel（‘x‘）;?ylabel（‘y‘）;?zlabel（‘z‘）;
p_whole?=?‘pdata_whole.txt‘;
t_whole?=?‘tdata_whole.txt‘;
save（p_whole‘pdata‘‘-ascii‘）;
save（t_whole‘tdata‘‘-ascii‘）;
%%?MoM計算目標RCS
RCS_plane?=?1;
dz?=?[-1.0?0.0?0.0];
Pol?=?[0.0?0.0?1.0];
fid?=?fopen（‘InputPara_MoM.txt‘‘w‘）;
fprintf（fidp_whole）;fprintf（fid‘?‘）;fprintf（fidt_whole）;fprintf（fid‘\r\n‘）;
fprintf（fid‘freq:?‘）;fprintf（fidnum2str（f））;fprintf（fid‘\r\n‘）;
fprintf（fid‘RCS_plane:?‘）;fprintf（fidnum2str（RCS_plane））;fprintf（fid‘\r\n‘）;
fprintf（fid‘dz:?‘）;fprintf（fidnum2str（dz））;fprintf（fid‘\r\n‘）;
fprintf（fid‘Pol:?‘）;fprintf（fidnum2str（Pol））;fprintf（fid‘\r\n‘）;
fclose（fid）;
system（‘MoM_EFIE.exe‘）;
%%?結果可視化
figure;
theta?=?-180:180;
rcs_mom?=?importdata（‘BiRCS_xoz_MoM.txt‘）;
plot（thetarcs_mom‘b‘）;hold?off;
grid?on;?xlabel（‘theta/deg‘）;ylabel（‘BiRCS/dBsm‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????1562112??2013-08-20?21:36??Matlab+cpp三維矩量法（MoM）通用計算程序\blas_win64_MT.dll

?????文件????8342016??2013-08-20?21:44??Matlab+cpp三維矩量法（MoM）通用計算程序\lapack_win64_MT.dll

?????文件????????432??2012-08-30?10:57??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\armadillo_icon.png

?????文件?????203838??2012-08-30?10:57??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\armadillo_nicta_2010.pdf

?????文件??????16074??2014-09-18?17:08??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\CMakeLists.txt

?????文件????????690??2011-03-22?13:51??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\cmake_aux\InstallFiles\ArmadilloConfig.cmake.in

?????文件????????381??2011-03-22?13:51??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\cmake_aux\InstallFiles\ArmadilloConfigVersion.cmake.in

?????文件???????1647??2009-06-23?23:04??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\cmake_aux\Modules\ARMA_CheckMathProto.cmake

?????文件???????1633??2009-06-23?23:04??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\cmake_aux\Modules\ARMA_CheckProto.cmake

?????文件???????1174??2014-09-16?21:58??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\cmake_aux\Modules\ARMA_FindACML.cmake

?????文件???????1268??2014-09-16?21:58??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\cmake_aux\Modules\ARMA_FindACMLMP.cmake

?????文件????????660??2014-09-16?21:58??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\cmake_aux\Modules\ARMA_FindARPACK.cmake

?????文件????????991??2014-09-16?21:58??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\cmake_aux\Modules\ARMA_FindBLAS.cmake

?????文件???????1374??2014-09-16?21:58??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\cmake_aux\Modules\ARMA_FindCBLAS.cmake

?????文件???????1512??2014-09-16?21:58??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\cmake_aux\Modules\ARMA_FindCLAPACK.cmake

?????文件???????1055??2014-09-16?21:58??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\cmake_aux\Modules\ARMA_FindLAPACK.cmake

?????文件???????1664??2013-04-18?12:50??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\cmake_aux\Modules\ARMA_FindMKL.cmake

?????文件???????1101??2014-09-16?21:58??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\cmake_aux\Modules\ARMA_FindOpenBLAS.cmake

?????文件????????425??2013-12-26?13:10??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\configure

?????文件?????366596??2014-11-06?10:38??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\docs.html

?????文件????????344??2015-01-06?20:11??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\examples\B.txt

?????文件???????4612??2014-10-27?15:27??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\examples\example1.cpp

?????文件????????143??2014-01-01?15:59??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\examples\example1_win64.README.txt

?????文件???17367040??2015-03-05?08:54??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\examples\example1_win64.sdf

?????文件???????1276??2014-01-01?15:08??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\examples\example1_win64.sln

????..A..H.?????16896??2015-03-05?08:54??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\examples\example1_win64.v11.suo

?????文件???????5809??2013-12-15?20:29??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\examples\example1_win64.vcxproj

?????文件???????1814??2015-01-06?20:07??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\examples\example2.cpp

?????文件???????1305??2010-01-24?22:14??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\examples\example_lsq.cpp

?????文件????1562112??2013-08-20?21:36??Matlab+cpp三維矩量法（MoM）通用計算程序\matrix_library\armadillo\armadillo-4.450.0\examples\lib_win64\blas_win64_MT.dll

............此處省略572個文件信息