頻譜分析的仿真實例

信號的頻譜分析

a=10; %波信號振幅

fb1=3e6; %頻率fb1

fb2=4e6; %頻率fb2

fs=2.4*fb2; %抽樣頻率fs

T=1e-5; %調頻周期，也可設置抽樣點數

t=0:(1/fs):(T-1/fs);

N=length(t);

df=fs/N;

f=0:df:(N-1)*fs/N;

x=1/fs*a*(cos(2*pi*fb1.*t)+cos(2*pi*fb2.*t)); %線性調頻連續波信號

subplot(4,1,1);

plot(t,x);

title('波信號')

subplot(4,1,2);

stem(t,x);

title('抽樣後');

subplot(4,1,3);

stem(f,abs(fft(x,N)));

title('頻譜幅度');

subplot(4,1,4);

stem(f,angle(fft(x,N)));

title('頻譜相位');