function [tfd, t, f] = margenau_hill1(x, fs, nfreq);% margenau_hill1 -- Compute samples of the type I Margenau-Hill %                    distribution.%%  Usage%    [tfd, t, f] = margenau_hill1(x, fs, nfreq)%%  Inputs%    x     signal vector%    fs    sampling frequency of x (optional, default is 1 sample/second)%    nfreq number of samples to compute in frequency (optional, default%          is the length of x)%%  Outputs%    tfd  matrix containing the Margenau-Hill distribution of signal x.  %         If x has length N, then tfd will be nfreq by N. (optional)%    t    vector of sampling times (optional)%    f    vector of frequency values (optional)%% If no output arguments are specified, then the Margenau-Hill distribution % is displayed using ptfd(tfd, t, f).% Copyright (C) -- see DiscreteTFDs/Copyrighterror(nargchk(1, 3, nargin));if (nargin == 1)  [tfd, t, f] = rihaczek1(x);elseif (nargin == 2)  [tfd, t, f] = rihaczek1(x, fs);elseif (nargin == 3)  [tfd, t, f] = rihaczek1(x, fs, nfreq);endtfd = real(tfd);if (nargout == 0)  ptfd(tfd, t, f);  clear tfdend