资源描述：
《统计与自适应信号处理作业》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、.山东大学实验报告题目统计与自适应信号处理作业信息科学与工程院（系）电子信息工程专业课程名称统计与自适应信号处理...CHAPTER4一．实验要求二．实验步骤与程序clc,clearall,closeall[y,Fs,bits]=wavread('sddx_normal.wav');%读取wav文件，并返回采样数据到向量y中，Fs表示采样频率，bits表示采样的位数Fs%显示采样频率bits%显示采样位数ft=y(:,1);sigLength=length(ft)%获取声音长度sound(y,Fs,bits);%%播放声音%绘制波形图t=(0:sigLength-1)/Fs;figure(
2、1);plot(t,ft),title('theentirespeechwaveform'),grid;xlabel('Time(s)');%横轴是时间,不是采样数...ylabel('Amplitude');%求每帧信号的绝对值%enframe分帧求出的各个帧放于矩阵的一行f1=enframe(ft,hamming(160),80)%对信号进行分帧加汉明窗,窗长度为160，重叠为80E=abs(mean(f1'));%求均值的绝对值，mean函数对矩阵的列求均值，所以求转置figure(2)plot(E)title('每帧信号的绝对值')xlabel('帧数')%绘制信号的方差VAR=v
3、ar(f1');%求方差的绝对值figure(3)plot(VAR)title('方差')xlabel('帧数')%绘制某一帧的波形图figure(4)f2=f1(55,:);%提取第28帧plot(f2)title('第55帧信号波形')xlabel('采样点数')%并不是真正的时间%绘制一帧信号的自相关figure(5)Acorr=xcorr(f2);%求这一帧信号的自相关plot(Acorr);title('这一帧图像的自相关');%绘制一帧信号的功率谱figure(6)spec=spectrum(f2)%求这一帧信号的功率谱plot(spec);title('这一帧图像的功率谱')
4、;...一．实验结果与分析运行上述程序后，结果如下图所示：上图是语音信号的波形图，有wavread可知信号的采样频率为8000Hz，信号采样点数为30001，图中横坐标是时间，单位秒（s）。上图是每帧信号的绝对值，很明显，“山东大学”四个字对应的帧的绝对值较大。...上图是每帧信号的方差，很明显，“山东大学”四字对应的帧的方差较大。上图是取自第55帧处的波形信号，横轴代表采样点。...上图是第55帧的自相关上图是第55帧的功率谱。...CHAPTER5一．实验要求二．实验步骤和程序clc,clear,closeallsigmaW=0.19;sigmaV=1.4;%w和v的方差w=sqrt(
5、sigmaW)*randn(100,1);v=sqrt(sigmaV)*randn(100,1);%假设对输入输出进行100次采样s=filter([1],[10.8],w);%s的表达式x=filter([1],[1],s)+v;%x的表达式%期望输出与实际输入的互相关rs0=1;rs1=0.8;rs2=0.8*0.8;%自相关rx0=2.5;rx1=0.8;rx2=0.8*0.8;rxd=[rs0,rs1]';%互相关Rx=[rx0,rx1;rx1,rx0];%自相关矩阵%由正则方程Rc=d求解预测系数c=Rdc=Rxrxdy=filter(c,[1],x);%绘制s和x的波形图.
6、..figure,subplot(211),plot(s),holdon,plot(x,'--r')axis([0100-33])legend('s','x')title('s和x的波形图')subplot(212)plot(s),holdonplot(y,'.-r')legend('s','y')axis([0100-33])title('s和y的波形图')三．实验结果与分析运行上述程序后，结果如下图所示：上面的一幅图中的蓝色实线是我们感兴趣的信号s(n)(取100个点)，红色的点划线是被噪声干扰的信号x(n)；为了从x(n)中恢复s(n)，设计一个二阶滤波器，得到y(n)就是对s(n)
7、的估计，下面的一幅图中蓝色实线是我们感兴趣的信号s(n)，红色的点划线是估计信号y(n)，由图可以看出，y(n)能够估计s(n)的趋势，适当增大滤波器的阶数，可是提高预测的准确度。...CHAPTER7一．实验要求二．实验步骤与程序本实验中，不妨采用前加窗的FIRLSE滤波器，程序如下clc,clearall,closeallsigmaW=0.36;sigmaV=1.5;%w和v的方差w=sqrt(sigmaW)*ran
