随机信号分析期末考核模板.docx

2016-2017第一学期随机信号分析期末考核 姓名： 吴炜 学号： 1407030540061 班级： 14级电本2班 按随机过程X(n)=0.8X(n-1)+W（n），其中W（n）为均值为0方差为4的高斯白噪声序列。（1）模拟产生X（n）序列的500观测样本序列绘出波形；（2）用观测点估计信号的均值和方差；（3）估计该过程的自相关函数和功率谱密度，绘出波形；画出样本序列的概率直方图。分析：（1）random格式：random(x)；（注释：x是个正整数），random功能：即可随机产生0~x之间的一切整数。而且它是一个函数，不能单独使用。（2）xcorr，是指互相关函数。xcorr函数只能估计有限序列。c = xcorr(x,y) 返回矢量长度为2*N-1互相关函数序列，其中x和y的矢量长度均为N，如果x和y的长度不一样，则在短的序列后补零直到两者长度相等。(3)利用rand()函数，rand（M,N）函数，random（‘unif’，a，b，N,M）函数产生随机序列；（4）利用mean（）函数求均值，var（）函数求方差；（5）先用xcorr函数求自相关函数，再用fft函数求功率谱；（6）用函数hist（）画出概率直方图,对所产生的随机数进行希尔伯特变换：MATLAB语言中有希尔伯特变换的函数hilbert（），使用时，如果用Z=hilbert（X），给出的Z是以X做实部的解析信号，相当于Z=X+jH[X];如果只想获得X的希尔伯特变换的结果，应再使用image（Z）函数。程序： clear all; b=[1]; a=[1 -0.8]; % 由描述的差分方程，得到系统传递函数 h=impz(b,a,20); % 得到系统的单位冲激函数，在20点处已经可以认为值是0 randn(state,0); w=normrnd(0,2,1,500); % 产生题设的白噪声随机序列，标准差为2 x=filter(b,a,w); % 通过线形系统，得到输出就是题目中要求的AR 过程 plot(x,r); ylabel(X(n)); title(产生的X(n) 随机序列); grid; x_mean=mean(x) % 进行时间平均，求均值 x_var=var(x) % 求序列的方差 N=500;plot(x);hist(x,20);Mlag=20; % 定义最大自相关长度 Rx=xcorr(x,Mlag,coeff); m=-Mlag:Mlag;subplot(121); stem(m,Rx,r.);grid on;window=hamming(20); % 采用hanmming 窗，长度为20 noverlap=10; % 重叠的点数 Nfft=512; % 做FFT 的点数 Fs=1000; % 采样频率，为1000Hz [Px,f]=pwelch(x,window,noverlap,Nfft,Fs, onesided); % 估计功率谱密度 f=[-fliplr(f) (f(2:end))]; % 构造一个对称的频率，范围是[-Fs/2, Fs/2] Py=[-fliplr(Px) (Px(2:end))]; % 对称的功率谱 subplot(122); plot(f,10*log10(Py), b); grid on;仿真结果：（1）（2）x_mean = -0.6026x_var =9.6052（3）（4）已知一线性系统为最大频率为8KHZ，带宽为3KHZ的带通滤波器。（1）产生一组均匀分布的白噪声序列，让其通过上述滤波器，画出输出信号的自相关函数，功率谱密度和概率直方图；（2）产生一组高斯分布的白噪声序列，让其通过上述滤波器，画出输出信号的自相关函数，功率谱密度和概率直方图。 分析：带通滤波器（band-pass filter）是一个允许特定频段的波通过同时屏蔽其他频段的设备。比如RLC振荡回路就是一个模拟带通滤波器。random格式：random(x)；（注释：x是个正整数），random功能：即可随机产生0~x之间的一切整数。而且它是一个函数，不能单独使用。xcorr，是指互相关函数。xcorr函数只能估计有限序列。c = xcorr(x,y) 返回矢量长度为2*N-1互相关函数序列，其中x和y的矢量长度均为N，如果x和y的长度不一样，则在短的序列后补零直到两者长度相等。利用random(unif,0,1,1,N)产生一组均匀分布