資源描述:
《典型函數(shù)的頻譜》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在工程資料-天天文庫。
1、典型函數(shù)的頻譜(矩形窗函數(shù),漢寧窗函數(shù),直線,階躍函數(shù),6函數(shù),方波,三角波等),如圖13?18所示。矩形窗函數(shù)的時域波形圖°010.550100150200250矩形窗函數(shù)頻域波形圖圖13躍數(shù)的時域波形圖0.8111iiii111L-111110.61111?mhb1IIII10.4IIII-IL1LIIII10.2nIIII.」L11IIIIIIII1ii050100150200250300豳數(shù)的頻域波形圖圖14方波的時域波形圖圖15漢寧窗函數(shù)的時域波形圖150漢寧窗函數(shù)頻域波形圖10050100頻
2、率150圖163002階躍函數(shù)的時域波形圖1.510.550100150200250階躍函數(shù)的頻域波形圖200100050100150頻率0三角波的時域波形圖00.020.040.060.080.10.120.140.160.180.2圖1710.50-0.5三角波的頻域波形圖050100150頻率-18060S40200圖18此部分MATLAB代碼如下:t=0:0.001:0.2;N=256;FS=300;w=boxcar(N);%產(chǎn)生信號figure;plot(w);titleC矩形窗函數(shù)的時域波形
3、圖axis([0,260,0,21);gridon;y=ffi(w,N);%FFT運算mag=abs(y);%lft幅值f=(O:length(y)-1)*FS/length(y);figure;plot(f(l:N/2),mag(l:N/2));%輸出FS/2點幅頻譜圖titleC矩形窗函數(shù)頻域波形圖*);grid;xlabel(瀕率J;ylabelC幅值);t=0:0.001:0.2;N=256;FS=300;w=hanning(N);%產(chǎn)生信號figure;plot(w);titlcC7X寧窗函數(shù)
4、的吋域波形圖*);gridon;y=ffi(w,N);%FFT運算mag=abs(y);%取幅值f=(0:length(y)-1)*FS/length(y);figure;plot(f(l:N/2),mag(l:N/2));%輸出FS/2點幅頻譜圖titled漢寧窗函數(shù)頻域波形圖工gridon;xlabef頻率J;ylabel(幅值');t=0:0.001:0.2;N=256;FS=300;w=l;%產(chǎn)牛信號y=fft(w,N);%FFT運算mag=abs(y);%取幅值f=(O:lcngth(y)-l
5、)*FS/lcngth(y);figure;plot(f(l:N/2),mag(l:N/2));%輸出FS/2點幅頻譜圖title?直線頻域波形圖);gridon;xlabel(,Frcquency(Hz),);ylabel(,Magnitude');%階躍函數(shù)的頻域波圖clc;clf;t=0:0.001:0.2;N=256;FS=300;w=oncs(l,N);%產(chǎn)生信號figure;plot(w);title。階躍函數(shù)的時域波形圖J;gridon;y=fft(w,N);%FFT運算mag=abs(y
6、);%取幅值f=(O:length(y)-1)*FS/length(y);figure;plot(f(l:N/2),mag(l:N/2));%輸出FS/2點幅頻譜圖titlg階躍函數(shù)的頻域波形圖);gridon;xlabcf頻率lylabcl(嘲值t=0:0.001:0.2;N=256;FS=3OO;w=zeros(l,N);w(l)=l;%產(chǎn)生信號figure;plot(w);gridon;title?函數(shù)的時域波形圖);y=fft(w,N);%FFT運算mag=abs(y);%取幅值f=(O:lcn
7、gth(y)-l)*FS/lcngth(y);figure;plot(f(l:N/2),mag(l:N/2));%輸出FS/2點幅頻譜圖title(f8函數(shù)的頻域波形圖);gridon;xlabel(瀕率J;ylabelC幅值工t=0:0.001:0.2;N=256;FS=300;w=square(2*pi*50*t);%產(chǎn)生信號figure;plot(t,w);titled方波的時域波形圖工axis([0,0.2,-0.2,1.2]);gridon;y=fft(w,N);%FFT運算mag=abs(y
8、);%取幅值f=(0:length(y)-1)*FS/1ength(y);figure;plot(f(l:N/2),mag(l:N/2));%輸出FS/2點幅頻譜圖titled方波的頻域波形圖工gridon;xlabef頻率J;ylabel(幅值');t=0:0.001:0.2;N=256;FS=300;w=sawtooth(2*pi*50*t,0.5);figure;plot(t,w);gridon;title(三角波的吋域波形圖);%產(chǎn)