离散时间随机过程建模实验报告

《离散时间随机过程建模实验报告》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、WORD文档可编辑实验报告姓名：实验名称：离散时间随机过程建模学号：课程名称：统计信号处理基础班级：实验室名称：组号：实验日期：2012.10.10一、实验目的、要求本实验的目的是在了解了Matlab编程语言的编程和调试的基础上，利用Matlab本身自带的函数来验证随机信号建模，并掌握子函数的编写方法。计算机根据理论模型生成随机数，学生需要根据观测的数据编程来计算随机过程的参数。本实验主要是为了让学生在充分理解不同的随机过程建模的理论方法的基础上，用计算机来认识理论和仿真模型之间的差异。要求包括以下几个部分：1．要求独立完成实验的内容所要求的各项功能，编制完整的M

2、atlab程序，并在程序中注释说明各段程序的功能。2．要填写完整的实验报告，报告应包含程序、图形和结论。要求记录在实验过程中碰到的问题，以及解决的方法和途径。3．实验报告是现场用Word填写并打印完成。个人或组必须在报告上署名。二、实验环境验所要求的设备：每组包含完整的计算机1台；可共用的打印机1台，A4纸张若干；计算机上安装的软件包括：Matlab6.5以上（应包含SignalProcessingToolbox,FilterDesignToolbox）；Word2000以上；三、实验原理实验内容包括2个，实验1．本实验主要是采用FIR最小二乘逆滤波器来实现反卷积

3、。假定观测的数据是由信号通过脉冲响应为技术资料专业分享WORD文档可编辑的滤波器而生成的。如果从中恢复的信号是一组脉冲序列，其中的取值为25405565859511013014015510.80.70.50.70.20.90.50.60.3a.根据上面的关系，画出观测数据，并看看是否能通过的峰值来确定的幅度和位置。（需要调用conv函数）程序n=[0:50];g=cos((n-25)/5).*exp(-(n-25).*(n-25)/100);g(51:200)=0;x=zeros(200,1);x(25)=1;x(40)=0.8;x(55)=0.7;x(65)=0

4、.5;x(85)=0.7;x(95)=0.2;x(110)=0.9;x(130)=0.5;x(140)=0.6;x(155)=0.2;y=conv(x,g);figure(1)subplot(3,1,1),plot(g);title('滤波器冲击响应');xlabel('n');ylabel('响应幅值');subplot(3,1,2),plot(x);title('输入序列x');xlabel('n');ylabel('幅值');subplot(3,1,3),plot(y);title('滤波器输出');xlabel('n');ylabel('幅值');技术资料

5、专业分享WORD文档可编辑图形b.用教材中给出的spike.m函数来设计长度的最小二乘逆滤波器，并确定最佳的延迟。程序err1=10;N=50;n1=0;forn0=0:200;[h,err]=spike(g,n0,N);iferr

6、H);n2=[-37:length(x1)-38];figure(3)subplot(2,1,1),plot(x);title('输入序列x(n)');xlabel('n');ylabel('幅值');subplot(2,1,2),plot(n2,x1);axis([0,200,0,1]);title('逆滤波器输出y');xlabel('n');ylabel('幅值');图形技术资料专业分享WORD文档可编辑d.如果观测数据中还包含噪声，即观测数据为，其中是方差为的高斯白噪声，分别取，重复b和c中的计算分析。评论这时获得的结果。0.0001时程序v=0.0001

7、;y1=y(1:205)+v.*randn(1,205);x2=conv(y1,H);n3=[-n1:length(x2)-n1-1];figure(4)subplot(2,1,1),plot(x);title('输入序列x(n)');xlabel('n');ylabel('幅值');subplot(2,1,2),plot(n3,x2);axis([0,200,0,5]);title('逆滤波器输出y1');xlabel('n');ylabel('幅值');图形0.001时程序v=0.001;y1=y(1:205)+v.*randn(1,205);x2=conv(

8、y1,H)