Spectrale Analyse en SALP tollbox

S

seifer

Guest
nu iets over SALP TOOLBOX?
dit is de code van het signaal:
[y, t] = sg_cissoid (bèta, f, Dim, <OSF>, <SNR>)
%
% Parameters: - beta: Vector met complexe amplitudes
% - F: Vector met een frequentie waarden (-0,5 <f <0,5)
% - Afm: Aantal monsters
% - OSF: Oversampling factor (standaard: OSF = 1)
% - SNR: Signaal-ruisverhouding in [dB] (standaard: SNR = inf)
%
% Output: - y: Vector met signaal monsters
% - T: tijd-as
en dit is de code van Welch bijvoorbeeld:

[PSD, f] = psd_welch (y, <Win>, <N_Overlap>, <N_FFT>, <fs>)
%
% Parameters: - y: Vector met signaal monsters
% - Win: Het signaal is opgedeeld in K_seg overlappende
% Segmenten met een lengte die gelijk is aan de lengte van
%
In het venster en ramen elk segment
% (Standaard: 8 segmenten met Hamming venster)
% - N_Overlap: Aantal overlappende monsters uit segment
% Segment (standaard: 50% overlap)
% - N_FFT: Aantal FFT punten (standaard: N_FFT = 256)
% - Fs: Bemonstering frequentie in Hz (standaard: fs = 1 Hz)
%
% Output: - PSD: spectrale vermogensdichtheid
% - F: Frequentie as
kan iemand mij vertellen ... welke parameters moet ik zien dat de resolutie?

 
Om te onderzoeken hoe de resolutie verandert u de parameters van het signaal wordt geanalyseerd, namelijk de toon frequenties, hun amplitudes en SNR in sg_cissoid functie.Door het veranderen van deze parameters vindt u het minimale verschil in frequenties die het mogelijk maakt om te zien harmonischen in de output spectrum raming (PSD perceel) als twee verschillende lijnen.In het echte aanvraag moet u de resolutie voor de PSD-analyse - de waarschijnlijkheid van de juiste resolutie is sterk afhankelijk van de drempel.Bij het werken met parametrische ramingsmethoden het mogelijk de complexe wortels van PSD spectrum schatten in plaats van de PSD perceel zelf.

 
hartelijk dank voor uw reactie eerst ... Ik zou graag willen weten een ding ... ik nodig om dit te zien in 64 gegevenswaarden dan 128 etc. .. dus ik denk dat ik behoefte aan een loop voor parameter f toch? .. . en ik moet een lus die de passende willekeurige getallen ... kunt u mij vertellen hoe ik ll dat doen? ... en als dit is het juiste ding om te doen

 
Eigenlijk
heb ik nooit gewerkt met SALP Toolbox, maar ik maakte de soortgelijke dingen met Signal Processing Toolbox.Ik denk dat U moet drie lussen:
de eerste is voor de dataset lengte (Dim in de definitie van de sg_cissoid functie),
de tweede is voor het verschil tussen toon frequenties (als je wilt onderzoeken hoe de resultaten veranderingen op verschillende tonen verplaatsing) - is het verschil tussen de twee waarden (uitgaande van 2-tone-signaal) in vector f,
en de derde is voor de SNR (als je wilt onderzoeken of het geluidsniveau invloed als goed).
Ik niet begrijp wat u bedoelt met willekeurige getallen.Als U bedoelt lawaai dan denk ik dat sg_cissoid functie maakt gebruik van een functie zoals WGN uit Signal Processing Toolbox, zodat u niet hoeft te storen over het genereren van geluid handmatig.U zal slechts tot het leveren van de sg_cissoid functie met een passende waarde van de SNR parameter.

 
Nog een keer thx veel voor uw reactie ... ik kon gebruiken ook het signaal proccesing gereedschapskist ... maar het heeft weinig functies voor psd ... kan u vertellen hoe kan ik meer van dit? Bedoel ik dat bijvoorbeeld de eerste keer dat ik kom onderzoeken psd voor 64 gegevenswaarden ... na dat voor de 128 en vervolgens aan de resolutie ... dus ik kom moeten 64 verschillende f voor mijn code ... toch? en als im right ... hoe ik ll doen lus?Toegevoegd na 3 minuten:en nog een ding ... dit is de C-code voor de gegeven cissoid functie
functie varargout = sg_cissoid (bèta, f, Dim, varargin);

