參與評論

評論

登錄后參與評論

全部評論(1條)

• 

  qxxheart 2009-07-12 00:00:00

  給個例題，自己變一下就好 fs=1000;%設(shè)定采樣頻率 N=1024;%設(shè)定數(shù)據(jù)長度 i=0:N-1; t=i/fs; f=100;%設(shè)定正弦信號頻率 %生成正弦信號 x=sin(2*pi*f*t); subplot(231); plot(t，x);%作正弦信號的時域波形 axis([0，0.01，-1，1]); xlabel('t'); ylabel('y'); title('正弦信號時域波形'); grid; %進(jìn)行FFT變換并做頻譜圖 y=fft(x，N);%進(jìn)行fft變換 mag=abs(y);%求幅值 f=(0:N-1)*fs/N;%橫坐標(biāo)頻率的表達(dá)式為f=(0:M-1)*Fs/M; subplot(232); plot(f，mag);%做頻譜圖 axis([0，100，0，80]); xlabel('頻率(Hz)'); ylabel('幅值'); title('正弦信號幅頻譜圖'); grid; %求均方根譜 sq=abs(y); figure(1); subplot(233); plot(f，sq); xlabel('頻率(Hz)'); ylabel('均方根譜'); title('正弦信號均方根譜'); grid; %求功率譜 power=sq.^2; figure(1); subplot(234); plot(f，power); xlabel('頻率(Hz)'); ylabel('功率譜'); title('正弦信號功率譜'); grid; %求對數(shù)譜 ln=log(sq); figure(1); subplot(235); plot(f，ln); xlabel('頻率(Hz)'); ylabel('對數(shù)譜'); title('正弦信號對數(shù)譜'); grid; %用IFFT恢復(fù)原始信號 xifft=ifft(y); magx=real(xifft); ti=[0:length(xifft)-1]/fs; subplot(236); plot(ti，magx); xlabel('t'); ylabel('y'); title('通過IFFT轉(zhuǎn)換的正弦信號波形'); grid;

  贊(2)

  回復(fù)(0)

  評論

  評論

  登錄后參與評論