毕业设计&课设-频域系统辨识Matlab工具箱.zip

function [cx,crx,crxmax,cry,crymax]=msinclip(varargin) %MSINCLIP - Multisine design with mimimum crest factor. % % Algorithm: time domain <--> frequency domain swapping, clipping in time % domain. The output is an object (structure) containing the complex Fourier % amplitudes; use MSINPREP to generate a tiime domain signal from these. % % Output arguments: % cx = fiddata object of the designed multisine. In the property 'input' % it contains the coefficients of the complex Fourier series: % the absolute values are equal to halves of the real amplitudes). % For the ZOH case (see below for crestmode), precompensation % with the inverse of the ZOH transfer function is applied, and the % calculated coefficients are those of the ZOH-generated excitation % signal. % crinfo = information about the crest factors. Fields of this structure: % crx = crest factor of the generated multisine, calculated with the % given oversampling factor. % crxmax = worst case crest factor of the multisine % cry = crest factor of the multisine at the output of the linear % system, calculated with the given oversampling factor. % crymax = worst case output crest factor of the multisine % % Input arguments: % Fdat = fiddata object (see 'help fiddata' and 'help(fiddata)', % with the basic quantities (in stand-alone applications % Fdat may be a structure with these fields): % freqpoints = vector of frequencies where the nonzero amplitudes are % given: the elements must be integer multiples of a df value, minimum % number of sines: 2. The frequency vector must be strictly % monotonously increasing. % input = absolute values of the desired nonzero complex amplitudes % (coefficients of the complex Fourier series at the corresponding % frequencies: halves of the real coefficients). % If any element is complex, the phases of the input vector will be % used as starting values, otherwise the Schroeder multisine is used. % For elements given with value NaN, auxiliary sine amplitudes are % returned ("snowing"): these are not included into the % calculation of the effective value of the useful excitation % signal, only help to decrease the peak value. % Default for these amplitudes: ones(length(fv),1) % Initial values for the snowing case can be given by a special % field of runmod: 'snowinginput' is the NaN-containing vector, the % real values are neglected. % output = (optional) outputs: input amplitudes multiplied by the complex % transfer function values at the given frequencies. % If the output is given, input-output optimization will be performed. % runmod = structure of run modifiers. Usually not necessary. % For the details, execute 'msinclip runmod'. % % Usage: [cx,crinfo]=msinclip(Fdat,runmod); % Examples: % [cx,crestx]=msinclip(fiddata([],[],[1:15]'/256)); % %Schroeder starting phases: % cx=msinclip(fiddata([],ones(12,1),4:15),struct('initset','Schroeder')); % %Prepare time series, 2048 points, fs = 1024 Hz: % timeobject=msinprep(cx,2048,1024); % %ZOH design: % cx=msinclip(fiddata([],ones(21,1)/sqrt(42),[0.2:0.01:0.4]*1e3),... % struct('crestmode','ZOH','N',100,'itno',350)); % Copyright (c) I. Kollar and VUBrussel, ELEC, 1991-2004 % % See also: DIBS, MSINPREP, CRESTMIN. if nargin==0, error('No input argument'), end if (nargin==1)&isstr(varargin{1}) if strncmpi(varargin{1},'preload',4) return %preload only elseif strcmp(varargin{1},'oldhelp') oldhelp msinclip return elseif strcmp(varargin{1},'runmod')|strcmp(varargin{1},'runpar') msincrmt return elseif strcmp(varargin{1},'secret') fprintf([' ''Secret'' way of setting the number of iterations:\n',... ' global msinclip_max_iteration_number\n msinclip_max_iteration_number=<iterations>;\n']) return end error(['First argument ''',varargin{1},''' is not allowed']) end % global msinclip_max_iteration_number c=computer; MatlV=version; isPC=strcmp(c(1:2),'PC'); isPlot=0; plotdone=0; if isPC %Matlab5 under Win95 still greys out calldrawnow=1; calliterctrl=1; else %regular case, e.g. Matlab5 on any platform, still may grey out calldrawnow=1; calliterctrl=1; end %calldrawnow=0; calliterctrl=0; %Force to eliminate drawnows %calldrawnow=1; calliterctrl=1; %Force to call drawnows % itctrllastchecked='Continue'; termcond=0; if calldrawnow, drawnow, end % crdef=1; %crdef=1: cr=max(abs(multisine))/effval, crdef=2:multisinepp/2/effval %mbint(crdef) %assertion (neglected by Matlab 5.3, but not by 5.2) yes=1; no=0; global halffiginmsinclip if ~exist('halffiginmsinclip'), halffiginmsinclip=''; end if ~isempty(halffiginmsinclip) halffig=halffiginmsinclip; halffiginmsinclip=''; else halffig='n'; end %for halffig='y', the plot of input minimization is subplot(1,2,1) cldef='JS'; %clipping level; JS: Johan Schoukens, PG: Patrick Guillaume, % VDO: Edwin van der Ouderaa % Fdat=varargin{1}; OK=1; if isstruct(Fdat) if ~isfield(Fdat,'freqpoints'), OK=0; fprintf(['Fdat field ''freqpoints'' is missing']) elseif ~isfield(Fdat,'input'), OK=0; fprintf(['Fdat field ''input'' is missing']) elseif ~isfield(Fdat,'output'), Fdat.output=''; end end if OK==0, error('First input argument is an improper structure'), end initset=''; showmsgs=0; Nmax=inf; if exist('fdguidev.mat'), showmsgs=1; end %messages for development if isa(Fdat,'fiddata')|isstruct(Fdat) if nargout>2, error('Too many output arguments'), end error(nargchk(1,2,nargin)) initset='random'; %default for the new call fv=Fdat.freqpoints; ampv=Fdat.input; tf=Fdat.output; if ~isempty(tf) ind=find(~isnan(Fdat.input)); tf(ind)=tf(ind)./Fdat.input(ind); end gmod=''; itno=[]; ovs=[]; N=[]; cl0=[]; sr=[]; pld=[]; crestmode='BL'; if (nargin>=2)&~isempty(varargin{2})&isstruct(varargin{2}) runmod=varargin{2}; fn=fieldnames(runmod); for ii=1:length(fn) if strncmp(fn{ii},'initset',4) initset=getfield(runmod,fn{ii}); if strncmpi(initset,'Schroeder',3) elseif strncmpi(initset,'random',4) elseif strncmpi(initset,'input',2) elseif isempty(initset) else error(['initset ''',initset,''' is invalid']) end elseif strncmp(fn{ii},'graphs',5), gmod=getfield(runmod,fn{ii}); elseif strncmp(fn{ii},'showmsgs',5) showmsgs=getfield(runmod,fn{ii}); showmsgs=~strncmp(showmsgs,'n',1); elseif strncmp(fn{ii},'crestmode',6) crestmode=getfield(runmod,fn{ii}); if strcmpi(crestmode,'discrete'), gmod=[gmod,'ZOHd']; elseif strcmpi(crestmode,'ZOH')|strcmpi(crestmode,'DAC'), gmod=[gmod,'ZOHc']; elseif strncmpi(crestmode,'filtered-DAC',8)|strcmpi(crestmode,'fDAC') gmod=[gmod,'fDAC']; elseif strcmpi(crestmode,'BL') %default else error(['crestmode is not allowed: ''',crestmode,'''']) end elseif strncmp(fn{ii},'itno',2), itno=getfield(runmod,fn{ii}); elseif strncmp(fn{ii},'pmin',4) %for crestmin only warning('pmin has no meaning in msinclip') elseif strncmp(fn{ii},'relvarmax',4) %for crestmin only warning('relvarmax has no meaning in msinclip') elseif strncmp(fn{ii},'pmax',4) %for crestmin only warning('pmax has no meaning in msinclip') elseif strncmp(fn{ii},'oversampling',5), ovs=getfield(runmod,fn{ii}); elseif strcmpi(fn{ii},'N'), N=getfield(runmod,fn{ii}); elseif strcmpi(fn{ii},'Nmax'), Nmax=getfield(runmod,fn{ii}); elseif strncmp(fn{ii},'cliplevel',5), cl0=getfield(runmod,fn{ii}); elseif strncmp(fn{ii},'slewrate