Φασματική Ανάλυση και SALP tollbox

[seifer]

τώρα τίποτα για SALP Εργαλειοθήκη;
αυτός είναι ο κώδικας του σήματος:
[y, t] = sg_cissoid (β, στ, Δημ., <OSF>, <SNR>)
%
% Parameters: - βήτα: Vector με σύνθετες πλάτη
% - Στ: Vector με συχνότητα αξίες (-0,5 <f <0,5)
% - Δημ.: Αριθμός δειγμάτων
% - OSF: υπερδειγματοληψίας παράγοντα (προεπιλογή: OSF = 1)
% - SNR: σήματος προς θόρυβο σε [dB] (προεπιλογή: SNR = inf)
%
% Output: - y: Vector με σήμα δείγματα
% - Τ: Ώρα άξονα
και αυτός είναι ο κωδικός του Welch για παράδειγμα:

[PSD, f] = psd_welch (y, <Win>, <N_Overlap>, <N_FFT>, <fs>)
%
% Parameters: - y: Vector με σήμα δείγματα
% - Win: Το σήμα χωρίζεται σε K_seg αλληλεπικαλυπτόμενες
% Τμήματα μήκους ίσου με το μήκος της
% Από το παράθυρο και τα παράθυρα κάθε τμήμα
% (Προεπιλογή: 8 τμήματα με Hamming παράθυρο)
% - N_Overlap:
Ο αριθμός των δειγμάτων από επικάλυψη τμήματος
Τμήμα% (προεπιλογή: 50% επικάλυψη)
% - N_FFT: Αριθμός σημείων FFT (προεπιλογή: N_FFT = 256)
% - Fs: συχνότητα δειγματοληψίας σε Hz (προεπιλογή: fs = 1 Hz)
%
% Output: - PSD: Power φασματική πυκνότητα
% - Στ: Συχνότητα άξονα
ο καθένας μπορεί να μου πει ... παραμέτρους που θα πρέπει να αλλάξω για να δείτε το ψήφισμα;;;

 

[Guest]

Για να εξετάσει πώς οι αλλαγές που έχετε ψηφίσματος για την αλλαγή των παραμέτρων του σήματος που αναλύονται, δηλαδή ο τόνος συχνότητες, την πλάτη και SNR στο sg_cissoid λειτουργία.Με την αλλαγή αυτών των παραμέτρων, μπορείτε να βρείτε την ελάχιστη διαφορά στη συχνότητα που να επιτρέπει να δείτε τις αρμονικές στην έξοδο φάσμα εκτίμηση (PSD οικόπεδο) ως δύο διαφορετικές γραμμές.Στην πραγματική εφαρμογή που θα πρέπει να παρέχουν το κατώτατο όριο για την επίλυση του PSD ανάλυση - η πιθανότητα για σωστή ανάλυση εξαρτάται σε μεγάλο βαθμό από το όριο.Κατά την εργασία με παραμετρικές μεθόδους εκτίμησης που είναι δυνατόν να εξετάσει τις πολύπλοκες ρίζες του PSD φάσμα εκτίμηση αντί του PSD ίδιο οικόπεδο.

 

[seifer]

σας ευχαριστώ πολύ για την απάντησή σας πρώτα ... θα ήθελα να μάθω κάτι άλλο ... Πρέπει να δω αυτό σε 64 δεδομένα τιμές τότε 128 κλπ. .. έτσι νομίζω ότι χρειαζόμαστε ένα βρόχο για την παράμετρο στ δικαίωμα; .. . i και χρειάζονται ένα βρόχο που να δίνει τις κατάλληλες τυχαίων αριθμών ... μπορεί και να μου πει πώς θ ιβιβ κάνει αυτό; ... και αν αυτό είναι το σωστό πράγμα που πρέπει να κάνουμε

 

[Guest]

Στην πραγματικότητα,
ποτέ δεν συνεργάστηκε με SALP Εργαλειοθήκη, αλλά έκανα τα παρόμοια πράγματα με Signal Processing Toolbox.Φαντάζομαι ότι θα χρειαστείτε τρεις βρόχους:
Η πρώτη είναι για το σύνολο δεδομένων μήκος (Δημ. στον ορισμό του sg_cissoid λειτουργία),
το δεύτερο είναι η διαφορά μεταξύ των ήχων συχνότητες (αν θέλετε να εξετάσουν με ποιο τρόπο τα αποτελέσματα της αλλαγής σε διάφορα ήχους εκτοπίσματος) - είναι η διαφορά των δύο τιμών (με
την προϋπόθεση 2-Ήχος σήμα) στο διάνυσμα f,
και η τρίτη είναι για το SNR (αν θέλετε να εξετάσει το επίπεδο θορύβου και επίδραση).
Δεν κατάλαβα τι εννοούσατε με τυχαίους αριθμούς.Αν σήμαινε θόρυβο από ό, τι φαντάζομαι ότι sg_cissoid λειτουργία χρησιμοποιεί fuction όπως wgn από Signal Processing Toolbox, οπότε δεν χρειάζεται να κάνετε τον κόπο να παράγει θόρυβο σχετικά με το χέρι.Εσείς θα πρέπει να παρέχουν τη λειτουργία sg_cissoid με μια κατάλληλη τιμή της παραμέτρου SNR.

 

[seifer]

