正弦序列FFT頻譜分析程序問題！！

beibei91 2010-07-11 14:44:33 461 瀏覽

1樓就是正弦包含頻率是20hz，20.5hz，40hz，采樣頻率fs是100hz，分析柵欄效應(yīng)，先是128個(gè)點(diǎn)fft，補(bǔ)零到512個(gè)點(diǎn)進(jìn)行fft，再512個(gè)點(diǎn)fft。程序是這樣的： N1=128;N2=512; fs=100;f1=20;f2=20.5;f3=40; n1=0:N1-1;n2=0:N2-1; xn1=sin(2*pi*f1*n1/fs)+sin(... 1樓就是正弦包含頻率是20hz，20.5hz，40hz，采樣頻率fs是100hz，分析柵欄效應(yīng)，先是128個(gè)點(diǎn)fft，補(bǔ)零到512個(gè)點(diǎn)進(jìn)行fft，再512個(gè)點(diǎn)fft。程序是這樣的： N1=128;N2=512; fs=100;f1=20;f2=20.5;f3=40; n1=0:N1-1;n2=0:N2-1; xn1=sin(2*pi*f1*n1/fs)+sin(2*pi*f2*n1/fs)+sin(2*pi*f3*n1/fs); xk11=fft(xn1，N1) mxk11=abs(xk11(1:N1/2)); figure(1); subplot(211);plot(n1，xn1); xlabel('n');title('x(n) 0<=n<127');axis([0，128，-3，3]); k1=(0:N1/2-1)*fs/N1; subplot(212) plot(k1，mxk11); xlabel('頻率單位Hz');title('X1(k)的幅度譜'); xn2=[xn1，zeros(1，N2-N1)]; xk12=fft(xn2，N2); mxk12=abs(xk12(1:N2/2)); figure(2); subplot(211);plot(n2，xn2); xlabel('n');title('x(n) 0<=n<=511');axis([0，512，-3，3]); k2=(0:N2/2-1)*fs/N2; subplot(212); plot(k2，mxk12); xlabel('頻率單位Hz');title('x1(k)補(bǔ)零后的幅度譜'); xn3=sin(2*pi*f1*n2/fs)+sin(2*pi*f2*n2/fs)+sin(2*pi*f3*n2/fs); xk2=fft(xn3，N2); mxk3=abs(xk2(1:N2/2)); figure(3); subplot(211);plot(n2，xn3); xlabel('n');title('x(n) 0<=n=511');axis([0，512，-3，3]); k3=(0:N2/2-1)*fs/N2; subplot(212); plot(k3，mxk3); xlabel('頻率單位Hz');title('512點(diǎn)有效數(shù)據(jù)的幅度譜'); 我看不懂的是 xk11=fft(xn1，N1) mxk11=abs(xk11(1:N1/2));（這個(gè)是什么意思？）和k1=(0:N1/2-1)*fs/N1;（為什么是二分之一得N1呢？）展開

參與評(píng)論

評(píng)論

登錄后參與評(píng)論

全部評(píng)論(2條)

萌貍仔甴 2010-07-12 00:00:00

diyi行是對(duì)XN1做fft的N1點(diǎn)變換第二條是求FFt變化后的幅值，括號(hào)中說的是采樣點(diǎn)，也就是說求對(duì)應(yīng)采樣點(diǎn)的幅值第三條是求頻率軸，除以二取整是求正半軸，當(dāng)然也可以取負(fù)值，這就看你對(duì)結(jié)果的要求了。

贊(19)

回復(fù)(0)

評(píng)論

評(píng)論
登錄后參與評(píng)論
夢(mèng)色未央 2016-07-28 14:05:04

因?yàn)镹個(gè)樣點(diǎn)的信號(hào)經(jīng)過fft以后變成N個(gè)樣點(diǎn)的頻譜，這個(gè)頻譜是關(guān)于第N/2+1樣點(diǎn)左右對(duì)稱的，所以真正有用的頻譜數(shù)據(jù)只有前面一半，后面一半是鏡像。mxk11是對(duì)前N/2個(gè)樣點(diǎn)取幅度譜，其實(shí)應(yīng)該是取1：N1/2+1，你這里少取了一個(gè)點(diǎn)。具體為什么會(huì)鏡像請(qǐng)看數(shù)字信號(hào)處理DFT章節(jié)。

贊(14)

回復(fù)(0)

評(píng)論

