С просьбой

[fatma1000]

это MATLAB кода для upml
Мне нужно преобразовать его в C
D = 8;% глубина PML регионе (D количества клеток)
NXA% = 100 х 50 = NYA Original расчетной области
NXA = 100; NYA = 50;
NX NXA = 2 * D; аспекта% Y
NY NYA = 2 * D;% X размерность
EY (1: NX, 1: Нью-Йорк) = 0; компонент% EY
DY (1: NX, 1: Нью-Йорк) = 0; компонент% Dy
DY1 (1: NX, 1: Нью-Йорк) = 0;% предыдущем шаге
EX (1: NX, 1: Нью-Йорк) = 0;% Ex компонента
DX (1: NX, 1: Нью-Йорк) = 0;% Dx компонента
DX1 (1: NX, 1: Нью-Йорк) = 0;% предыдущем шаге
HZ (NX, NY) = 0; компонент% Гц
BZ (NX, NY) = 0; компонент% Bz
BZ1 (NX, NY) = 0;% предыдущем шагеIS = NX / 2; JS = NY / 2;CC = 2.99792458e8;m0=4.0*pi*1.0e-7; %permability of free space E0 = 1,0 / (CC * CC * M0)% проницаемость свободного пространства
DX = 1.5E-2; Dy = 1.5E-2, размер% клетке
DT = DX / (2,0 * CC); шагом% времени размера

% C = 1/sqrt (M0 * E0);% скорости света
N = 700;% число шагов по времени
H (N) = 0; функцию возбуждения%
W (N) = 0;% смотрят точка
% Матрицами коэффициентов
cadx (1: NX, 1: Нью-Йорк) = 1;
cbdx (1: NX, 1: Нью-Йорк) = DT;
Кади (1: NX, 1: Нью-Йорк) = 1;
cbdy (1: NX, 1: Нью-Йорк) = DT;
caex (1: NX, 1: Нью-Йорк) = 1;
cbex (1: NX, 1: Нью-Йорк) = 1/e0;
ccex (1: NX, 1: Нью-Йорк) = 1/e0;
caey (1: NX, 1: Нью-Йорк) = 1;
cbey (1: NX, 1: Нью-Йорк) = 1/e0;
ccey (1: NX, 1: Нью-Йорк) = 1/e0;
cabz (1: NX, 1: Нью-Йорк) = 1;
cbbz (1: NX, 1: Нью-Йорк) = DT;
cahz (1: NX, 1: Нью-Йорк) = 1;
cbhz (1: NX, 1: Нью-Йорк) = 1/m0;
cchz (1: NX, 1: Нью-Йорк) = 1/m0;
sigma_x (1: NX, 1: Нью-Йорк) = 0;
sigma_y (1: NX, 1: Нью-Йорк) = 0;

% исходного
для N = 1: N
если N * DT <1E-9
H (N) = (1 / 320) * (10-15 * COS (2 * PI * 1E9 * DT * N) 6 * COS (4 * пи * 1E9 * Л * N) ...
-COS (6 * PI * 1E9 * Л * N));
конец
конец
% caculation для coefficints в ПМЛ
M = 2;% вариации (заказ)
R = 1e-5;% reflecion коэффициента

SM =- (SQRT (e0/m0 ))*(( M 1) * LOG (R) / (2 * D * DX));

для I = 1: D
Sigma (I) = Sm * ((I-1) / D) ^ M;
сигм (I) = Sm * ((I-0.5) / D) ^ M;
конец
%
% УИК за слоем ПМЛ
EX (1: NX, NY) = 0;% УИК
EX (1: NX, 1) = 0;% УИК
Е.Я. (1,1: Нью-Йорк) = 0;% УИК
Е.Я. (NX, 1: Нью-Йорк) = 0;% УИКза II = 1: NX
для JJ = 1: Нью-Йорк

