**3.白噪聲****************%
fs=10;%設定采樣頻率
t=-5:0.1:5;
x=zeros(1,100);
x(50)=100000;
figure(3);
subplot(231);
plot(t(1:100),x);%作白噪聲的時域波形
xlabel('t');
ylabel('y');
title('白噪聲時域波形');
grid;
%進行FFT變換并做頻譜圖
y=fft(x);%進行fft變換
mag=abs(y);%求幅值
f=(0:length(y)-1)'*fs/length(y);%進行對應的頻率轉換
figure(3);
subplot(232);
plot(f,mag);%做頻譜圖
xlabel('頻率(Hz)');
ylabel('幅值');
title('白噪聲幅頻譜圖');
grid;
%求均方根譜
sq=abs(y);
figure(3);
subplot(233);
plot(f,sq);
xlabel('頻率(Hz)');
ylabel('均方根譜');
title('白噪聲均方根譜');
grid;
%求功率譜
power=sq.^2;
figure(3);
subplot(234);
plot(f,power);
xlabel('頻率(Hz)');
ylabel('功率譜');
title('白噪聲功率譜');
grid;
%求對數譜
ln=log(sq);
figure(3);
subplot(235);
plot(f,ln);
xlabel('頻率(Hz)');
ylabel('對數譜');
title('白噪聲對數譜');
grid;
%用IFFT恢復原始信號
xifft=ifft(y);
magx=real(xifft);
ti=[0:length(xifft)-1]/fs;
figure(3);
subplot(236);
plot(ti,magx);
xlabel('t');
ylabel('y');
title('通過IFFT轉換的白噪聲波形');
grid;