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

实验报告：实验名称：离散时间随机过程建模学号：课程名称：统计信号处理基础班级：实验室名称：组号：实验日期：2012.10.10一、实验目的、要求本实验的目的是在了解了Matlab编程语言的编程和调试的基础上，利用Matlab本身自带的函数来验证随机信号建模，并掌握子函数的编写方法计算机根据理论模型生成随机数，学生需要根据观测的数据编程来计算随机过程的参数本实验主要是为了让学生在充分理解不同的随机过程建模的理论方法的基础上，用计算机来认识理论和仿真模型之间的差异要求包括以下几个部分：1．要求独立完成实验的内容所要求的各项功能，编制完整的Matlab程序，并在程序中注释说明各段程序的功能2．要填写完整的实验报告，报告应包含程序、图形和结论要求记录在实验过程中碰到的问题，以及解决的方法和途径3．实验报告是现场用Word填写并打印完成个人或组必须在报告上署名二、实验环境验所要求的设备： 每组包含完整的计算机 1 台；可共用的打印机1台，A4纸张若干； 计算机上安装的软件包括： Matlab 6.5以上（应包含Signal Processing Toolbox, Filter Design Toolbox）； Word 2000以上；三、实验原理实验内容包括2个，实验1．本实验主要是采用FIR最小二乘逆滤波器来实现反卷积。

假定观测的数据是由信号通过脉冲响应为的滤波器而生成的如果从中恢复的信号是一组脉冲序列，其中的取值为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.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('幅值');图形b. 用教材中给出的spike.m函数来设计长度的最小二乘逆滤波器，并确定最佳的延迟。

程序err1=10;N=50;n1=0;for n0=0:200; [h,err]=spike(g,n0,N); if err

评论这时获得的结果程序r=0.001/12.*rand(1,length(g));g1=g+r;y2=conv(x,g1);H1=[];err0=1;for n0=0:200; [h,err]=spike(g1,n0,N); if err

程序A=[1,-1.978,2.853,-1.877,0.904];B=[1,-0.9,0.18];v=randn(1,100);x=filter(B,A,v);p=length(A)-1;q=length(B)-1;rx=xcorr(x,p+q,'biased');pw=abs(fft(rx))/100;figure(1)plot(pw);title('理论功率谱');xlabel('n');ylabel('幅值');图形c. 用a中编制的程序根据观测数据来求解模型的参数，把计算结果与理论模型的系数相比，有什么结论重复10次不同的样本实现，并计算10次的模型参数再取平均，与理想的系数相比，平均是否有效果程序sa=zeros(5,1);sb=zeros(1,3);for i=1:10v=randn(1,100);x=filter(B,A,v);rx=xcorr(x,length(A)+length(B)-2,'biased');[b ,a]=arma1(rx,p,q);sa=sa+a;sb=sb+b;endb1=sb./10a1=sa./10err1=sum((b-B).*(b-B))+sum((a-A').*(a-A'))err10=sum((b1-B).*(b1-B))+sum((a1-A').*(a1-A'))图形A= 1 -1.978 2.853 -1.877 0.904B= 1 -0.9 0.18b = 1.0000 -0.5106 0.3082a = 1.0000 -1.8961 2.7162 -1.7549 0.8683b1 = 1.0000 -0.9898 0.4565a1 = 1.0000 -2.2047 3.2952 -2.3071 1.1074平均效果不佳指导教师评语:成绩： 指导教师签名： 批阅日期：. v。