spm_dcm_fmri_csd_gen.m
8a186cd82567master

spm_dcm_fmri_csd_gen.m
View Options

			function [y] = spm_dcm_fmri_csd_gen(x,v,P)
			% Conversion routine for DEM inversion of DCM for CSD (fMRI)
			% FORMAT [y] = spm_dcm_fmri_csd_gen(x,v,P)
			%
			% This routine computes the spectral response of a network of regions
			% driven by endogenous fluctuations and exogenous (experimental) inputs.
			% It returns the complex cross spectra of regional responses as a
			% three-dimensional array. The endogenous innovations or fluctuations are
			% parameterised in terms of a (scale free) power law, in frequency space.
			%
			%
			%__________________________________________________________________________
			% Copyright (C) 2008 Wellcome Trust Centre for Neuroimaging

			% Karl Friston
			% $Id: spm_dcm_fmri_csd_gen.m 5736 2013-11-10 13:17:10Z karl $


			% global DCM and evaluate orginal generative model
			%==========================================================================
			global GLOBAL_DCM

			% conditional prediction
			%--------------------------------------------------------------------------
			y = feval(GLOBAL_DCM.M.IS,v,GLOBAL_DCM.M,GLOBAL_DCM.U);
			y = feval(GLOBAL_DCM.M.FS,y,GLOBAL_DCM.M);
			y = spm_vec(y);