% Check parameters
fout (nargchk (3,5, nargin));
[msg, OSF, SNR] = chk_param (bèta, f, Dim, varargin {:});
fout (msg);

% Bemonstering parameters
Fmax = 0,5;
fn = 2 * Fmax;
T = 1 / (OSF * fn);

% Setup as
t = T * [0: Dim-1];

% Preprocess gegevens
beta = beta (

<img src="http://www.edaboard.com/images/smiles/icon_smile.gif" alt="Lachten" border="0" />

. ";
f = f (

<img src="http://www.edaboard.com/images/smiles/icon_smile.gif" alt="Lachten" border="0" />

;

Genereer signaal%
x = beta * exp (j * 2 * pi * f * t);

% Toevoegen lawaai indien nodig
if (SNR ~ = inf)
y = add_noise (x, SNR);
anders
y = x;
eindigen

% Zichtbare resultaten als er geen output argument is gedefinieerd
if (nargout == 0)
figuur ( 'naam', 'Signal Generator: Bovenop cissoids in witte ruis',' NumberTitle ',' off ');
subplot (2,2,1);
plot (t, abs (y));
xlabel ( 'Tijd t ');
ylabel ( 'Magnitude - absolute waarde');
subplot (2,2,2);
plot (t, hoek (y));
xlabel ( 'Tijd t ');
ylabel ( 'Magnitude - fase');
subplot (2,2,3);
plot (t, real (y));
xlabel ( 'Tijd t ');
ylabel ( 'Magnitude - reëel deel');
subplot (2,2,4);
plot (t, IMAG (y));
xlabel ( 'Tijd t ');
ylabel ( 'Magnitude - imaginaire deel');
anders
Wijs% output argumenten
varargout (1) =
y ';
varargout (2) = t. ";
eindigen

terugkeer;% Controleren parameters en setup defaults indien nodig
functie [msg, OSF, SNR] = chk_param (bèta, f, Dim, varargin);

% Geen foutmelding
msg ='';

% Get optionele parameters
if (nargin> 3 & ~ isempty (varargin (1)))
OSF varargin = (1);
anders
OSF = 1;
eindigen
if (nargin> 4 & ~ isempty (varargin (2)))
SNR = varargin (2);
anders
SNR = inf;
eindigen

% Check parameter
if (OSF <1)
msg = 'Over-bemonstering factor''OSF''moet ten minste waarde''1''.';
terugkeer;
eindigen
if (lengte (OSF)> 1)
msg = 'Parameter''OSF''moet een scalaire. ";
terugkeer;
eindigen
if (lengte (SNR)> 1)
msg = 'Parameter''SNR''moet een scalaire. ";
terugkeer;
eindigen
als (enige (Dim <= 1) | Lengte (Dim) ~ = 1)
msg = 'Geen eendimensionale signaal genoemd. ";
terugkeer;
eindigen
if ((lengte (bèta) ~ = lengte (f)))
msg = 'Parameter vectoren bèta''''en''F''moeten van dezelfde lengte. ";
terugkeer;
eindigen
Fmax = 0,5;
als (enige (abs (f)> Fmax))
msg = 'Don''t meer dan Nyquist-frequentie van | Fmax | = 0,5 Hz. ";
terugkeer;
eindigen

terugkeer;
U vertelde mij dat alleen f en dimmen van invloed zijn op de psd ... maar als ik een andere OSF i merkt verschillen in de plot .. kunt u mij vertellen waarom?

 
Citaat:

kon ik gebruik ook het signaal proccesing gereedschapskist ... maar het heeft weinig functies voor psd ... kan u vertellen hoe kan ik meer van dit?
 
You must log in or register to reply here.

Welcome to EDABoard.com

Sponsor

Back
Top