%产生色噪声%
N=2048;
x1_rand=rand(1,N);
x2_rand=rand(1,N);
xn=0:N-1;
x=x1_rand-x2_rand;
figure(1)
plot(xn,x,'r')
title('色噪声')
%******统计色噪声*******************
E=mean(x) ; %色噪声的均值E = 0.0125
S=var(x) ; %色噪声的方差S = 0.1669
X=E*E+S ; %色噪声的均方值X = 0.1670
%**************求色噪声概率密度函数*****************%
eachi=linspace(min(x),max(x),42);
yyi=hist(x,eachi); %计算各个区间的个数
yyi=yyi/length(x); %对各个区间的个数归一化处理
figure(2); %绘制色噪声的概率密度函数
plot(eachi,yyi,'-b')
title('色噪声的概率密度函数')
%***************色噪声自相关函数 *****%
maxlags=2*N-1;
Rx=xcorr(x,x,maxlags); %色噪声的自相关函数Rx
tau=-maxlags:1:maxlags;
figure(3);
plot(tau,Rx,'-r')
title('色噪声的自相关函数'); %色噪声的自相关函数波形
本内容试读结束,登录后可阅读更多
下载后可阅读完整内容,剩余2页未读,立即下载