Tulis program dahulu:
clc
jelas
CU=[]; %%Import data anda di sini
CU_f=fft(CU); %% melakukan transformasi FFT kepada domain frekuensi
angka(1)
subplot(2,1,1)
plot(CU);
subplot(2,1,2)
plot(abs(CU_f))
**********************************
CU_f ialah spektrogram, tetapi titik yang anda nyatakan sepadan dengan 50Hz perlu dicari berdasarkan proses percubaan anda. Input data ke MATLAB hanyalah beribu-ribu mata, yang hanya mempunyai makna matematik dan tiada makna fizikal. Pertama sekali, anda telah mencuba beribu-ribu mata ini, jadi apakah kekerapan pensampelan anda? Hanya orang yang melakukan eksperimen yang mengetahui perkara ini, bukan MATLAB sendiri. Andaikan bahawa kekerapan pensampelan ialah Fs. Kemudian dalam gambar CU_f, frekuensi sebenar sepadan dengan titik terakhir ialah Fs (frekuensi sebenar sepadan dengan titik pertama ialah 0Hz). Kemudian anda boleh mencari titik yang sepadan dengan 50Hz berdasarkan hubungan berkadar linear.
Gunakan MATLAB untuk melaksanakan transformasi Fourier:
Pengguna memasukkan fungsi secara sewenang-wenangnya, dan kemudian mengeluarkan fungsi transformasi Fourier bagi fungsi tersebut, dan kemudian mengeluarkan frekuensi amplitud.
x=sin(2*pi*t); %Masukkan sebarang fungsi.
y=fft(x); %Fourier transform function.
plot(abs(y)); % kekerapan amplitud.
Fungsi mewakili kesesuaian antara setiap nilai input dan nilai output unik. Hubungan ini menyebabkan setiap elemen dalam satu set sepadan dengan elemen unik dalam set lain (mungkin sama). Simbol piawai untuk nilai keluaran x sepadan dengan nilai input dalam fungsi f ialah f(x). Set yang mengandungi semua nilai input fungsi dipanggil domain fungsi, dan set yang mengandungi semua nilai output dipanggil julat. Jika kita mentakrifkan konsep pemetaan dahulu, kita boleh mentakrifkan fungsi sebagai .
Transformasi Fourier boleh menyatakan fungsi yang memenuhi syarat tertentu sebagai gabungan linear fungsi trigonometri (fungsi sinus dan/atau kosinus) atau kamirannya. Terdapat banyak variasi transformasi Fourier dalam bidang penyelidikan yang berbeza, seperti transformasi Fourier berterusan dan transformasi Fourier diskret.
% Jangan lupa beri saya mata, [senyum lebar]
kosongkan semua
ts=0.001; % Tempoh pensampelan
t=0:ts:20; % Urutan masa
y=sin(t)+0.5*sin(2*t)+0.2*sin(6*t);
figura
plot(t,y)
tajuk('Original Singal')
xlabel('Masa (s)')
ylabel('Magnitud')
Fs=1/ts; % Kekerapan pensampelan
L=panjang(y);
NFFT = 2^nextpow2(L); % Kuasa seterusnya 2 dari panjang y
Y = fft(y,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);
% Plot spektrum amplitud satu sisi.
figura
plot(f,2*abs(Y(1:NFFT/2+1)))
tajuk('Spektrum Amplitud Satu Sisi y(t)')
xlim([0,3])
xlabel('Frekuensi (Hz)')
ylabel('|Y(f)|')
Atas ialah kandungan terperinci Bagaimana pula dengan melakukan Fourier Transform menggunakan Matlab?. Untuk maklumat lanjut, sila ikut artikel berkaitan lain di laman web China PHP!