close all	%關(guān)閉所有的窗口
clear all	%清除所有的變量
%要求輸入k,從而計(jì)算總?cè)狱c(diǎn)數(shù)N=2^k
k=input('取樣點(diǎn)數(shù)=2^k,k=[14]   ');
if isempty(k),
   k=14;
end;
N=2^k;
%初始化時(shí)域分辨率dt（單位ms），
%頻域分辨率df（單位kHz）
%并滿足df*dt=1/N
dt=0.01;   
df=1.0/(N*dt);  
%計(jì)算時(shí)間寬度T（單位ms）
%系統(tǒng)帶寬Bs（單位kHz）
T=N*dt;
Bs=N*df/2;
%產(chǎn)生時(shí)間取樣點(diǎn)距陣t
%和頻域取樣點(diǎn)距陣f
t=[-T/2+dt/2:dt:T/2];
f=[-Bs+df/2:df:Bs];
%設(shè)置sin函數(shù)的基頻f0（單位kHz）
f0=1;
%產(chǎn)生sin函數(shù)抽樣距陣s
s=sin(2/3*pi*f0*t);
%做傅立葉變換，產(chǎn)生頻域抽樣距陣S
S=t2f(s,dt);
%計(jì)算出功率譜距陣P
P=S.*conj(S)/T;
%畫時(shí)域圖
figure(1)	%開窗口1
%設(shè)置窗口位置和大小
set(1,'position',[10,50,300,200])	
plot(t,s,'r');
grid on		%顯示網(wǎng)格
%設(shè)置坐標(biāo)范圍
axis([-2^10*dt/2,2^10*dt/2,min(s),max(s)])	
xlabel('t(ms)')	%給X軸加標(biāo)簽
ylabel('s(t)')	%給Y軸加標(biāo)簽
%畫頻譜曲線
figure(2)
set(2,'position',[350,50,300,200])
plot(f,abs(S),'g')
grid on
axis([-2*f0,2*f0,min(abs(S)),max(abs(S))])
xlabel('f(kHz)')
ylabel('S(f)')
%畫功率譜曲線
figure(3)
set(3,'position',[750,50,300,200])
plot(f,P,'b')
grid
axis([-2*f0,2*f0,min(P),max(P)])
xlabel('f(kHz)')
ylabel('P(f)')