1.(5″)证明下面定理：任何一个无偏估计子方差的下界叫作cramer-rao下界。
定理：令为一样本向量，是的条件密度，若是的一个无偏估计子，且存在，则。
式中。其中是的某个不包含的正函数。
2.(10″)wiener滤波是信号处理中最常用和基础的波形估计工具之一，对其在自己研究领域的应用情况进行一个简单综述。
3.(5″)二阶滑动平均过程由。
定义，式中表示正态分布，其均值为零、方差为。求x(n)的功率谱。
4.(20″)信号的函数表达式为：
其中，a(t)为一随时间变化的随机过程，dn(t)为经过390-410hz带通滤波器后的高斯白噪声，n(t)为高斯白噪声，采样频率为1khz，采样时间为2.048s。
1) 利用现代信号处理知识进行信号的谱估计；
2) 利用现代信号处理知识进行信号的频率提取；
3) 分别利用wiener滤波和kalman滤波进行去噪；
4) 利用wigner-ville分布分析信号的时频特征。
5.(10″)附件中表sheet1 为某地2023年4月28日凌晨12点至2023年5月4日凌晨12点的电力系统负荷数据，采样时间间隔为1小时，利用arma方法**该地5月5日的电力系统负荷，并给出**误差（5月5日的实际负荷数据如表sheet2）。
定理：令为一样本向量，是的条件密度，若参数估计是真实参数的一个无偏估计子，且、存在，则的均方误差所能达到的下界（称为cramer-rao下界）等于fisher信息的导数，即：
不等式中等号成立的充分必要条件是：
其中是的某个正函数，与样本无关。
证明：由假设条件知，或，因此有。
对上式两边求关于的偏导，得。即有。
另一方面，由复合函数的求导法，又有。
由于是的条件概率密度，故。
将式（1-5）和式（1-6）带入式（1-4），得。或改写作。
由cauchy-schwartz不等式知，对于任意两个复函数和，恒有不定式：
成立，并且当且仅当，等号成立，将cauchy-schwartz不等式应用于式（1-7），则有。或等价为。
由cauchy-schwartz不等式等号成立的条件知：当且仅当，即式（1-2）成立时，不等式（1-8）才取等号。
注意到，故有。
另由公式知。
将式（1-9）和式（1-10）代入式（1-8），直接得到不等式（1-1）。根据前面的分析，不等式等号成立的充分必要条件是式（1-2）成立。
解：对取延迟形式：
于是有：展开上式得到：
即：对上式进行傅里叶变换，则有：
从而得到：即。
music算法原理：
music(multiple signal classification)算法是针对多元天线阵测向问题提出的。假定m元的均匀线阵，阵元间距为d，信号的工作波长为λ。空间信号源共有d个，各信号不相关，各阵元的噪声互不相关，噪声和信号也不相关。
因此，第m个阵元的输出为。
式中，为第ｋ个信号源的方向。将式（4-1）写成矩阵形式：
式中：、。求各阵元输出的相关矩阵，有：
式中：——噪声的方差。
对式（4-3）的相关矩阵r作特征分解，其各特征值及其相对应的特征向量分别为：
据式式（4-3），可得以下结论：
1）r的最小特征值等于，重数为(md)，即λｄ+1=…=据此，空间信号源的个数d可由下式得出：d=m-（r最小特征值的重数）。最小重数为1，因此，m阵元可测向的信号源数目的最大值为。
2）各特征向量相互正交。这些向量为矩阵ｒ列空间的基，由于最小特征值为噪声的贡献，因此与最小特征值对应的那些特征向量所张成的子空间也是噪声的贡献，称之噪声子空间，记为。这样ｒ的列空间被划分成两个子空间，即信号子空间和噪声子空间：
由于各特征向量相互正交，故有：。在信号源所在方向上，诸方向向量，均处于信号子空间中，故。构造矩阵：
显然有。music算法就是根据式（4-9）来求空间谱，有。
谱峰所对应θ值就是信号源方向的估值。
维纳滤波算法原理：
维纳(wiener)是用来解决从噪声中提取信号的一种过滤（或滤波）方法。这种线性滤波问题，可以看做是一种估计问题或一种线性估计问题。一个线性系统，如果它的单位样本响应为，当输入一个随机信号，且：
其中：表示信号，表示噪声，则输出为：
我们希望通过线性系统后得到的尽量接近于，因此称为的估计值，用表示，即：
则维纳滤波器的输入—输出关系可用下面图1表示。
图4-1 维纳滤波器的输入—输出关系。
实际上，式（4-2）所示的卷积形式可以理解为从当前和过去的观察值，来估计信号的当前值。因此，用进行过滤问题实际上是一种统计估计问题。
一般地，从当前的和过去的观察值，估计当前的信号值成为过滤或滤波；从过去的观察值，估计当前的或者将来的信号值称为外推或**；从过去的观察值，估计过去的信号值称为平滑或内插。因此维纳滤波器又常常被称为最佳线性过滤与**或线性最优估计。这里所谓的最佳与最优是以最小均方误差为准。
如果我们分别以与表示信号的真实值与估计值，而用表示他们之间的误差，即：
显然可能是正值，也可能是负值，并且它是一个随机变量。因此，用它的均方误差来表达误差是合理的，所谓均方误差最小即它的平方的统计期望最小。
卡尔曼滤波算法原理：
卡尔曼滤波是基于状态空间方法的一套递推滤波算法，在状态空间方法中，引入了状态变量的概念。卡尔曼滤波的模型包括状态空间模型和观测模型。状态模型是反映状态变化规律的模型，通过状态方程来描写相邻时刻的状态转移变化规律；观测模型反映了实际观测量与状态变量之间的关系。
kalman滤波问题就是联合观测信息及状态转移规律来得到系统状态的最优估计。
假设动态系统的状态空间模型为。
其中，x(t)为系统在时刻t的状态；y(t)为对状态的观测值；w(t)为系统噪声，方差阵为q；v(t)为观测噪声，方差阵为r；为状态转移矩阵；h为观测矩阵；为系统噪声驱动矩阵。
卡尔曼滤波的计算流程为：
计算状态估计值：
计算状态一步**：
计算新息：
计算卡尔曼滤波增益：
计算一步**均方误差：
计算一步**估计均方误差：
下面给出卡尔曼滤波的系统模型框图：
matlab图形：
图1 原始型号时域波形图。
图2 music方法对信号进行谱估计。
图3 信号频率提取。
图4 两种滤波方式时域结果图。
图5 wiener滤波去噪信号频域波形。
图6 kalman滤波去噪信号频域波形。
图7 信号wigner-ville分布。
图8 信号伪wigner-ville分布。
matlab程序：
clc; clear all
ff1=100;
ff2=300;
ff3=200;
fs = 1000采样频率。
t = 0:1/fs:2.047产生2048个采样时间点。
n=length(t);
at=normrnd(0,1,1,2048); 产生均值为0，方差为1的正态随机过程。
dnt1=wgn(1,2048,5产生高斯白噪声。
巴特沃斯带通滤波。
fs1=1000;
wp = 390 410]/(fs1/2归一化通带。
ws = 350 450]/(fs1/2);
rp = 3通带内的最大衰减。
rs = 100阻带内的最小衰减。
n,wn] =buttord(wp,ws,rp,rs); 计算阶数和截至频率。
b,a] =butter(n,wn计算h(z)
freqz(b,a,128,fs1画的幅频相频图。
dnt=filter(b,a,dnt1将产生的高斯白噪声进行巴特沃斯带通滤波器滤波。
现代信号处理作业
中国矿业大学。20 级硕士研究生课程考试试卷。考试科目现代信号处理基础 考试时间。学生姓名。学生学号。所在院系信电学院 任课教师。中国矿业大学研究生院培养管理处印制。一 自适应滤波。1 自适应df的工作原理 自适应滤波器是以最小均方误差为准则的最佳滤波器，它能自动调节其本身的单位脉冲响应h n 特性...
现代信号处理作业
中北大学。研究生 现代信号处理 作业。所在院 系 信息与通信工程学院 专业 测试计量技术及仪器 姓名王栋。学号 s20100329 班级 y100503 156 小波变换与matlab相关程序的学习 小波变换是在傅里叶变换的基础上发展起来的，它优于傅里叶分析的地方是它在时域与空域都是局部化的，其局部...
现代信号处理作业
1.总结学过的滤波器设计方法，用matlab 例子分析不同设计方法的滤波器的性能及适应场合。答 1.1模拟低通滤波器的设计方法。1.1.1 butterworth滤波器设计步骤 确定阶次n 1 已知 c s和as求butterworth df阶数n 由 求出n 2 已知 c s和 p 的衰减ap 求...