스펙트럼 분석과 SALP tollbox

S

seifer

Guest
SALP 도구에 대해, 지금은 아무것도?
이것은 신호의 코드입니다 :
[y를, T는 '] = (베타로, F, 침침한 <OSF>, <SNR>)
sg_cissoid%
% 매개 변수 : - 베타 : 벡터 복잡한 amplitudes
함께% - f :
시스템 벡터 주파수 값을 (0.5 <F <0.5)
% - Dim : 샘플의 개수
% - OSF : 오버 요인 (기본값 : OSF = 1)
% - SNR을 : 신호에 잡음 비율 [dB의] (기본값 : SNR을 = inf를)
%
% 출력 : - y를 : 벡터 신호를 표본
% - 안 : 시간 축
이 예를 들어 웰치의 코드입니다 :

[의 PSD로, F] = psd_welch는 (y, <Win>, <N_Overlap>, <N_FFT>, <fs>)
%
% 매개 변수 : - y를 : 벡터 신호를 표본
% - 승 : 그 신호 K_seg로 나뉘어집니다 중복
길이와 동일한 길이의 세그먼트의 %
%는 윈도우와 윈도우 각 세그먼트
% (기본 : Hamming 창으로 8 세그먼트)
% - N_Overlap : 중복 표본의 개수 세그먼트에 이르기까지
% 세그먼트 (기본 : 50 % 중복)
% - N_FFT : FFT 포인트의 번호 (기본값 : N_FFT = 256)
% - 때 fs : Hz에서에서 샘플링 주파수 (기본 : 때 fs = 1 Hz에서)
%
% 출력 : -의 PSD : 파워 스펙트럼 밀도
% - f :
시스템 주파수 축
매개 변수 i는 해상도 변경해야 볼 수있는 날 ... 말할 수있는 사람이???

 
해상도를 변경하기 위해서는 어떻게 전파의 매개 변수를 변경해야 할 검사, 톤 주파수, 즉 자신의 amplitudes와 SNR을 sg_cissoid 함수에서 분석하고있다.당신이 주파수에서 출력 스펙트럼 추정 별개의 두 개 라인으로 (의 PSD 음모)에서 볼 수있는 최소한의 차이를 찾을 수 고조파 이러한 매개 변수를 변경하여.실제 응용 프로그램에서 분석을위한 당신의 PSD 해상도 임계값 공급해야한다 - 올바른 해결의 가능성을 강하게 임계값에 따라 달라집니다.가능성 자체가 음모의 PSD 대신 스펙트럼의 PSD 추정의 복잡한 뿌리를 검토 파라메 트릭 추정 방법으로 노력하고있습니다.

 
대단히 감사합니다 귀하의 답변을 먼저 ... 내가 또 다른 한가지는 알고 싶다 ... 내가 64이 데이터를 볼 필요는 다음 128
등 ..
난 그렇게
값을 나는 F 오른쪽 매개 변수에 대한 루프를 필요가있는 것 같아? .. . 그리고 난
유 게요 어떻게 말해줄 수는 적절한 무작위 숫자를 부여하는 루프 ...가 필요합니까? ... 그리고
만약이 할 수있는 올바른 것은

 
사실, 난 SALP 제작팁과 함께 일한 적이 있지만 신호 처리 도구와 유사한 것들을했다.난 당신이 필요합니다 3 루프 추측 :
첫 번째는 데이터의 길이 (Dim
the sg_cissoid 함수의 정의에),
() 벡터에로,
F 2
- 톤 신호를 가정 번째 한 톤 주파수 사이의 차이를 (해당되는 경우 그 결과가 어떻게 다른 음색 변위에서 변경) -이 두 가치의 차이가 검사를 원하는 것은
그리고 제 3 SNR을위한 하나입니다 (해당되는 경우는 물론, 소음 영향 조사).

당신이 임의의 숫자가 무슨 의미를 이해하지 않았다.보다 그 기능을 추측 sg_cissoid 소음을 의미하는 경우, 그래서 당신이 수동으로 생성하는 소음에 대해 걱정을하지 않아도 신호 처리 도구에서 WGN처럼 fuction을 사용합니다.당신은 SNR을 매개 변수의 적절한 가치와 sg_cissoid 기능을 공급해야한다.

 
한 번 더 들으 귀하의 답변에
많은 ...
저도 ... proccesing 도구를 사용할 수 있지만 신호의 PSD를위한 몇 가지있다 기능 아주 ... 유에 대한보다 자세한 내용을보실 수있습니다 이걸 어떻게 알 수 있나요? 내 말은 그 예를 들면 처음으로
줄게 내가 64의 PSD 데이터 값에 대한 검토 ... 후에 128을 누른
F ...
내 코드에 대한 권리를 다른 해상도를 ... 그래서 난 곧 64 필요 관찰? 그리고 만약 메신저 바로 ... 어떻게 곧 그 루프합니까?올린날짜3 분 후 :그리고 한 가지만 더 ... cissoid이 주어진 기능을 위해 C 코드입니다
기능 varargout = sg_cissoid (베타로, F, 침침한 varargin);

% 매개 변수를
확인하기오류 (nargchk (3,5, nargin));
[메시지, OSF, SNR을 '] = chk_param (베타로, F, 침침한
{:}); varargin
오류 (메시지);

% 샘플링 매개 변수
fmax = 0.5;
fn = 2 * fmax;
T는 = 1 / (OSF * fn);

% 설치 축
마 = T는 * [0 : Dim - 1];

% Preprocess 데이터
베타 = 베타 (

<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 = 베타 * 뉴스 (j를 * 2 * 파이 * F * T)에;

필요하다면 % 노이즈 추가
만약 (SNR을 ~ = inf를)
y를 = add_noise (엑스, SNR을);
그 밖의
y를 = x;


% 표시 결과가 없을 경우 출력 인자 정의됩니다
만약 (nargout == 0)
그림 ( '이름', '신호 생성기 :
백색 소음'에 겹쳐 cissoids, 'NumberTitle', '을');
subplot (2,2,1);
줄거리 (마, 복근는 (y));
xlabel ( '타임 마는 ');
ylabel ( 'Magnitude - 절대적인 가치');
subplot (2,2,2);
줄거리 (마, 각도는 (y));
xlabel ( '타임 마는 ');
ylabel ( 'Magnitude - 차');
subplot (2,2,3);
줄거리 (,는
(y)) 진짜 안;
xlabel ( '타임 마는 ');
ylabel ( 'Magnitude - 진짜 부분이');
subplot (2,2,4);
줄거리 (마,는 (y)) 심상;
xlabel ( '타임 마는 ');
ylabel ( 'Magnitude - 상상 부분이');
그 밖의
% 할당 출력 인자
varargout (1) =
당 ';
varargout (2) =
t를 ';


반환;% 매개 변수를 확인하고 필요한 경우 설치 프로그램의 기본값
함수 [, OSF, SNR을 '] = chk_param (베타로, F, 침침한 varargin) 오류;

% 아니 오류 메시지
메시지 = '';

% 선택적 매개 변수!
만약 (nargin> 3 & ~ IsEmpty (varargin (1)))
OSF = varargin (1);
그 밖의
OSF = 1;

만약 (nargin> 4 & ~ IsEmpty (varargin (2)))
SNR을 = varargin (2);
그 밖의
SNR을 = inf를;


% 체크 매개 변수
만약 (OSF <1)
메시지 = '오버 샘플링 OSF - 팩터 '' '' ''에서 1 이상의 가치를 ''.';되어야합니다
반환;

면 (길이 (OSF)> 1)
메시지 = '매개 변수 '' '' OSF 스칼라이어야합니다.';
반환;

면 (길이 (SNR)가> 1)
메시지 = '매개 변수 '' '' SNR을 스칼라이어야합니다.';
반환;

(임의 (Dim <= 1) | 길이 (Dim) ~ = 1)
메시지 = '한 - 차원 신호를 지정합니다.';
반환;

만약 ((길이 (베타) ~ = 길이 (바)))
메시지 = '매개 변수'베타 '와'F '로 벡터 동일한 길이의해야합니다.';
반환;

fmax = 0.5;
만약 (임의 (복근 (바)> fmax))
메시지 = 'don''t |
fmax | = 0.5 Hz의 Nyquist - 주파수 제한을 초과합니다.';
반환;


반환;
유에만 영향을 희미의 PSD F와 ...하지만 만약 내가 다른 OSF 제가 말해줬 음모에 차이가 관찰 유 .. 왜 그러는지 말해 줄 수 있나요?

 
인용구 :

나도 ... proccesing 도구를 사용할 수 있지만 신호의 PSD를위한 몇 가지있다 기능은 매우 ... 유에 대한보다 자세한 내용을보실 수있습니다 이걸 어떻게 알 수 있나요?
 

Welcome to EDABoard.com

Sponsor

Back
Top