Szerintetek helyes a következő kód?

% Vezeteknev_V.m

% 1. Adatok betöltése

data = load('Teljes Nev.txt'); % vagy readtable ha .csv-t használsz

wavelength = data(:, 1); % hullámhossz adatok

spectrum = data(:, 2); % spektrum adatok

% 2. A háttér (a + b*x) illesztése

background = @(params, x) params(1) + params(2)*x;

% A két Gauss függvény összegének illesztése

gauss_sum = @(params, x) params(3)*exp(-((x-params(4))/params(5)).^2) + ...

params(6)*exp(-((x-params(7))/params(8)).^2);

% Az illesztett modell teljes függvénye (háttér + Gauss)

model = @(params, x) background(params, x) + gauss_sum(params, x);

% 3. Paraméterek kezdő értékeinek meghatározása

% a, b, a1, b1, c1, a2, b2, c2

% Kezdeti értékek tippelése:

a_init = min(spectrum); % háttér érték (alapvonal)

b_init = 0; % háttér lineáris

a1_init = max(spectrum) - a_init; % első Gauss amplitúdó

b1_init = wavelength(spectrum == max(spectrum)); % első csúcs helye

c1_init = 2; % első Gauss szélessége (kb.)

a2_init = a1_init * 0.5; % második Gauss amplitúdó

b2_init = b1_init + 5; % második csúcs helye (kb. 5 nm különbség)

c2_init = c1_init; % második Gauss szélessége (kb. ugyanakkora)

% Kezdő paraméterek vektora

init_params = [a_init, b_init, a1_init, b1_init, c1_init, a2_init, b2_init, c2_init];

% 4. Nemlineáris illesztés

options = optimset('Display', 'off');

params_opt = lsqcurvefit(@(params, x) model(params, x), init_params, wavelength, spectrum, [], [], options);

% 5. Illesztett adatok és maradványok számítása

fitted_spectrum = model(params_opt, wavelength);

residuals = spectrum - fitted_spectrum;

% 6. Grafikonok ábrázolása

% 6.1 Mért és illesztett spektrum

figure;

subplot(2,1,1);

plot(wavelength, spectrum, 'b-', 'DisplayName', 'Mért spektrum');

hold on;

plot(wavelength, fitted_spectrum, 'r--', 'DisplayName', 'Illesztett spektrum');

xlabel('Hullámhossz (nm)');

ylabel('Fotonszám (counts)');

title('Mért és illesztett spektrum');

legend('show');

% 6.2 Maradványok grafikonja

subplot(2,1,2);

plot(wavelength, residuals, 'k-', 'DisplayName', 'Maradványok');

xlabel('Hullámhossz (nm)');

ylabel('Maradványok (counts)');

title('Maradvány függvény');

legend('show');

% 7. c1 és c2 paraméterek kiíratása

disp(['c1: ', num2str(params_opt(5)), ' nm']);

disp(['c2: ', num2str(params_opt(8)), ' nm']);