|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本帖最后由 chonwang 于 2010-12-22 14:57 编辑
Fluent acoustic 求得某點的聲壓
再利用MATLAB FFT分析
圖:(上)原始聲壓訊號
(中)FFT後聲壓訊號
(下)將FFT訊號轉換成dB(SPL)
MATLAB語法
----------------------------------------------------------------------------
clc;clear;
Timeset=0.008725; %時間間隔
Fstar=0; %頻譜圖頻率起點
Fend=10; %頻譜圖頻率終點
Fscal=0.01;%頻譜圖頻率最小尺度
SR=1/Timeset;
data = xlsread('test120steps_from885(0.008725s)x80_20101221.xls'); %讀取DATA
t = data(:,1);
y = data(:,2);
subplot(3,1,1)
plot(t,y)
title('原始訊號')
xlabel('Time(s)')
L = length(y);
NFFT = SR/Fscal;
Y = fft(y,NFFT)/L;
f = SR/2*linspace(0,1,NFFT/2+1);
subplot(3,1,2)
plot(f,2*abs(Y(1:NFFT/2+1)))
axis([Fstar,Fend,-inf,inf])
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')
P = 2*abs(Y(1:NFFT/2+1));
dB = 20*(log10(P/(2*0.00001)));
subplot(3,1,3)
plot(f,(dB))
axis([Fstar,Fend,-inf,inf])
xlabel('Frequency (Hz)')
ylabel('SPL(dB)')
----------------------------------------------------------------------------
問題:
1.聲壓訊號 只抓風機約轉三圈的訊號,FFT後是否不夠準確?
2.圖二,利用複製40倍,是不是會出現很多雜訊?
|
-
圖一,擷取自ANSYS fluent的聲壓訊號
-
圖二,將原始訊號複製x40
|