clear all;										% 清所有的變量
clf;											% 清屏
N=1224;										% 設(shè)置數(shù)據(jù)點(diǎn)數(shù)
f1=200;										% 設(shè)置信號(hào)1頻率
f2=700;										% 設(shè)置信號(hào)2頻率
f3=1100;										% 設(shè)置信號(hào)3頻率
dt=1/1024;										% 設(shè)置采樣時(shí)間
sample=1/dt;									% 設(shè)置采樣頻率
randn('state',sum(100*clock));						% 產(chǎn)生隨機(jī)噪聲
dis=2*rand(1,N);
dis=dis-mean(dis);								% 將隨機(jī)噪聲的均值設(shè)置成0
dis=2*dis;										% 擴(kuò)大隨機(jī)噪聲的幅度
for k=1:N;										% 產(chǎn)生待檢測(cè)的信號(hào)
    w(k)=1*sin(2*pi*f1*k*dt)+1*sin(2*pi*f2*k*dt)+1*sin(2*pi*f3*k*dt)+dis(k);
end

figure(1);										% 計(jì)算待檢測(cè)信號(hào)的均勻采樣頻譜
fw=fft(w);
s=abs(fw).^2;
f=(1:N-1)*sample/N;
plot(f,s(2:N));
xlabel('Frequency /Hz');							% 橫坐標(biāo)說(shuō)明
ylabel('Amplitude');								% 縱坐標(biāo)說(shuō)明
title('DFT of uniformly sampling signal');				% 標(biāo)題說(shuō)明

t=tnonunif(N,1,1);								% 調(diào)用非均勻采樣時(shí)間函數(shù)
for k=1:N;										% 產(chǎn)生非均勻采樣數(shù)據(jù)
    w2(k)=1*sin(2*pi*f1*t(k)*dt)+1*sin(2*pi*f2*t(k)*dt)+1*sin(2*pi*f3*t(k)*dt)+dis(k);
end

for k=1:N;										% 計(jì)算非均勻采樣數(shù)據(jù)的頻譜
    s2(k)=0;
    for n=1:N;									% 計(jì)算非均勻采樣數(shù)據(jù)的傅立葉變換
        s2(k)=s2(k)+w2(n)*(cos(2*pi*t(n)*(k-1)/N)-j*sin(2*pi*t(n)*(k-1)/N));
    end
    s2(k)=abs(s2(k)).^2;
end

figure(2);										% 畫出非均勻采樣數(shù)據(jù)的頻譜
f=(1:N-1)*sample/N;
plot(f,s2(2:N));
xlabel('Frequency /Hz');
ylabel('Amplitude');
title('DFT of additive random sampling signal');

f=(1:N-1)*sample/N;
plot(f,s2(2:N));
xlabel('Frequency /Hz');
ylabel('Amplitude');
title('DFT of additive random sampling signal');