加噪声信号频谱图

matlab关于加均匀噪声、高斯噪声、粉红噪声的信号，分别画出其频谱图、对数谱、功率谱

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

x=randn（1n）;
figure（1）;
subplot（221）;
plot（x）;%作白噪声的时域波形
xlabel（‘t‘）;
ylabel（‘x‘）;
title（‘高斯白噪声时域波形‘）;
grid;
%进行FFT变换并做频谱图
X=fft（x）;%进行fft变换
magx=abs（X）;%求幅值
fx=（0:length（X）-1）‘*fs/length（X）;%进行对应的频率转换
figure（1）;
subplot（222）;
plot（fxmagx）;%做频谱图
xlabel（‘频率（Hz）‘）;
ylabel（‘幅值‘）;
title（‘高斯白噪声幅频谱图‘）;
grid;
lnx=log（magx）;%求对数谱
figure（1）;
subplot（223）;
plot（fxlnx）;
xlabel（‘频率（Hz）‘）;
ylabel（‘对数谱‘）;
title（‘高斯白噪声对数谱‘）;
grid;                   
powerx=magx.^2; %求功率谱
figure（1）;
subplot（224）;
plot（fxpowerx）;
xlabel（‘频率（Hz）‘）;
ylabel（‘功率谱‘）;
title（‘高斯白噪声功率谱‘）;
grid;                    

y=rand（1n）;
figure（2）;
subplot（221）;
plot（y）;%作白噪声的时域波形
xlabel（‘t‘）;
ylabel（‘y‘）;
title（‘均匀噪声时域波形‘）;
grid;
%进行FFT变换并做频谱图
Y=fft（y）;%进行fft变换
magy=abs（Y）;%求幅值
fy=（0:length（Y）-1）‘*fs/length（Y）;%进行对应的频率转换
figure（2）;
subplot（222）;
plot（fymagy）;%做频谱图
xlabel（‘频率（Hz）‘）;
ylabel（‘幅值‘）;
title（‘均匀噪声幅频谱图‘）;
grid;
lny=log（magy）;%求对数谱
figure（2）;
subplot（223）;
plot（fylny）;
xlabel（‘频率（Hz）‘）;
ylabel（‘对数谱‘）;
title（‘均匀噪声对数谱‘）;
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）;%滤波产生粉红色噪声
figure（3）;
subplot（221）;
plot（z）;%作噪声的时域波形
xlabel（‘t‘）;
ylabel（‘z‘）;
title（‘粉红噪声时域波形‘）;
grid;
%进行FFT变换并做频谱图
Z=fft（z）;%进行fft变换
magz=abs（Z）;%求幅值
fz=（0:length（Z）-1）‘*fs/length（Z）;%进行对应的频率转换
figure（3）;
subplot（222）;
plot（fzmagz）;%做频谱图
xlabel（‘频率（Hz）‘）;
ylabel（‘幅值‘）;
title（‘粉红噪声幅频谱图‘）;
grid;
lnz=log（magz）;%求对数谱
figure（3）;
subplot（223）;
plot（fzlnz）;
xlabel（‘频率（Hz）‘）;
ylabel（‘对数谱‘）;
title（‘粉红噪声对数谱‘）;
grid;                   
powerz=magz.^2; %求功率谱
figure（3）;
subplot（224）;
plot（fzpowerz）;
xlabel（‘频率（Hz）‘）;
ylabel（‘功率谱‘）;
title（‘粉红噪声功率谱‘）;
grid; 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2372  2015-12-17 22:28  pinputu.m