Hi, i want to get a gaussian curve fitting plot, function is y = Aexp-((x-mu)^2/2sig^2). y axis must be between -0.2 and 1 with counting by 0.2 per each step, any tips?

Question

Baris Unsal 2022-1-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1620840-hi-i-want-to-get-a-gaussian-curve-fitting-plot-function-is-y-a-exp-x-mu-2-2-sig-2-y-axis-mu

回答： Image Analyst 2022-1-2

f = @(x,A,mu,sig) A*exp(-(x-mu).^2)/(2*sig^2);
A = 1; mu = 7; sig = 1.5;
x = 0:2:10;
y = -0.2:0.2:;
y = f(x,a,b,c);
plot(x,y)

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Torsten 2022-1-2

编辑：Torsten 2022-1-2

在 MATLAB Online 中打开

What do you want to fit against what ?

f = @(x,A,mu,sig) A*exp(-(x-mu).^2)/(2*sig^2);
A = 1; mu = 7; sig = 1.5;
x = 0:2:10;
y = f(x,A,mu,sig);
plot(x,y)

Arrange the y-axis by using appropriate axes properties:

https://de.mathworks.com/help/matlab/ref/matlab.graphics.axis.axes-properties.html

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Matt J 2022-1-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1620840-hi-i-want-to-get-a-gaussian-curve-fitting-plot-function-is-y-a-exp-x-mu-2-2-sig-2-y-axis-mu#answer_866635

编辑：Matt J 2022-1-2

在 MATLAB Online 中打开

f = @(x,A,mu,sig) A*exp(-(x-mu).^2)./(2.*sig^2);

A = 1; mu = 7; sig = 1.5;

x = (0:2:10);

y = f(x,A,mu,sig);

plot(x,y);

ylim([-0.2,1])

yticks(-0.2:0.2:1)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Image Analyst 2022-1-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1620840-hi-i-want-to-get-a-gaussian-curve-fitting-plot-function-is-y-a-exp-x-mu-2-2-sig-2-y-axis-mu#answer_866655

在 MATLAB Online 中打开

See my well-commented code below (and attached demos) to fit a Guassian to any data. The code below uses your data points.

% Uses fitnlm() to fit a non-linear model (a single Gaussian curve) through noisy data.
% Requires the Statistics and Machine Learning Toolbox, which is where fitnlm() is contained.
% Initialization steps.
clc;    % Clear the command window.
close all;  % Close all figures (except those of imtool.)
clear;  % Erase all existing variables. Or clearvars if you want.
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 20;
% Create the X coordinates from 0 to 20 every 0.5 units.
A = 1; 
X = 0:2:10;
mu = 7; % Mean, center of Gaussian.
sigma = 1.5; % Standard deviation.
% Define function that the X values obey.
Y = A * exp(-(X - mu) .^ 2 / sigma); % Get a vector.  No noise in this Y yet.
% Now we have noisy training data that we can send to fitnlm().
% Plot the noisy initial data.
plot(X, Y, 'b*', 'LineWidth', 2, 'MarkerSize', 15);
grid on;
% Convert X and Y into a table, which is the form fitnlm() likes the input data to be in.
tbl = table(X(:), Y(:));
% Define the model as Y = a + exp(-b*x)
% Note how this "x" of modelfun is related to big X and big Y.
% x((:, 1) is actually X and x(:, 2) is actually Y - the first and second columns of the table.
modelfun = @(b,x) b(1) * exp(-(x(:, 1) - b(2)).^2/b(3));  
beta0 = [1,6,2]; % Guess values to start with.  Just make your best guess.
% Now the next line is where the actual model computation is done.
mdl = fitnlm(tbl, modelfun, beta0);
% Now the model creation is done and the coefficients have been determined.
% YAY!!!!
% Extract the coefficient values from the the model object.
% The actual coefficients are in the "Estimate" column of the "Coefficients" table that's part of the mode.
coefficients = mdl.Coefficients{:, 'Estimate'}
% Let's do a fit, but let's get more points on the fit, beyond just the widely spaced training points,
% so that we'll get a much smoother curve.
X = linspace(0, 20, 1920); % Let's use 1920 points, which will fit across an HDTV screen about one sample per pixel.
% Create smoothed/regressed data using the model:
yFitted = coefficients(1) * exp(-(X - coefficients(2)).^2 / coefficients(3));
% Now we're done and we can plot the smooth model as a red line going through the noisy blue markers.
hold on;
plot(X, yFitted, 'r-', 'LineWidth', 2);
grid on;
title('Exponential Regression with fitnlm()', 'FontSize', fontSize);
xlabel('X', 'FontSize', fontSize);
ylabel('Y', 'FontSize', fontSize);
legendHandle = legend('Noisy Y', 'Fitted Y', 'Location', 'northeast');
legendHandle.FontSize = 25;
% Set up figure properties:
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0 0 1 1]);
% Get rid of tool bar and pulldown menus that are along top of figure.
% set(gcf, 'Toolbar', 'none', 'Menu', 'none');
% Give a name to the title bar.
set(gcf, 'Name', 'Demo by ImageAnalyst', 'NumberTitle', 'Off') 