matlab底层代码可视化高斯分布
为了回顾忘的差不多的matlab，使用底层代码实现一些界面功能。随便可视化下一元和多元高斯分布。(加深机器学习中 吴恩达 - 15.7的多元高斯分布 的印象 )
为什么不用GUIDE实现界面？虽然GUIDE做界面简单、好用，但即便界面多么简单、其生成的程序都会有两个文件，一个界面文件fig一个逻辑文件m，界面与逻辑分离在这些简单的小程序里面优势并不大，所以我还是倾向使用底层代码，它可以只使用一个m文件完成界面与逻辑。
1. 可视化一元高斯分布：
中文wiki：一元高斯分布
MATLAB 底层代码可视化一元高斯分布：
function [] = UnivarGaussianDistribution()
% Plot Univariate Gaussian distribution
% Parameter
[x_min, x_max] = deal(-2, 10);
[sigma_min, sigma_max] = deal(-2, 10);
[mu_min, mu_max] = deal(-2, 10);
% Create data
x = linspace(x_min, x_max, 1000);
global sigma mu
% Default sigma, mu is average of maximum and minimum.
sigma = (sigma_min + sigma_max) / 2;
mu = (mu_min + mu_max) / 2;
% Now create the other GUI
% position: [x, y, width, height]
S.fh = figure('units','pixels',...
'position',[500 300 600 500],...
'menubar','none',...
'name','GaussianDistributionVisulization-Univariate ',...
'numbertitle','off',...
'resize','off');
% axes
ha = axes('Parent', S.fh, ...
'Units', 'pixels', ...
'Position', [50 120 500 320]);
% plot data
line = plot(x, gaussianfunc(x, mu, sigma), 'r', 'linewidth', 1);
xlim([x_min, x_max]); ylim([-0.5, 2.5]);
title('y = f(x|\mu, \sigma) = \frac{1}{\sigma \sqrt{2\pi}}e^\frac{-(x - \mu)^2}{2 \sigma^2}',...
'Interpreter','latex', 'FontSize', 18)
% slider_1
S.tx(1) = uicontrol('style','text',...
'unit','pixel',...
'position',[150 10 40 20],...
'string','mu',...
'fontsize',10);
S.sl(1) = uicontrol('style','slide',...
'unit','pixel',...
'position',[200 10 200 20],...
'min', mu_min, 'value', mu,'max', mu_max);
S.ed(1) = uicontrol('style','edit',...
'unit','pixel',...
'position',[410 10 40 20],...
'fontsize',10,...
'string',num2str(mu)); % Displays the value.
% slider_2
S.tx(2) = uicontrol('style','text',...
'unit','pixel',...
'position',[150 50 40 20],...
'string','sigma',...
'fontsize',10);
S.sl(2) = uicontrol('style','slide',...
'unit','pixel',...
'position',[200 50 200 20],...
'min', sigma_min, 'value', sigma, 'max', sigma_max);
S.ed(2) = uicontrol('style','edit',...
'unit','pixel',...
'position',[410 50 40 20],...
'fontsize',10,...
'string',num2str(sigma)); % Displays the value.
set([S.sl(:);S.ed(:)],'call',{@h_call, S, line, x});
% shared Callback.
end
function gaus_y = gaussianfunc(x, mu, sigma)
gaus_y = 1/(sqrt(2*pi) * sigma) * exp(-(x - mu).^2 / (2 * sigma.^2));
end
function [] = h_call(varargin)
% Callback for the slider.
[h, ~, S, line, x] = deal(varargin{:}); % Get the calling handle and structure.
global sigma mu
switch h % Who called?
case S.sl(1)
mu = get(S.sl(1),'val');
set(S.ed(1),'string', mu) % Set edit to current slider.
case S.sl(2)
sigma = get(S.sl(2),'val');
set(S.ed(2),'string', sigma) % Set edit to current slider.
case S.ed(1)
E = str2double(get(S.ed(1),'string')); % Numerical edit string.
SL1 = get(S.sl(1),{'min','value','max'}); % Get the slider's info.
if E >= SL1{1} && E <= SL1{3}
mu = E;
set(S.sl(1),'value',E) % E falls within range of slider.
else
set(S.ed(1),'string',SL1{2}) % User tried to set slider out of range.
end
case S.ed(2)
E = str2double(get(S.ed(2),'string')); % Numerical edit string.
SL2 = get(S.sl(2),{'min','value','max'}); % Get the slider's info.
if E >= SL2{1} && E <= SL2{3}
sigma = E;
set(S.sl(2),'value',E) % E falls within range of slider.
else
set(S.ed(2),'string',SL2{2}) % User tried to set slider out of range.
end
end
gaus_y = gaussianfunc(x, mu, sigma);
set(line, 'ydata', gaus_y);
end
运行样图：
滑动 slider 减小 sigma
增大 mu
2. 可视化多元高斯分布
中文wiki：多元高斯分布
function [] = MultiGaussianDistribution()
% Plot Multivariate Gaussian distribution
close
% Create data
[xy_min, xy_max] = deal(-2, 10);
var_range = 5; % mu, Sigma range +/-5
global Sigma mu line tog_value z_shaped
mu = [1, 3];
Sigma = [1, 0; 0, 4];
% Sigma must be a square, symmetric, positive definite matrix.
tog_value = 1;
[mu_min, mu_max] = deal(mu - var_range, mu + var_range);
[Sigma_min, Sigma_max] = deal(Sigma - var_range, Sigma + var_range);
[x, y] = meshgrid(linspace(xy_min, xy_max, 100),linspace(xy_min, xy_max, 100));
X = [x(:) y(:)];
z = mvnpdf(X, mu, Sigma);
z_shaped = reshape(z, 100, 100);
% Now create the other GUI
% position: [x, y, width, height]
fg = figure('units','pixels',...
'position',[500 200 665 500],...
'name','GaussianDistributionVisulization-Multivariate ',...
'numbertitle','off',...
'resize','off');
% axes
S.ax = axes('Parent', fg, ...
'Units', 'pixels', ...
'Position', [100 120 500 320]);
% plot data
line = surf(x, y, z_shaped,...
'edgecolor', 'interp');
xlim([xy_min, xy_max]); ylim([xy_min, xy_max]);
latex_text = ['{f_x(x_1,\ldots ,x_k)={\frac {1}{\sqrt {(2\pi)^k|{',...
'\Sigma}|}}} \exp \left(-{\frac {1}{2}}({{X} }-{{\mu }})^{{T}',...
' }{{\Sigma }}^{-1}({{X} }-{{\mu }})\right)}'];
title(latex_text,...
'Interpreter', 'latex', 'FontSize', 12);
xlim([xy_min, xy_max]);ylim([xy_min, xy_max]);zlim([-0.02, 0.09]);
xlabel('x'); ylabel('y');
%% check box for surf or contourf
S.togbtn = uicontrol('style', 'togglebutton',...
'unit', 'pixel',...
'position', [20, 450, 50, 30],...
'string', 'surf',...
'val', tog_value,...
'fontsize', 10,...
'TooltipString', 'Switch figure style between 3-D surface and contour.');
%% slider 1-2: mu(1), mu(2)
position_mu = {[5 50 75 20], [80 50 140 20]; [5 10 75 20], [80 10 140 20]};
for ii = 1:2
S.tx(ii) = uicontrol('style', 'text',...
'unit', 'pixel',...
'position', position_mu{ii, 1},...
'string', sprintf('mu(%d): %d', ii, mu(ii)),...
'fontsize', 10);
S.sl(ii) = uicontrol('style', 'slide',...
'unit', 'pixel',...
'position', position_mu{ii, 2},...
'min', mu_min(ii), 'value', mu(ii),'max', mu_max(ii));
end
%% slider 3-6: Sigma(1, 1)-Sigma(2, 2)
position_Sigma = {
[225 40 75 40], [300 50 140 20];
[225 0 75 40], [300 10 140 20];
[445 40 75 40], [520 50 140 20];
[445 0 75 40], [520 10 140 20];
};
for ii = 1:4
[row, column] = ind2sub(size(Sigma), ii);
S.tx(ii + 2) = uicontrol('style', 'text',...
'unit', 'pixel',...
'position', position_Sigma{ii, 1},...
'string', sprintf('Sigma(%d, %d): %d', row, column, Sigma(ii)),...
'fontsize', 10);
S.sl(ii + 2) = uicontrol('style', 'slide',...
'unit', 'pixel',...
'position', position_Sigma{ii, 2},...
'min', Sigma_min(ii), 'value', Sigma(ii),'max', Sigma_max(ii));
end
set([S.togbtn; S.sl(:)],'call',{@h_call, S, X, x, y}); % shared Callback.
end
function [] = h_call(varargin)
% Callback for the slider.
[h, ~, S, X, x, y] = deal(varargin{:}); % Get the calling handle and structure.
global Sigma mu line tog_value z_shaped
switch h % Who called?
case S.togbtn
tog_value = get(S.togbtn, 'value'); % Get toggle button's state
if tog_value
line = surf(x, y, reshape(z_shaped, 100, 100),...
'edgecolor', 'interp');
set(S.togbtn, 'string', 'surf');
else
line = contour(x, y, reshape(z_shaped, 100, 100),...
'ShowText', 'on', 'LineWidth', 1);
set(S.togbtn, 'string', 'contour');
end
case S.sl(1)
mu(1) = get(S.sl(1),'val'); % Get slider's value.
set(S.tx(1), 'string',['mu(1): ', num2str(mu(1))]); % update text box.
case S.sl(2)
mu(2) = get(S.sl(2),'val'); % Get slider's value.
set(S.tx(2), 'string',['mu(2): ', num2str(mu(2))]); % update text box.
case S.sl(3)
Sigma(1) = get(S.sl(3),'val'); % Get slider's value.
set(S.tx(3), 'string',['Sigma(1,1): ', num2str(Sigma(1))]); % update text box.
case S.sl(4)
Sigma(2) = get(S.sl(4),'val');
set(S.tx(4), 'string',['Sigma(2,1): ', num2str(Sigma(2))]); % update text box.
Sigma(3) = Sigma(2); % make sure Sigma is symmetric matrix.
set(S.sl(5), 'value', Sigma(2)); % Adjust the slider to match the value.
set(S.tx(5), 'string',['Sigma(1,2): ', num2str(Sigma(2))]);
case S.sl(5)
Sigma(3) = get(S.sl(5),'val');
set(S.tx(5), 'string',['Sigma(1,2): ', num2str(Sigma(3))]); % update text box.
Sigma(2) = Sigma(3); % Make sure Sigma is symmetric matrix.
set(S.sl(4), 'value', Sigma(3)); % Adjust the slider to match the value.
set(S.tx(4), 'string',['Sigma(2,1): ', num2str(Sigma(3))]); % update text box.
case S.sl(6)
Sigma(4) = get(S.sl(6),'val');
set(S.tx(6), 'string',['Sigma(2,2): ', num2str(Sigma(4))]);
end
try
z_shaped = reshape(mvnpdf(X, mu, Sigma), 100, 100); % if Sigma not square, symmetric, positive definite matrix.
if tog_value
%{
I want to improve script performance. Casue I think update
axes's zdata faster than plot a new axes.
%}
% line = surf(x, y, z_shaped,...
% 'edgecolor', 'interp');
% % line = surf(x, y, z_shaped,...
% % 'edgecolor', 'interp');
% set(S.togbtn, 'string', 'surf');
set(line, 'zdata', z_shaped)
else
h_contour = allchild(S.ax);
set(h_contour, 'zdata', z_shaped)
%{
I want to improve script performance. Casue I think update
axes's zdata faster than plot a new axes.
%}
% line = contour(x, y, z_shaped,... % I want to improve script performance.
% 'ShowText', 'on', 'LineWidth', 1);
% set(S.togbtn, 'string', 'contour');
end
catch ME
if (strcmp(ME.identifier, 'stats:mvnpdf:BadMatrixSigma'))
disp('Incorrect Sigma: ');disp(Sigma);
warndlg('Sigma must be a square, symmetric, positive definite matrix.')
else
errordlg(ME.message)
end
end
end
图例：
点击左上角按钮切换成等势线-contour
