Помощь с научной работой

clc; 
clear; 
close all;
 
% Настройка параметров графиков
set(groot, 'defaultFigureColor', 'w');
 
% Параметры
f0 = 10e9;                % Несущая частота
c = 3e8;                  % Скорость света
lambda = c / f0;          % Длина волны             
N1 = 15; N2 = 35;         % Количество элементов
k1 = 0.34; k2 = 0.69;     % Коэффициенты для Хэмминга
d0 = lambda / 2;          % Базовое расстояние между элементами
 
% Логарифмическое распределение
sym_log_spacing = @(n) sign(n) * d0 * log(abs(n) + 1);
 
% Окно Хэмминга
hamming_spacing = @(n, N, k) sign(n) * k * (N-1) * d0 * (0.54 - 0.46 * cos(2 * pi * n / (N-1)));
 
% Индексы элементов
n1 = -(N1-1)/2:(N1-1)/2;
n2 = -(N2-1)/2:(N2-1)/2;
 
% Расчеты расстояний
sym_log1 = arrayfun(sym_log_spacing, n1);
ham1 = arrayfun(@(n) hamming_spacing(n, N1, k1), n1);
uniform1 = n1 * d0;
 
sym_log2 = arrayfun(sym_log_spacing, n2);
ham2 = arrayfun(@(n) hamming_spacing(n, N2, k2), n2);
uniform2 = n2 * d0;
 
spacing_ratio = @(x) abs(diff(x)) / lambda;
 
% Построение графиков
figure;
subplot(2,2,1);
plot(n1(1:end-1) + 1, spacing_ratio(sym_log1), 'k-o'); hold on;
plot(n1(1:end-1) + 1, spacing_ratio(uniform1), 'b-s'); hold on;
plot(n1(1:end-1) + 1, spacing_ratio(ham1), 'r-^'); 
grid on; legend('Uniform','Sym-log','Hamming');
xlabel('Element number'); ylabel('Spacing to wavelength ratio');
ylim([0 1])
xlim([-8 8])
 
subplot(2,2,2);
plot(n2(1:end-1) + 1, spacing_ratio(sym_log2), 'k-o'); hold on;
plot(n2(1:end-1) + 1, spacing_ratio(uniform2), 'b-s'); hold on;
plot(n2(1:end-1) + 1, spacing_ratio(ham2), 'r-^'); 
grid on; legend('Uniform','Sym-log','Hamming');
xlabel('Element number'); ylabel('Spacing to wavelength ratio');
ylim([0 1])
xlim([-18 18])
 
subplot(2,2,3);
plot(n1(1:end-1) + 1, spacing_ratio(sym_log1), 'k-o'); hold on;
plot(n1(1:end-1) + 1, spacing_ratio(uniform1), 'b-s'); hold on;
plot(n1(1:end-1) + 1, spacing_ratio(ham1 * 2), 'r-^'); 
grid on; legend('Uniform','Sym-log','Hamming');
xlabel('Element number'); ylabel('Spacing to wavelength ratio');
ylim([0 1])
xlim([-8 8])
 
subplot(2,2,4);
plot(n2(1:end-1) + 1, spacing_ratio(sym_log2), 'k-o'); hold on;
plot(n2(1:end-1) + 1, spacing_ratio(uniform2), 'b-s'); hold on;
plot(n2(1:end-1) + 1, spacing_ratio(ham2 * 2), 'r-^'); 
grid on; legend('Uniform','Sym-log','Hamming');
xlabel('Element number'); ylabel('Spacing to wavelength ratio');
ylim([0 1])
xlim([-18 18])

 hamming_spacing = @(n, N, k) sign(n) * k * (N-1) * d0 * (0.54 - 0.46 * cos(2 * pi * n / (N-1)));