MATLAB

Question

içindeki sinyal verisinden kaynaklanan frekansı belirleme Bir sensörden veri aldım ve frekansını bulmalıyım. Görünüşe göre, fft() gibi görünüyor, ancak MATLAB dokümanları sadece frekansların bir grafiğini nasıl göstereceklerini gösteriyor, oradan ne yapacağımı bilmiyorum.

İşte benim veri gibi görünüyor: gitmek

Answer 1

bir yolu bir fft kullanmak gerçekten de öyle. Fft, sinyalin frekans gösterimini sağladığından, maksimumı aramak istersiniz ve fft karmaşık bir sinyal olduğundan, önce mutlak değeri almak isteyeceksiniz. Endeks maksimum enerjiyle normalize frekansa karşılık gelecektir. Son olarak, eğer sinyaliniz bir ofsete sahipse, gösterdiğiniz ile olduğu gibi, fft'yi almadan önce bu ofsetden kurtulmak istersiniz, böylece DC bileşenini temsil eden bir menzil elde edemezsiniz. Ben olurdu bir satırda koymak açıklanan

Herşey: indexMax max FFT değeri bulunabilir endeksi ise

[maxValue,indexMax] = max(abs(fft(signal-mean(signal)))); 

. Not: indexMax öğesinden ilgili gerçek frekansa geçmek için, fft'in L uzunluğunu (sinyalin uzunluğu ile aynı) ve örnekleme frekansını (Fs) bilmeniz gerekir. Sinyal frekansı sonra olacak:

frequency = indexMax * Fs/L; 

Alternatif olarak, daha hızlı ve sahip sinyale bağlı olarak da oldukça iyi çalışıyor, senin sinyalin otokorelasyonunu atın:

autocorrelation = xcorr(signal); 

ve birinci maksimum sonrasında ortaya bulmak Otokorelasyonun merkez noktası. (Otokorelasyon, maksimum ortasındaki simetrik olacaktır.) En yüksek değeri bularak, kaymış sinyalin kendisinin daha çok ya da daha az göründüğü ilk yeri bulursunuz. Yani sinyalinin periyodunu bulursun. Döngünün bir katı tarafından değiştirilen sinyal her zaman kendi gibi görüneceğinden, bulduğunuz maksimum değerin sinyalin süresine karşılık geldiğinden ve katlarından biri olmadığından emin olmalısınız.

Sinyalinizdeki parazit nedeniyle, mutlak maksimum, periyodun kendisinin değil, döneminizin çoğunda meydana gelebilir. Bu gürültüyü hesaba katabilmek için, otokorelasyonun mutlak maksimumunu (otokorelasyon (uzunluk (otokorelasyon)/2 + 1) alırsınız ve daha sonra otokorelasyonun, ilk önce bu maksimum değerin% 95'inden daha büyük olduğunu bulursunuz. sinyalin ikinci yarısında geçen süre% 95,% 99 veya başka bir sayı, ne kadar gürültünün sinyalinizi bozduğuna bağlıdır.

GÜNCELLEME: Sinyalinizin "frekansı" ile kastettiğinizi varsaydığımı fark ettim. en enerji ile zift veya baz harmonik veya frekans, ancak bakmak istiyorum.Frekansla, sinyalin frekans gösterimi anlamına geldiyse, daha sonra bir ilk yaklaşım için, sadece bir FFT'nin absini çizmek için Enerjinin nerede olduğu fikri:

Bir abs'in neden olduğunu anlamak veya fft'in fazını temsil etmeyerek hangi bilgileri kaybettiğinizi anlamak istiyorsanız, tam olarak ne aldığınızı anlamak için DFT dönüşümü hakkında biraz daha fazla bilgi edinmek isteyebilirsiniz.

Answer 2

Ben

(indexMax-1) * Fs/L 

abs ilk öğesi (fft (x)) doğru akım (DC), ya da önyargı, olabilir ya sinyal ya X0'ın anlamına gerektiğini düşünüyorum. İkinci elementten (X1) saydık. Lütfen yanılıyorsam bildirin. Teşekkürler. Ben

clear all 
clc 
close all 
Fs = 1; 
T = 11 % Note this T is deliberately chosen , so that we have about 1.7 cycle of cosine singal 
t = 0:Fs:T; % T seconds 
L = length(t); % L is the length of sample sequence 
bias = 4 
signal = sin(t) + bias; 

[maxValue,indexMax] = max(abs(fft(signal-mean(signal)))); 

frequency_method1 = (indexMax-1) * Fs/(L-1); 
frequency_method2 = (indexMax-1) * Fs/L; 


number_of_cycles_method1 = frequency_method1*T 

number_of_cycles_method2 = frequency_method2*T 


subplot(2,1,1) 
plot(t,signal,'-or') ; grid on; 
legend('about 1.7 cycles of cosine signal') 
subplot(2,1,2) 
plot(abs(fft(signal-mean(signal))),'-xb'); grid on 
legend('abs of fft') 

number_of_cycles_method1 = 

    2 


number_of_cycles_method2 = 

    1.8333