matlab波形的thd波形,Matlab图片处理函数小记

function ha = tight_subplot (Nh, Nw, gap, marg_h, marg_w)

%TIGHT_SUBPLOT suplot子图间距调整，自适应调整间距、空白

% 代码参考: http://blog.csdn.net/lanchunhui/article/details/49820721

%

% 输入：行数Nh,列数Nw 用法同subplot(row,col)

% gap 归一化坐标轴间间距 (0...1)

% or [gap_h gap_w] 垂直、水平 for different gaps in height and width

% marg_h 归一化上下边距 (0...1)

% or [lower upper] for different lower and upper margins

% marg_w 归一化左右边距 (0...1)

% or [left right] for different left and right margins

% 输出：ha 子图的句柄数组

% starting from upper left corner, going row-wise as in

% going row-wise as in

%

% Example: ha = tight_subplot(3,2,[.01 .03],[.1 .01],[.01 .01])

% for ii = 1:6; axes(ha(ii)); plot(randn(10,ii)); end

% set(ha(1:4),'XTickLabel',''); set(ha,'YTickLabel','')

if nargin<3; gap = .02; end

if nargin<4 || isempty(marg_h); marg_h = .05; end

if nargin<5; marg_w = .05; end

if numel(gap)==1

gap = [gap gap];

end

if numel(marg_w)==1

marg_w = [marg_w marg_w];

end

if numel(marg_h)==1

marg_h = [marg_h marg_h];

end

axh = (1-sum(marg_h)-(Nh-1)*gap(1))/Nh;

axw = (1-sum(marg_w)-(Nw-1)*gap(2))/Nw;

py = 1-marg_h(2)-axh;

k=0;

ha = zeros(Nh*Nw,1);

for ih = 1:Nh

px = marg_w(1);

for ix = 1:Nw

k = k+1;

ha(k) = axes('Units','normalized','Position',[px py axw axh], ...

'XTickLabel','', 'YTickLabel','');

px = px+axw+gap(2);

end

py = py-axh-gap(1);

end

由fig仿真/实验图分析电流THD

%BYN_FFTBOX Compute Current THD,