--- a
+++ b/code/preprocessing/lib/my_PSD.m
@@ -0,0 +1,46 @@
+function [Fxx,Pxx]=my_PSD(data,NFFT,Fs,win,overlap,minMaxHz,flag)
+
+
+over=round(NFFT-NFFT*(overlap/100));
+
+
+resolution=Fs/NFFT;
+
+minMaxHz(1)=round(minMaxHz(1)*(1/resolution)+1);
+minMaxHz(2)=round(minMaxHz(2)*(1/resolution)+1);
+
+%commento fine --------------------------------
+
+if(overlap==0)
+    if ((nargin>=6)&(nargin<=7))
+        [r,c]=size(data);
+        
+        j=1;
+        for(i=1:r)
+           % subplot(r/nCol,nCol,i);
+            [Pxx(:,i),Fxx(:,i)]=psd(data(j,:),NFFT,Fs,window(win,NFFT),'NOVERLAP');
+            if strcmpi(flag,'on')
+            plot(Fxx(minMaxHz(1):minMaxHz(2),:),sqrt(Pxx(minMaxHz(1):minMaxHz(2),:)));
+            end
+            j=j+1;
+            %pause;
+        end
+    end
+end
+
+if(overlap>0)
+    if ((nargin>=6)&(nargin<=7))
+        [r,c]=size(data);
+        
+        j=1;
+        for(i=1:r)
+            %subplot(r/nCol,nCol,i);
+            [Pxx(:,i),Fxx(:,i)]=psd(data(j,:),NFFT,Fs,window(win,NFFT),over);
+            if strcmpi(flag,'on')
+            plot(Fxx(minMaxHz(1):minMaxHz(2),:),sqrt(Pxx(minMaxHz(1):minMaxHz(2),:)));
+            end
+            j=j+1;
+            %pause;
+        end
+    end
+end
\ No newline at end of file