評(píng)論
登錄后參與評(píng)論

登錄或新用戶注冊(cè)

微信登錄
密碼登錄
短信登錄

請(qǐng)用手機(jī)微信掃描下方二維碼
快速登錄或注冊(cè)新賬號(hào)

微信掃碼，手機(jī)電腦聯(lián)動(dòng)

注冊(cè)登錄即表示同意《儀器網(wǎng)服務(wù)條款》和《隱私協(xié)議》

熱門問答

正弦序列FFT頻譜分析程序問題??！: 1樓就是正弦包含頻率是20hz，20.5hz，40hz，采樣頻率fs是100hz，分析柵欄效應(yīng)，先是128個(gè)點(diǎn)fft，補(bǔ)零到512個(gè)點(diǎn)進(jìn)行fft，再512個(gè)點(diǎn)fft。程序是這樣的： N1=128;N2=512; fs=100;f1=20;f2=20.5;f3=40; n1=0:N1-1;n2=0:N2-1; xn1=sin(2*pi*f1*n1/fs)+sin(... 1樓就是正弦包含頻率是20hz，20.5hz，40hz，采樣頻率fs是100hz，分析柵欄效應(yīng)，先是128個(gè)點(diǎn)fft，補(bǔ)零到512個(gè)點(diǎn)進(jìn)行fft，再512個(gè)點(diǎn)fft。程序是這樣的： N1=128;N2=512; fs=100;f1=20;f2=20.5;f3=40; n1=0:N1-1;n2=0:N2-1; xn1=sin(2*pi*f1*n1/fs)+sin(2*pi*f2*n1/fs)+sin(2*pi*f3*n1/fs); xk11=fft(xn1，N1) mxk11=abs(xk11(1:N1/2)); figure(1); subplot(211);plot(n1，xn1); xlabel('n');title('x(n) 0<=n<127');axis([0，128，-3，3]); k1=(0:N1/2-1)*fs/N1; subplot(212) plot(k1，mxk11); xlabel('頻率單位Hz');title('X1(k)的幅度譜'); xn2=[xn1，zeros(1，N2-N1)]; xk12=fft(xn2，N2); mxk12=abs(xk12(1:N2/2)); figure(2); subplot(211);plot(n2，xn2); xlabel('n');title('x(n) 0<=n<=511');axis([0，512，-3，3]); k2=(0:N2/2-1)*fs/N2; subplot(212); plot(k2，mxk12); xlabel('頻率單位Hz');title('x1(k)補(bǔ)零后的幅度譜'); xn3=sin(2*pi*f1*n2/fs)+sin(2*pi*f2*n2/fs)+sin(2*pi*f3*n2/fs); xk2=fft(xn3，N2); mxk3=abs(xk2(1:N2/2)); figure(3); subplot(211);plot(n2，xn3); xlabel('n');title('x(n) 0<=n=511');axis([0，512，-3，3]); k3=(0:N2/2-1)*fs/N2; subplot(212); plot(k3，mxk3); xlabel('頻率單位Hz');title('512點(diǎn)有效數(shù)據(jù)的幅度譜'); 我看不懂的是 xk11=fft(xn1，N1) mxk11=abs(xk11(1:N1/2));（這個(gè)是什么意思？）和k1=(0:N1/2-1)*fs/N1;（為什么是二分之一得N1呢？）展開

頻譜分析儀和FFT信號(hào)分析儀的相關(guān)問題: 我應(yīng)該用頻譜分析儀去測(cè)某喇叭在單一頻率聲源輸入下，是否會(huì)產(chǎn)生高次諧波；但是我沒有頻譜分析儀，有FFT信號(hào)分析儀，我今天研究半天，哪位高手幫幫忙，可以用后者代替前者來進(jìn)行實(shí)驗(yàn)嗎...盡量詳細(xì)一點(diǎn)，真的很急，謝謝啦

matlab中如何對(duì)正弦函數(shù)的采樣序列進(jìn)行頻譜分析: clear; Fs=100; N=64; t=(0:1/Fs:1); y=sin(2*pi*10*t); stem(t(1:N)，y(1:N)); 用上面的語句采樣，之后用Y=fft（）進(jìn)行FFT計(jì)算頻譜，那括號(hào)中應(yīng)該是什么？還有FFT之后如何將橫坐標(biāo)轉(zhuǎn)換成頻率。謝謝！

fft頻譜分析頻率誤差:

FFT和頻譜分析的關(guān)系:

示波器怎么用fft實(shí)現(xiàn)頻譜分析:

求助，F(xiàn)FT頻譜分析，怎樣將頻譜分析振幅轉(zhuǎn)為真實(shí)振幅:

FFT進(jìn)行頻譜分析，可以應(yīng)用在什么方面:

如何對(duì)wav文件進(jìn)行fft頻譜分析: 現(xiàn)在正用matlab做一wav文件分析（數(shù)據(jù)時(shí)從一個(gè)milling mchine上采集的），對(duì)信號(hào)處理這里完全不懂。。從整個(gè)文件聽感覺聲音是由很多階段的，我想應(yīng)該頻率都不同吧。。我自己把這個(gè)大文件分成了8份，可是用wavread分析出來的Fs都是44100，這怎么可能呢？我想畫... 現(xiàn)在正用matlab做一wav文件分析（數(shù)據(jù)時(shí)從一個(gè)milling mchine上采集的），對(duì)信號(hào)處理這里完全不懂。。從整個(gè)文件聽感覺聲音是由很多階段的，我想應(yīng)該頻率都不同吧。。我自己把這個(gè)大文件分成了8份，可是用wavread分析出來的Fs都是44100，這怎么可能呢？我想畫出各段的時(shí)間頻率圖。。。這樣如果一旦機(jī)器出現(xiàn)故障，這個(gè)時(shí)間段的頻率必定會(huì)發(fā)生波動(dòng)就可以知道機(jī)器是不是在正常工作了。。不知道我的想法對(duì)不對(duì)，還有這種曲線要怎么模擬呢？請(qǐng)高手指教！非常感謝！展開

快速fft采集音頻信號(hào)，如何進(jìn)行頻譜分析: 用單片機(jī)采集音頻信號(hào)，將音頻信號(hào)進(jìn)行傅立葉轉(zhuǎn)化后存在一個(gè)數(shù)組內(nèi)，那然后怎么處理這些數(shù)據(jù)呢，才能得到對(duì)這個(gè)音頻信號(hào)的進(jìn)行頻譜分析呢，想知道它是不是失真了，，請(qǐng)各位大俠幫助!!!

示波器FFT進(jìn)行一段聲音的頻譜分析，如圖所示: 如何判定在300HZ-4000Hz之間的諧波分量的相對(duì)聲壓級(jí)在脈沖頻率處幅度的正負(fù)15dB范圍內(nèi)？

怎么用FFT對(duì)正余弦函數(shù)進(jìn)行頻譜分析:

用DFT/FFT進(jìn)行頻譜分析會(huì)有那些問題，如何減少？:

VHDL序列信號(hào)發(fā)生器問題: 編好的程序如下：LIBRARYIEEE;USEIEEE.STD_LOGIC_1164.ALL;ENTITYMY_100110ISPORT(CLK:INSTD_LOGIC;Z:OUTSTD_LOGIC);ENDMY_100110;ARCHITECTURERTLOFMY_100110ISTYPESTATE_TYPEIS(S0，... 編好的程序如下： LIBRARY IEEE; USE IEEE.STD_LOGIC_1164.ALL; ENTITY MY_100110 IS PORT( CLK:INSTD_LOGIC; Z :OUTSTD_LOGIC ); END MY_100110; ARCHITECTURE RTL OF MY_100110 IS TYPE STATE_TYPE IS(S0，S1，S2，S3，S4，S5); SIGNAL CURRENT_STATE，NEXT_STATE:STATE_TYPE; BEGIN SYNCH: PRECESS BEGIN WAIT UNTIL CLK'EVENT AND CLK='1'; CURRENT_STATE<=NEXT_STATE; END PROCESS; STATE_TRANS:PROCESS(CURRENT_STATE) BEGIN CASE CURRENT_STATE IS WHEN S0=> NEXT_STATE<=S1; Z<='1'; WHEN S1=> NEXT_STATE<=S2; Z<='0'; WHEN S2=> NEXT_STATE<=S3; Z<='0'; WHEN S3=> NEXT_STATE<=S4; Z<='1'; WHEN S4=> NEXT_STATE<=S5; Z<='1'; WHEN S5=> NEXT_STATE<=S0; Z<='0'; END CASE; END PROCESS; END RTL; 用quartusⅡ運(yùn)行有5個(gè)錯(cuò)誤，限于水平有限，求大神幫助修改！展開