西南交大现代信号处理部分答案.doc

题1：(1) 是随变化的随机信号，因此=. 所以谐波信号的均值为 = 由于谐波信号的均值等于零，故其方差等于二阶矩，既有 所以x(t)的方差为 谐波信号的自相关函数 又 所以 由于x(t)的均值为0，故所以 (2) y(t)是随B变化的随机信号，因此 B是标准高斯随机变量，所以，所以. 由于统计独立，故有 而x(t)和y(t)的均值均为0，所以 题2：令，由于是零均值、方差为的高斯随机过程，和是确定的过程，所以x(n)也是一高斯随机过程，其均值是时间的函数.所以x(n)的概率密度函数是 = 在多个未知参数的情况下，Cramer-Rao不等式变为矩阵不等式 : 其中无偏估计子的协方差矩阵，而是Fisher信息矩阵J的逆矩阵，而信息矩阵J的构成元素为 本题中，计算得 = = = 求逆矩阵可得 所以可得 所以估计子的估计方差的Cramer-Rao下界分别是和 题3.信号的函数表达式为：， 其中，A(t)为一随时间变化的随机过程，为经过390-410Hz带通滤波器后的高斯白噪声，为高斯白噪声，采样频率为1kHz，采样时间为2.048s。分别利用Wiener滤波和Kalman滤波进行去噪。 解：维纳滤波原理： 维纳(Wiener)是用来解决从噪声中提取信号的一种过滤（或滤波）方法。这种线性滤波问题，可以看做是一种估计问题或一种线性估计问题。一个线性系统，如果它的单位样本响应为，当输入一个随机信号，且： 其中：表示信号，表示噪声，则输出为： 我们希望通过线性系统后得到的尽量接近于，因此称为的估计值，用表示，即： 则维纳滤波器的输入—输出关系可用下面图1表示。 图3-1 维纳滤波器的输入—输出关系 实际上,上式所示的卷积形式可以理解为从当前和过去的观察值，来估计信号的当前值。因此，用进行过滤问题实际上是一种统计估计问题。 一般地，从当前的和过去的观察值，估计当前的信号值成为过滤或滤波；从过去的观察值，估计当前的或者将来的信号值称为外推或预测；从过去的观察值，估计过去的信号值称为平滑或内插。因此维纳滤波器又常常被称为最佳线性过滤与预测或线性最优估计。这里所谓的最佳与最优是以最小均方误差为准。 如果我们分别以与表示信号的真实值与估计值，而用表示他们之间的误差，即： 显然可能是正值，也可能是负值，并且它是一个随机变量。因此，用它的均方误差来表达误差是合理的，所谓均方误差最小即它的平方的统计期望最小。 卡尔曼滤波算法原理： 卡尔曼滤波是基于状态空间方法的一套递推滤波算法，在状态空间方法中，引入了状态变量的概念。卡尔曼滤波的模型包括状态空间模型和观测模型。状态模型是反映状态变化规律的模型，通过状态方程来描写相邻时刻的状态转移变化规律；观测模型反映了实际观测量与状态变量之间的关系。Kalman滤波问题就是联合观测信息及状态转移规律来得到系统状态的最优估计。 假设动态系统的状态空间模型为 其中，X(t)为系统在时刻t的状态；Y(t)为对状态的观测值；W(t)为系统噪声，方差阵为Q；V(t)为观测噪声，方差阵为R；为状态转移矩阵；H为观测矩阵；为系统噪声驱动矩阵。 卡尔曼滤波的计算流程为： 计算状态估计值： 计算状态一步预测： 计算新息： 计算卡尔曼滤波增益： 计算一步预测均方误差： 计算一步预测估计均方误差： 下面给出卡尔曼滤波的系统模型框图： 根据上述要求，编写程序如下： fs = 1000; %采样频率 t = 0:1/fs:2.047; N = length(t); %采样点数 randn(state,0); W1=randn(1,length(t)); %高斯白噪声 %%%%%%% 一维数字滤波 %%%%%%% f2=390;f3=410; %带通滤波带宽设置 wc1=2*f2/fs; wc2=2*f3/fs; f1=[0 wc1-0.05 wc1 wc2 wc2+0.05 1]; %归一化频率f1 A=[0 0 1 1 0 0]; %设置带通1或带阻0 weigh=[1 1 1 ]; %权重 b=remez(50,f1,A,weigh); %传函分子 W2=filter(b,1,W1); %一维数字滤波 B=normrnd(0,1,1,N); %随机过程A（t） x0= sin(2*pi*t*100) +1.5*sin(2*pi*t*300)+B.*sin(2*pi*t*200); x1 = sin(2*pi*t*100) +1.5*sin(2*pi*t*300)+B.*sin(2*pi*t*200)+W1+W2; % 原始信号 figure(1) su