先寫程式:
clc
clear
CU=[]; %%這裡導入你的資料
CU_f=fft(CU);%%做FFT轉換到頻域
figure(1)
subplot(2,1,1)
plot(CU);
subplot(2,1,2)
plot(abs(CU_f))
***************************
#CU_f就是頻譜圖,但你說的對應50Hz的點,需要結合你的實驗過程來找了。 MATLAB輸入的資料只不過是上千個點,它們只有數學意義,沒有物理意義。首先,這上千個點是你取樣採到的,那你的取樣頻率是多少?這只有你做實驗的人才知道,MATLAB本身不知道。假設取樣頻率是Fs。那麼CU_f這張圖裡,最後一個點對應的實際頻率就是Fs,(第一個點對應的實際頻率是0Hz)。然後你再根據線性比例的關係,找到50Hz對應的那個點就行了。
#用MATLAB 實作傅立葉變換:
使用者任意輸入一個函數,然後,輸出函數的傅立葉變換函數,然後輸出振幅頻率 。
x=sin(2*pi*t); %任意輸入一個函數。
y=fft(x); % 傅立葉轉換函數。
plot(abs(y)); %振幅頻率。
傅立葉轉換能將滿足某一條件的某個函數表示成三角函數(正弦和/或餘弦函數)或它們的積分的線性組合。在不同的研究領域,傅立葉變換具有多種不同的變體形式,如連續傅立葉變換和離散傅立葉變換。
用matlab驗證傅立葉變換性質怎麼寫程式
% 不要忘記給我分, [一個大寫的微笑]
clear all
ts=0.001; % Sampling period
t=0:ts:20; % Time sequence
y=sin(t) 0.5*sin(2*t) 0.2*sin(6*t);
figure
plot(t,y)
title('Original Singal')
xlabel('Time (s)')
ylabel('Magnitude')
Fs=1/ts; % Sampling frequency
L=length(y);
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(y,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2 1);
% Plot single-sided amplitude spectrum.
figure
plot(f,2*abs(Y(1:NFFT/2 1)))
title('Single-Sided Amplitude Spectrum of y(t)')
xlim([0,3])
###xlabel('Frequency (Hz)')### ###ylabel('|Y(f)|')###以上是使用Matlab進行傅立葉變換的方法如何?的詳細內容。更多資訊請關注PHP中文網其他相關文章!