Μία ακόμη φορά thx πολύ για την απάντησή σας ... Θα κάνω και εγώ το σήμα proccesing εργαλειοθήκη ... αλλά έχει πολύ λίγες λειτουργίες για PSD ... και μπορεί να πει πώς μπορώ να βρω περισσότερα για αυτό; Εννοώ ότι για παράδειγμα την πρώτη φορά i ll εξετάσει PSD για 64 τιμές δεδομένων μετά ... ότι για 128 και στη συνέχεια να τηρούν το ψήφισμα ... έτσι ιβιβ ανάγκη για 64 διαφορετικές στ κωδικό μου ... έτσι; im δικαίωμα, εάν και πώς ... ll ότι κάνουμε βρόχο;Προστέθηκε μετά από 3 λεπτά:και κάτι ακόμα ... αυτό είναι το γ κώδικας για τη δεδομένη cissoid λειτουργία
λειτουργία varargout = sg_cissoid (β, στ, Δημ., varargin)?

Ελέγξτε τις παραμέτρους%
σφάλμα (nargchk (3,5, nargin))?
[msg, OSF, SNR] = chk_param (β, στ, Δημ., varargin {:});
σφάλμα (MSG)?

Δειγματοληψία% παράμετροι
fmax = 0.5?
fn = 2 * fmax?
T = 1 / (OSF * fn)?

% Setup άξονα
t = T * [0: Δημ.-1]?

% Preprocess δεδομένων
β = beta (

<img src="http://www.edaboard.com/images/smiles/icon_smile.gif" alt="Χαμόγελο" border="0" />

. "?
f = f (

<img src="http://www.edaboard.com/images/smiles/icon_smile.gif" alt="Χαμόγελο" border="0" />

?

Δημιουργήστε% σήμα
x = β * EXP (j * 2 * pi * f * t)?

Προσθήκη% θορύβου εάν είναι απαραίτητο
εάν (SNR ~ = inf)
y = add_noise (x, SNR)?
άλλος
y = x?
τέλος

% Εμφάνιση αποτελεσμάτων εξόδου επιχείρημα εφόσον δεν ορίζεται
εάν (nargout == 0)
σχήμα ( «όνομα», «Signal Generator: υπερτιθέμενους cissoids στο λευκό θόρυβο», «NumberTitle ',' off ')?
subplot (2,2,1)?
οικόπεδο (t, abs )?
xlabel ( 'Time t ')?
ylabel ( 'Μέγεθος - απόλυτη τιμή »)?
subplot (2,2,2)?
οικόπεδο (t, γωνία )?
xlabel ( 'Time t ')?
ylabel ( 'Μέγεθος - φάση »)?
subplot (2,2,3)?
οικόπεδο (t, Real )?
xlabel ( 'Time t ')?
ylabel ( 'Μέγεθος - πραγματικό μέρος ")?
subplot (2,2,4)?
οικόπεδο (t, IMAG )?
xlabel ( 'Time t ')?
ylabel ( 'Μέγεθος - φανταστικό μέρος »)?
άλλος
Αντιστοιχίστε% εξόδου επιχειρήματα
varargout (1) = y. "?
varargout (2) = t. '?
τέλος

επιστροφή?% Ελέγξετε παραμέτρους και setup defaults αν χρειάζεται
λειτουργία [msg, OSF, SNR] = chk_param (β, στ, Δημ., varargin)?

Δεν% μήνυμα λάθους
msg =''?

Λάβετε% προαιρετικών παραμέτρων
if (nargin> 3 & ~ isempty (varargin (1)))
OSF = varargin (1)?
άλλος
OSF = 1?
τέλος
if (nargin> 4 & ~ isempty (varargin (2)))
SNR = varargin (2)?
άλλος
SNR = inf?
τέλος

Ελέγξτε% παράμετρο
εάν (OSF <1)
msg = 'Over-δειγματοληψία παράγοντα OSF''''πρέπει να είναι τουλάχιστον της αξίας''1''.';
επιστροφή?
τέλος
εάν (μήκος (OSF)> 1)
msg = 'Παράμετρος OSF''''πρέπει να είναι μια βαθμωτή.'?
επιστροφή?
τέλος
εάν (μήκος (SNR)> 1)
msg = 'Παράμετρος''''SNR πρέπει να είναι μια βαθμωτή.'?
επιστροφή?
τέλος
εάν (κάθε (Δημ. <= 1) | μήκος (Δημ.) ~ = 1)
msg = 'Δεν μονοδιάστατες σήμα ορίζεται. "?
επιστροφή?
τέλος
if ((μήκος (β) ~ = μήκος (στ)))
msg = 'Παράμετρος φορείς''Β''και''στ''πρέπει να είναι ίδιο μήκος. "?
επιστροφή?
τέλος
fmax = 0.5?
εάν (κάθε (abs (στ)> fmax))
msg = 'Don''t υπερβαίνει Nyquist συχνότητας όριο | fmax | = 0.5 Hz. "?
επιστροφή?
τέλος

επιστροφή?
και μου είπε ότι μόνο στ και ισχνές επηρεάσει το PSD ... αλλά αν μου δοθεί μια διαφορετική OSF i παρατηρεί διαφορές σε οικόπεδο .. μπορεί και να μου πείτε γιατί;

 

[Guest]

Απόσπασμα:

Θα κάνω εγώ το σήμα proccesing εργαλειοθήκη ... αλλά έχει πολύ λίγες λειτουργίες για PSD ... και να πείτε πώς μπορώ να βρω περισσότερα για αυτό;