Για ποιες διαστάσεις του h σου βγαίνει αυτός ο χρόνος εκτέλεσης; Η πολυπλοκότητα σου είναι O(N^4) ή O(N^2);

Tην fft2 θα χρησιμοποιήσεις. Υπάρχει και κάπου παραπάνω νομίζω πως το κανεις.

edit:

Κι εγώ έχω θέμα με τα διπλά γραφήματα. Μου βγαίνει ο μέσος χρόνος πιο ψηλά από τη θεωρητική πολυπλοκότητα, ακόμη και στο πρώτο με το fft2. Το έκανα με plotyy και όταν όρισα τον κάθε άξονα στο max του έχω αυτό το αποτέλεσμα της εικόνας.

Την ίδια εικόνα έχω ακόμη κι οταν ορίσω κατάλληλο συντελεστή A για τη θεωρητική πολυπλοκότητα.

Η πολυπλοκότητα δε θα είναι Ο(My∙Ny∙log(My∙Ny )) και αφού My=Mx+Mh -1 και  Νy=Νx+Νh-1 και h σταθερή θα γίνει Ο(Mx∙Nx∙log(My∙Ny )) ;
Μες στο λογαριθμο τι γίνεται; Θεωρώ κι εκεί log(Mx∙Nx); και μετά αυτό που λες Mx=Nx=N άρα γίνεται Ο((Ν^2)logN);

Ναι το y απ το data.m είναι το h. Άρα θα έχεις ως όρισμα στις συναρτήσεις σου τα x,y (για να υπολογίσεις τη συνέλιξη των x*y)

παιδιά πόση ώρα κάνει η myconv2 σας να τρέξει και με τι μέγεθος εικόνας h; Εγώ που την έκανα με 4 for-loops κάνει γυρω στις 3 ωρες να τρέξει η κάθε επανάληψη για x[128x128] και h[128x128]..   

H αλήθεια είναι πως με τόσο μικρή εικόνα η conv2 είναι αρκετά πιο γρήγορη από τη myconv2freq, οπότε θα προτιμήσω μεγαλύτερη 
Απλά δεν εχω καταλάβει αν πρέπει να χρησιμοποιήσουμε ένα μεγεθος για όλα τα Ν. Δηλαδή για μια εικόνα [4x4] είναι οκ να χρησιμοποιήσουμε εικόνα h διαστάσεων [256x256];

Δηλαδή θα μπορούσαμε να θεωρήσουμε εικόνα h διαστάσεων όσο και η μικρότερη εικόνα x [4x4] για k=2 ή ακόμη και [2x2] για παράδειγμα;

Στο ζητούμενο 3, η εικόνα x θα είναι ένας τυχαίος πίνακας διαστάσεων [NxN] όπου Ν=2^k, k=2,3,...,10
Για την εικόνα h που ζητάει να είναι σταθερών διαστάσεων τι θεωρούμε; Θα είναι κι αυτή [NxN] σαν τη x;

Σε ευχαριστώ πολύ, βοήθησε αρκετά να καταλάβω τι παίζει.
Ο τρόπος υλοποίησης που περιγράφει βασίζεται στο κέντρο του πίνακα h. Αν όμως ο πίνακας h δεν έχει κέντρο (πχ. είναι 2x2 ή 3x4) με βάση ποιο στοιχείο του πίνακα κινείσαι;

Έχω κολλήσει με τη myconv2. Στην ουσία είναι η υλοποίηση της εξίσωσης της συνέλιξης, που μας δίνει στην εκφώνηση, στο matlab, σωστά;
Επειδή δε το χω ακόμη με τα vectors, το έχω κάνει με for. Μου βγάζει όμως ότι κάποια στοιχεία είναι εκτός των ορίων των πινάκων που χρησιμοποιώ.
Οι πίνακες των πεπερασμένων xe και he έχουν διαστάσεις [My x Ny]

Θα γράψω ένα παράδειγμα μήπως δεν έχω καταλάβει κάτι σωστά στην υλοποίηση της εξίσωσης της συνέλιξης.
Έστω οι πίνακες των σημάτων εισόδου με διαστάσεις x [2x3] και h[2x2] άρα οι xe, he, y θα έχουν διαστάσεις [3,4].
Έστω ότι θέλω το στοιχείο y(1,1), θα έχω:

y(1,1)= xe(1,1) * he(0,0) + xe(1,2) * ye(0,-1) + xe(1,3) * he(0,-2) +  xe(2,1) * he(-1,0) + xe(2,2) * he(-1,-1) + xe(2,3) * he(-1,-2)

Tα στοιχεία  he(0,0) , he(0,-1) , he(0,-2),..  είναι εκτός των ορίων του πίνακα he και γι αυτό μου βγάζει το αντίστοιχο error. Αν μηδενίσω αυτά τα στοιχεία τότε το y(1,1)=0 που είναι λάθος. Τι δε βλέπω;