- 1、本文档共4页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验随机过程的模拟与特征估计
实验随机过程的模拟与特征估计
PAGE
实验随机过程的模拟与特征估计
实验 随机过程的模拟与特征估计
实验结果及分析:
实验 (1)估计x(n)=*x(n-1)+1+4.*randn(N,1)随机序列的自相关函数和功率谱
MATLAB仿真程序
%估计x(n)=*x(n-1)+1+4.*randn(N,1)随机序列的自相关函数和功率谱
%x(n)=*x(n-1)+1+4.*randn(N,1)随机序列的产生
a=;
N=500;
w=1+2.*randn(N,1);
x(1)=w(1);
for n=2:N
x(n)=a*x(n-1)+w(n);
end
subplot(3,2,1);
plot(x);
title(随机序列x(n)=*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);
periodogram(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 on
subplot(3,2,6);
pwelch(x,128,64,[],1000);
axis([0 500 -50 10]);
title(韦尔奇功率谱估计);
grid on;
实验 (2)x=sin(2*pi**n)+2*cos(2*pi**n)+randn(N,1)随机序列的自相关函数和功率谱
N=256时的结果:
N=1024时的结果:
MATLAB仿真程序
N=256:
%估计x=sin(2*pi**n)+2*cos(2*pi**n)+randn(N,1)随机序列的自相关函数和功率谱
%x=sin(2*pi**n)+2*cos(2*pi**n)+randn(N,1)随机序列的产生
N=256; %N=256或1024
w=randn(N,1);
for n=1:N
x(n)=sin(2*pi**n)+2*cos(2*pi**n)+w(n);
end
subplot(3,2,1);
plot(x);
axis([0 260 -8 8]);
title(随机序列x(N)=sin(2*pi**n)+2*cos(2*pi**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(周期图功率谱估计/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 on
subplot(3,2,6);
pwelch(x,128,64,[],1000);
axis([0 500 -50 10]);
title(韦尔奇功率谱估计/N=256);
grid on;
N=1024:
%估计x=sin(2*pi**n)+2*cos(2*pi**n)+randn(N,1)随机序列的自相关函数和功率谱
%x=sin(2*pi**n)+2*c
您可能关注的文档
最近下载
- 济南版(2024)初中生物学七年级上册《脊椎动物身体背部有脊柱》教学设计及反思.docx
- 中铁建工集团质量管理手册.pdf
- 二级中医医院评审细则解读院感部分.pptx
- 2024年迪瑞医疗分析报告:强化协同,仪器放量布局市场.pdf
- 部编版初中道德与法治九年级上册单元作业设计 (优质案例12页) .pdf
- 广东实验中学2023-2024学年八年级上学期期中考试语文试卷.docx VIP
- 2020年世界发展报告:全球价值链时代的贸易换发展.pdf VIP
- 湖北省水利工程重大设计变更报告编制大纲.pdf
- 2024高中地理教师课程标准考试模拟试卷及参考答案.docx VIP
- 《材料成型工艺学》全套教学课件.ppt
文档评论(0)