tarea-CI7411/plot_fa.m at master · ppizarror/tarea-CI7411 · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
function [w,faw] = plot_fa(fa, wmin, wmax, famin, famax, plot_title, yaxis_label, logx, NPOINTS, plot_colorw, show_legend)
%PLOT_FA Grafica la funcion del factor de amplificacion.
%
%   plot = plot_fa(fa, wmin, wmax, famin, famax, 'Plot title');
%   plot = plot_fa(fa, wmin, wmax, famin, famax, 'Plot title');
%   plot = plot_fa(fa, wmin, wmax, famin, famax, 'Plot title', false);
%   plot = plot_fa(fa, wmin, wmax, famin, famax, 'Plot title', true, 'FA');
%   plot = plot_fa(fa, wmin, wmax, famin, famax, 'Plot title', true, 'FA', 1000);
%   plot = plot_fa(fa, wmin, wmax, famin, famax, 'Plot title', true, 'FA', 100, 'k');
%   plot = plot_fa(fa, wmin, wmax, famin, famax, 'Plot title', true, 'FA', 100, 'k, false);
%
% Parametros:
%   fa              Funcion factor de amplificacion
%   wmin            Frecuencia minima de evaluacion, eje x
%   wmax            Frecuencia maxima de evaluacion, eje x
%   famax           Valor maximo de FA en el eje y
%   plot_title      Titulo del grafico
%   yaxis_label     Label de la leyenda en y
%   logx            Plot en semilogx
%   NPOINTS         Numero de puntos de evaluacion, por defecto es 100
%   plot_colorw     Color del grafico FA(w)
%   show_legend     Muestra la leyenda

%% Inicia variables
if ~exist('NPOINTS', 'var')
    NPOINTS = 100;
end
if ~exist('plot_colorw', 'var')
    plot_colorw = 'k';
end
if ~exist('show_legend', 'var')
    show_legend = false;
end
if ~exist('yaxis_label', 'var')
    yaxis_label = 'FA($\omega$)';
end
if ~exist('logx', 'var')
    logx = false;
end

%% Genera el grafico
plt = figure();

% Modifica la figura
movegui(plt, 'center');
set(gcf, 'name', plot_title);

%% Crea las variables
w = linspace(0.01, wmax, NPOINTS);
faw = zeros(NPOINTS, 1);
for i = 1:NPOINTS
    faw(i) = fa(w(i));
end

%% Grafica
if ~logx
    plot(w, faw, plot_colorw);
else
    semilogx(w, faw, plot_colorw);
end
hold on;
grid on;
title(plot_title);
xlabel('$\omega$', 'interpreter', 'latex');
ylabel(yaxis_label, 'interpreter', 'latex');
ylim([famin, famax]);
xlim([wmin, wmax]);
if show_legend
    legend({yaxis_label}, 'Location', 'northwest')
end
hold off;

end