clear?all；
close?all;
fs=10;
n=100

x=randn（1n）;
figure（1）;
subplot（221）;
plot（x）;%作白噪聲的時(shí)域波形
xlabel（‘t‘）;
ylabel（‘x‘）;
title（‘高斯白噪聲時(shí)域波形‘）;
grid;
%進(jìn)行FFT變換并做頻譜圖
X=fft（x）;%進(jìn)行fft變換
magx=abs（X）;%求幅值
fx=（0:length（X）-1）‘*fs/length（X）;%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換
figure（1）;
subplot（222）;
plot（fxmagx）;%做頻譜圖
xlabel（‘頻率（Hz）‘）;
ylabel（‘幅值‘）;
title（‘高斯白噪聲幅頻譜圖‘）;
grid;
lnx=log（magx）;%求對(duì)數(shù)譜
figure（1）;
subplot（223）;
plot（fxlnx）;
xlabel（‘頻率（Hz）‘）;
ylabel（‘對(duì)數(shù)譜‘）;
title（‘高斯白噪聲對(duì)數(shù)譜‘）;
grid;???????????????????
powerx=magx.^2;?%求功率譜
figure（1）;
subplot（224）;
plot（fxpowerx）;
xlabel（‘頻率（Hz）‘）;
ylabel（‘功率譜‘）;
title（‘高斯白噪聲功率譜‘）;
grid;????????????????????

y=rand（1n）;
figure（2）;
subplot（221）;
plot（y）;%作白噪聲的時(shí)域波形
xlabel（‘t‘）;
ylabel（‘y‘）;
title（‘均勻噪聲時(shí)域波形‘）;
grid;
%進(jìn)行FFT變換并做頻譜圖
Y=fft（y）;%進(jìn)行fft變換
magy=abs（Y）;%求幅值
fy=（0:length（Y）-1）‘*fs/length（Y）;%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換
figure（2）;
subplot（222）;
plot（fymagy）;%做頻譜圖
xlabel（‘頻率（Hz）‘）;
ylabel（‘幅值‘）;
title（‘均勻噪聲幅頻譜圖‘）;
grid;
lny=log（magy）;%求對(duì)數(shù)譜
figure（2）;
subplot（223）;
plot（fylny）;
xlabel（‘頻率（Hz）‘）;
ylabel（‘對(duì)數(shù)譜‘）;
title（‘均勻噪聲對(duì)數(shù)譜‘）;
grid;???????????????????
powery=magy.^2;?%求功率譜
figure（2）;
subplot（224）;
plot（fypowery）;
xlabel（‘頻率（Hz）‘）;
ylabel（‘功率譜‘）;
title（‘均勻噪聲功率譜‘）;
grid;??????????????????


a=[0.98443604?0.83392334?0.07568359];
b=[0.99572754?0.94790649?0.53567505];
Hz?=?zpk（ab11/44100）;
Hz2=get（tf（Hz））
num?=?cell2mat（Hz2.num（））;
den?=?cell2mat（Hz2.den（1））;
[H?F]=freqs（dennum）;
z=filter（numdenx）;%濾波產(chǎn)生粉紅色噪聲
figure（3）;
subplot（221）;
plot（z）;%作噪聲的時(shí)域波形
xlabel（‘t‘）;
ylabel（‘z‘）;
title（‘粉紅噪聲時(shí)域波形‘）;
grid;
%進(jìn)行FFT變換并做頻譜圖
Z=fft（z）;%進(jìn)行fft變換
magz=abs（Z）;%求幅值
fz=（0:length（Z）-1）‘*fs/length（Z）;%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換
figure（3）;
subplot（222）;
plot（fzmagz）;%做頻譜圖
xlabel（‘頻率（Hz）‘）;
ylabel（‘幅值‘）;
title（‘粉紅噪聲幅頻譜圖‘）;
grid;
lnz=log（magz）;%求對(duì)數(shù)譜
figure（3）;
subplot（223）;
plot（fzlnz）;
xlabel（‘頻率（Hz）‘）;
ylabel（‘對(duì)數(shù)譜‘）;
title（‘粉紅噪聲對(duì)數(shù)譜‘）;
grid;???????????????????
powerz=magz.^2;?%求功率譜
figure（3）;
subplot（224）;
plot（fzpowerz）;
xlabel（‘頻率（Hz）‘）;
ylabel（‘功率譜‘）;
title（‘粉紅噪聲功率譜‘）;
grid;?

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????2372??2015-12-17?22:28??pinputu.m