利用matlab繪制正弦信號的頻譜圖并做相關(guān)分析

1、利用Matlab繪制正弦信號的頻譜圖并做相關(guān)分析一、作業(yè)要求：1、信號可變（信號的賦值、相位、頻率可變）；2、采樣頻率fs可變；3、加各種不同的窗函數(shù)并分析其影響；4、頻譜校正；5、頻譜細(xì)化。二、采用matlab編寫如下程序：clearclffs=100N=1024%采樣頻率和數(shù)據(jù)點(diǎn)數(shù)A=20B=30C=0.38n=0:N1t=nfs%時(shí)間序列x=Asin(2piBtC)%信號y=fft(xN)%對信號進(jìn)行傅里葉變換yy=abs(y)%

2、求得傅里葉變換后的振幅yy=yy2N%幅值處理f=nfsN%頻率序列subplot(331)plot(fyy)%繪出隨頻率變化的振幅xlabel(頻率itHz)ylabel(振幅)title(圖1：fs=100，N=1024)gridon%兩種信號疊加，x=Asin(2piBtC)2Asin(2pi1.5Bt2.5C)%信號y=fft(xN)%對信號進(jìn)行傅里葉變換yy=abs(y)%求得傅里葉變換后的振幅yy=yy2N%幅值處理f=nf

3、sN%頻率序列subplot(332)plot(fyy)%繪出隨頻率變化的振幅xlabel(頻率itHz)ylabel(振幅)title(圖2：fs=100N=1024，兩種信號疊加)gridon%加噪聲之后的圖像x=Asin(2piBtC)28rn(size(t))y=fft(xN)yy=abs(y)yy=yy2N%幅值處理subplot(333)plot(f(1:N2.56)yy(1:N2.56))xlabel(頻率itHz)yla

4、bel(振幅)title(圖3：fs=100N=1024混入噪聲)y=fft(xN)%對信號進(jìn)行快速傅里葉變換yy=abs(y)%求取傅里葉變換的振幅yy=yy2N%幅值處理f=nfsNsubplot(337)plot(f(1:N2.56)1.852yy(1:N2.56))%繪出隨頻率變化的振幅xlabel(頻率itHz)ylabel(振幅)title(圖7：fs=100，N=1024加海明窗函數(shù))gridon%加漢寧窗函數(shù)后的頻譜fs

5、=100N=1024n=0:N1t=nfsx=Asin(2piBtC)%信號window=hanning(N)%生成漢寧窗函數(shù)x=x.window%加窗函數(shù)y=fft(xN)%對信號進(jìn)行快速傅里葉變換yy=abs(y)%求取傅里葉變換的振幅yy=yy2N%幅值處理f=nfsNsubplot(338)plot(f(1:N2.56)2yy(1:N2.56))%繪出隨頻率變化的振幅xlabel(頻率itHz)ylabel(振幅)title(圖