下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、精选优质文档-倾情为你奉上精选优质文档-倾情为你奉上专心-专注-专业专心-专注-专业精选优质文档-倾情为你奉上专心-专注-专业实验2.1 随机过程的模拟与特征估计实验结果及分析:实验2.1 (1)估计x(n)=0.8*x(n-1)+1+4.*randn(N,1)随机序列的自相关函数和功率谱MATLAB仿真程序%估计x(n)=0.8*x(n-1)+1+4.*randn(N,1)随机序列的自相关函数和功率谱%x(n)=0.8*x(n-1)+1+4.*randn(N,1)随机序列的产生a=0.8;N=500;w=1+2.*randn(N,1);x(1)=w(1); for n=2:Nx(n)=a*x
2、(n-1)+w(n); endsubplot(3,2,1);plot(x);title(随机序列x(n)=0.8*x(n-1)+1+4.*randn(N,1);grid on%估计自相关函数R=xcorr(x,coeff);subplot(3,2,2);axis(0 500 0 1);plot(R);title(自相关函数);grid on%估计功率谱%周期图功率谱估计subplot(3,2,3);periodogram(x,512,1000);axis(0 500 -50 0);title(周期图功率谱估计)%加汉宁窗window=hann(500);subplot(3,2,4);perio
3、dogram(x,window,512,1000);axis(0 500 -50 10);title(汉宁周期功率谱估计)%相关函数法R=xcorr(x)/15000;Pw=fft(R);subplot(3,2,5);f=(0:length(Pw)-1)*1000/length(Pw);plot(f,10*log10(abs(Pw);axis(0 500 -50 10);title(BT功率谱估计)grid onsubplot(3,2,6);pwelch(x,128,64,1000);axis(0 500 -50 10);title(韦尔奇功率谱估计);grid on;实验2.1 (2)x=s
4、in(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+randn(N,1)随机序列的自相关函数和功率谱N=256时的结果:N=1024时的结果:MATLAB仿真程序N=256:%估计x=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+randn(N,1)随机序列的自相关函数和功率谱%x=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+randn(N,1)随机序列的产生N=256; %N=256或1024w=randn(N,1); for n=1:N x(n)=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+
5、w(n); endsubplot(3,2,1);plot(x);axis(0 260 -8 8);title(随机序列x(N)=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+randn(N,1)/N=256);grid on%估计自相关函数R=xcorr(x,coeff);subplot(3,2,2);plot(R);axis(0 500 -1 1);title(自相关函数/N=256);grid on%估计功率谱%周期图功率谱估计subplot(3,2,3);periodogram(x,512,1000);axis(0 500 -50 0);title(周期图功率谱
6、估计/N=256)%加汉宁窗window=hann(256);subplot(3,2,4);periodogram(x,window,256,1000);axis(0 500 -50 10);title(汉宁周期功率谱估计)%相关函数法R=xcorr(x)/15000;Pw=fft(R);subplot(3,2,5);f=(0:length(Pw)-1)*1000/length(Pw);plot(f,10*log10(abs(Pw);axis(0 500 -50 10);title(BT功率谱估计/N=256)grid onsubplot(3,2,6);pwelch(x,128,64,1000
7、);axis(0 500 -50 10);title(韦尔奇功率谱估计/N=256);grid on;N=1024:%估计x=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+randn(N,1)随机序列的自相关函数和功率谱%x=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+randn(N,1)随机序列的产生N=1024; %N=256或1024w=randn(N,1); for n=1:N x(n)=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+w(n); endsubplot(3,2,1);plot(x);axis(
8、0 1030 -8 8);title(随机序列x(N)=sin(2*pi*0.05*n)+2*cos(2*pi*0.12*n)+randn(N,1)/N=1024);grid on%估计自相关函数R=xcorr(x,coeff);subplot(3,2,2);plot(R);axis(0 2000 -1 1);title(自相关函数/N=1024);grid on%估计功率谱%周期图功率谱估计subplot(3,2,3);periodogram(x,1024,1000);axis(0 500 -50 0);title(周期图功率谱估计/N=1024)%加汉宁窗window=hann(1024);subplot(3,2,4);periodogram(x,window,1024,1000);axis(0 500 -50 10);title(汉宁周期功率谱估计)%相关函数法R=xcorr(x)/15000;Pw=fft(R);subplot(3,2,5);f=(0:length(Pw)-1)*1000/length(Pw);plot(f,10*log10(abs(Pw
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 飞燕舞蹈艺术中心家长问卷调查表
- 中外合资企业合同书书模板
- 果园农业慈善活动承包合同
- 建设工程施工合同书好运
- 市场调研教案
- 冀少版三年级下册 第四单元 Do-Re-Mi教案
- 客车运输合同书
- 山东省巨野县麒麟镇第一中学八年级政治下册 12.1 感受大自然之美教案 鲁人版六三制
- 喷泉水池保洁合同
- 物业与业主合同书书物业管理合同书书
- 2024年西安中小学教师招聘真题
- 第7课 实践出真知-【中职专用】2024年中职思想政治《哲学与人生》金牌课件(高教版2023·基础模块)
- 《电工电子技术基础》高职全套教学课件
- 1.2 规划初中生活(课件)2024-2025学年七年级道德与法治上册 (统编版2024)
- 秋季儿童传染病预防知识讲座课件
- 2024年风力发电运维值班员(技师)技能鉴定考试题库-上(选择题)
- 临床研究项目申请书
- 2024-2030年中国低温储运罐市场运行前景分析与发展趋势预判研究报告
- 2024测绘成果保密协议
- 2023下半年广东深圳市大鹏新区大鹏办事处招聘社区专职工作者拟聘用人员(第二批)笔试历年典型考题及考点剖析附答案带详解
- 2024男女双方自愿离婚协议书
评论
0/150
提交评论