马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本人新手,在做关于气动噪声方面有点问题,希望能跟大家探讨一下。 1、 我在fluent 中取到了压力脉动,导入virtuallab中,做一次FFT变换,再做一次IFFT,结果两次的结果不一样。这个是什么原因引起的呢?(用matlab 做的两次结果基本一样的。) 2、 做两次变换,我原本的思路是,先把时域信号转化成频域(噪声计测量噪声,它的A计全是在频域计全,而不是在时域)。在频域中进行A计全。再把A计全的结果转化成时域信号。这样得到的时域噪声值一般是一个波动值。对这个波动值进行能量平均,最后得到的就是时域噪声。(我所理解的声级计测量噪声的原理,可能理解有误,勿怪) 3、 想法是美好的,现实是残酷的。。。在virtual lab 中,我把时域信号转化成频域之后,只能对这个结果进行A计全(就是作图那时能用下,对结果没影响),找不到能把它A计全之后,再转化为时域的方法。(有没有高手会这个的。。求助) 4、 基于这种情况,我就琢磨这用matlab编程试试,结果,做出来的频域结果与virtuallab得出的频域结果相差挺大的。(特别是高频的地方,matlab做出来的接近一条直线了)你所,virtual lab 有个什么经验修正,那也不至于差这么多吧。附上matlab程序。 5.傅里叶变换前面的系数根号2或者根号2/2,这个对傅里叶正正变换有何影响吗? file:///C:/Users/liu/AppData/Local/Temp/msohtmlclip1/01/clip_image002.jpg file:///C:/Users/liu/AppData/Local/Temp/msohtmlclip1/01/clip_image004.jpg fs=20000; x=VarName2; t=linspace(0.1,0.2,2000); plot(t,x); title('原时域噪声');xlabel('时间/s');ylabel('噪声/db');grid; N=length(x);%采样长度 f=[0:N-1]*fs/N;%FFT后各个频率 x1=fft(x); magx=abs(x1); magxj=abs(x1)*2/N;%幅值校正 magxjdb=20*log10(0.5*10^5*magxj); figure plot(f,magxjdb); title('频域噪声');xlabel('频率/Hz');ylabel('噪声/db');grid;
axis([0,10000,0,200]);
|