- 1、本文档共6页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
随机信号 产生独立的高斯随机序列 功率谱分析
随机信号分析基础
编程作业
学院:自动化与信息工程学院
班级:通信092
姓名:
学号:
题一:
编写一个产生独立的高斯随机序列的上机运行程序,用100个点数据估计它的均值,方差并与理论设定值相比较。
解:
程序如下:
%产生独立的高斯随机序列%
N=100;
S=randn(1,N);
plot(S)
a=mean(S) %求高斯随机序列的均值%
y=var(S) %求高斯随机序列的方差%
运行结果为:
a =
-0.0109
y =
0.9775
随机序列的图为:
题目二:
origin=sin(2*pi*0.01*(1:1000)+pi)+0.002*(1:1000)+randn(1,1000);
数据中含有趋势项和周期项,
1、画出原始曲线;
2、对原数据做功率谱分析;
3、拟合趋势项,对去掉趋势项后的数据做谱分析;
4、用自相关函数确定周期信号的周期和幅度,用互相关函数确定初相,给出每一步处理结果的图形。对去掉周期分量后的数据做谱分析。
解:
编写程序如下:
%画原始曲线
N=1000;
t=1:1000;
origin=sin(2*pi*0.01*t+pi)+0.002*t+randn(1,1000);
subplot(3,2,1);
plot(t,origin)
title(原始曲线)
xlabel(时间)
ylabel(幅值)
grid on
%用周期图法对原始数据进行功率谱分析
Fs=1000;
Px=abs(fft(origin,1024)).^2/length(origin);
subplot(3,2,2);
plot((0:1023)/1024*Fs,10*log10(Px))
title(原始曲线的周期图谱估值法)
xlabel(频率(Hz))
ylabel(功率谱(dB))
grid on
%用分段周期图法对原始数据进行功率谱分析
w=hanning(256);
Pxx=(abs(fft(w.*origin(1:256))).^2+abs(fft(w.*origin(129:384))).^2+abs(fft(w.*origin(257:512))).^2+abs(fft(w.*origin(385:640))).^2+abs(fft(w.*origin(513:768))).^2+abs(fft(w.*origin(641:896))).^2)/(norm(w)^2*6);
subplot(3,2,3);
plot((0:255)/256*Fs,10*log10(Pxx))
title(原始曲线的Welch谱估值法(K=6,L=256))
xlabel(频率(Hz))
ylabel(功率谱(dB))
grid on
nh=polyfit(t,origin,1) %拟合曲线
%画出去掉趋势项的曲线图
oriqu=origin-(nh(1)*t+nh(2));
subplot(3,2,4);
plot(t,oriqu)
title(去掉趋势项的信号曲线)
xlabel(时间)
ylabel(幅值)
grid on
%用周期图法对去掉趋势项数据进行功率谱分析
Pxqu=abs(fft(oriqu,1024)).^2/length(oriqu);
subplot(3,2,5);
plot((0:1023)/1024*Fs,10*log10(Pxqu))
title(去掉趋势项曲线的周期图谱估值法)
xlabel(频率(Hz))
ylabel(功率谱(dB))
axis([0,1000,-30,30])
grid on
%用分段周期图法对去掉趋势项数据进行功率谱分析
Pxxqu=(abs(fft(w.*oriqu(1:256))).^2+abs(fft(w.*oriqu(129:384))).^2+abs(fft(w.*oriqu(257:512))).^2+abs(fft(w.*oriqu(385:640))).^2+abs(fft(w.*oriqu(513:768))).^2+abs(fft(w.*oriqu(641:896))).^2)/(norm(w)^2*6);
subplot(3,2,6);
plot((0:255)/256*Fs,10*log10(Pxxqu))
title(去掉趋势项曲线的Welch谱估值法(K=6,L=256))
xlabel(频率(Hz))
ylabel(功率谱(dB))
grid on
figure
%周期信号的自相关函数求解与分析
orizhou=sin(2*pi*0.01*t+pi);
corizhou=xcorr(orizhou,200,unbias
您可能关注的文档
- 安全标识管理标准..doc
- 时间序列模型2.doc
- 安全检查人员试题..doc
- 离散时间序列(图).docx
- 数字信号处理实验 DFT分析连续信号频谱.doc
- 安全标志及其使用导则2008(必威体育精装版完整版)..doc
- 土方开挖施工手册.doc
- 安全检测技术..doc
- 安全流行语..doc
- 安全检查评分表(JGJ59-2011)..doc
- 第十一章 电流和电路专题特训二 实物图与电路图的互画 教学设计 2024-2025学年鲁科版物理九年级上册.docx
- 人教版七年级上册信息技术6.3加工音频素材 教学设计.docx
- 5.1自然地理环境的整体性 说课教案 (1).docx
- 4.1 夯实法治基础 教学设计-2023-2024学年统编版九年级道德与法治上册.docx
- 3.1 光的色彩 颜色 电子教案 2023-2024学年苏科版为了八年级上学期.docx
- 小学体育与健康 四年级下册健康教育 教案.docx
- 2024-2025学年初中数学九年级下册北京课改版(2024)教学设计合集.docx
- 2024-2025学年初中科学七年级下册浙教版(2024)教学设计合集.docx
- 2024-2025学年小学信息技术(信息科技)六年级下册浙摄影版(2013)教学设计合集.docx
- 2024-2025学年小学美术二年级下册人美版(常锐伦、欧京海)教学设计合集.docx
最近下载
- 17J008 挡土墙(重力式、衡重式、悬臂式)(必威体育精装版).pdf
- 造血干细胞移植的护理干预.pptx
- 布料车岗位安全规程.pptx
- YDT 5178-2017 通信管道人孔和手孔图集.docx VIP
- 精品解析:【区级联考】上海徐汇区2019届九年级学习能力诊断(二模)数学试题(解析版).pdf VIP
- 精品解析:广东省佛山市南海区,三水区2022-2023学年九年级上学期数学期末考试(原卷版).pdf VIP
- 一种护筒导向架结构.pdf VIP
- 老旧小区雨污分流改造要点与难点分析.docx VIP
- 鞍钢宪法及后福特主义.pdf
- 精品解析:广东省广州市2022-2023学年九年级上学期期末数学考前模拟试题(三)(解析版).pdf VIP
文档评论(0)