


BML_LOCAVGREF_MATRIX creates a local average referencing matrix Use as U = bml_locavgref_matrix(n,alpha)


0001 function U=bml_locavgref_matrix(n,alpha) 0002 0003 % BML_LOCAVGREF_MATRIX creates a local average referencing matrix 0004 % 0005 % Use as 0006 % U = bml_locavgref_matrix(n,alpha) 0007 % 0008 0009 if ~exist('alpha','var') 0010 alpha=1; 0011 end 0012 0013 U = eye(n); 0014 U(1,2)=-alpha; 0015 U(n,n-1)=-alpha; 0016 for i=2:(n-1) 0017 U(i,i+1)=-alpha/2; 0018 U(i,i-1)=-alpha/2; 0019 end 0020