sigma_x (II, JJ) = 0;
sigma_y (II, JJ) = 0;
sigmas_x (II, JJ) = 0;
sigmas_y (II, JJ) = 0;
если (II <= D)
sigma_x (II, JJ) = Sigma (D-II 1);
sigmas_x (II, JJ) = сигм (D-II 1);
конец
если (II> NX-D)
sigma_x (II, JJ) = Sigma (II-NX D);
sigmas_x (II, JJ) = сигм (II-NX D);
конецесли (JJ <= D)
sigma_y (II, JJ) = Sigma (D 1-й);
sigmas_y (II, JJ) = сигм (D 1-й);
конец
если (JJ> NY-D)
sigma_y (II, JJ) = Sigma (JJ-Йорк D);
sigmas_y (II, JJ) = сигм (JJ-NY D);
конец
cbex (II, JJ) = (2 * E0 sigma_x (II, JJ) * DT) / (2 * E0 * E0);
ccex (II, JJ) = (2 * E0-sigma_x (II, JJ) * DT) / (2 * E0 * E0);
caey (II, JJ) = (2 * E0-sigma_x (II, JJ) * DT) / (2 * E0 sigma_x (II, JJ) * DT);
cbey (II, JJ) = (2 * E0 sigma_y (II, JJ) * DT) / (2 * E0 * E0 sigma_x (II, JJ) * E0 * DT);
ccey (II, JJ) = (2 * E0-sigma_y (II, JJ) * DT) / (2 * E0 * E0 sigma_x (II, JJ) * E0 * DT);
cadx (II, JJ) = (2 * E0-sigma_y (II, JJ) * DT) / (2 * E0 sigma_y (II, JJ) * DT);
cbdx (II, JJ) = 2 * E0 * DT / (2 * E0 sigma_y (II, JJ) * DT);
cahz (II, JJ) = (2 * E0-sigmas_y (II, JJ) * DT) / (2 * E0 sigma_y (II, JJ) * DT);
cbhz (II, JJ) = (2 * E0) / (2 * E0 sigmas_y (II, JJ) * DT) / M0;
cchz (II, JJ) = (2 * E0) / (2 * E0 sigmas_y (II, JJ) * DT) / M0;
cabz (II, JJ) = (2 * E0-sigmas_x (II, JJ) * DT) / (2 * E0 sigmas_x (II, JJ) * DT);
cbbz (II, JJ) = (2 * E0 * DT) / (2 * E0 sigmas_x (II, JJ) * DT);

конец
конецдля N = 1: N

DX1 = DX;
DY1 = Dy;
за II = 1: NX
для JJ = 1: Нью-Йорк
если (JJ> 1)
DX (II, JJ) = cadx (II, JJ) * DX (II, JJ) cbdx (II, JJ) * (HZ (II, JJ)-HZ (II, JJ-1)) / Dy;
еще
DX (II, JJ) = 0;
конец
конец
конецза II = 1: NX
для JJ = 1: Нью-Йорк
если (II> 1)
DY (II, JJ) = Кади (II, JJ) * DY (II, JJ)-cbdy (II, JJ) * (HZ (II, JJ)-HZ (II-1, JJ)) / DX;
еще
DY (II, JJ) = 0;
конец

конец
конец

за II = 1: NX
для JJ = 1: Нью-Йорк
EX (II, JJ) = caex (II, JJ) * EX (II, JJ) cbex (II, JJ) * DX (II, JJ)-ccex (II, JJ) * DX1 (II, JJ);

конец
конец
за II = 1: NX
для JJ = 1: Нью-Йорк
Е.Я. (II, JJ) = caey (II, JJ) * Е.Ю. (II, JJ) cbey (II, JJ) * DY (II, JJ)-ccey (II, JJ) * DY1 (II, JJ);
конец
конецBZ1 = BZ;
за II = 1: NX-1
для JJ = 1: НГ-1
BZ (II, JJ) = cabz (II, JJ) * BZ (II, JJ)-cbbz (II, JJ) * ((EY (II 1, JJ)-EY (II, JJ)) / DX-( EX (II, JJ 1)-EX (II, JJ)) / ду);
конец
конец

за II = 1: NX-1
для JJ = 1: НГ-1
HZ (II, JJ) = cahz (II, JJ) * Гц (II, JJ) cbhz (II, JJ) * BZ (II, JJ)-cchz (II, JJ) * BZ1 (II, JJ);
конец
конец
Гц (это, JS) = Гц (это, JS) H (N);
% W (N) = Гц (это, D 2);
(П) = Гц (это, D 2);

конец
N = 1: N;
участка (N, W (N))