




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
大连理工大学实验预习报告学院〔系〕:信息与通信工程学院专业:电子信息工程班级:电子1401姓名:******学号:******实验时间:实验室:c221指导教师:郭成安实验II:系统对随机信号响应的统计特性分析、功率谱分析及应用实验实验目的和要求掌握直接法估计随机信号功率谱的原理和实现方法;掌握间接法估计随机信号功率谱的原理和实现方法;掌握系统对随机信号响应的统计特性分析及仿真实现方法。熟悉MATLAB信号处理软件包的使用。实验原理和内容〔一〕实验原理:1.直接法估计随机信号功率谱原理直接法又称为周期图法,它是把随机信号x(n)的N点观察数据xN(n)视为一能量有限信号,直接取xN(n)的傅里叶变换,得到XN(ejω),然后取其模值的平方,并除以N,作为对x(n)真实的功率谱P(ejω)的估计。工程上,常使用离散Fourier变换〔DFT,编程上使用其快速算法FFT〕,即PX(k)=进行计算。2.间接法估计随机信号功率谱间接法的理论根底是Wiener-Khintchine定理,具体的实现方法是先由xN(n)估计出自相关函数,然后对求傅里叶变换得到xN(n)的功率谱,记之为XN(ejω),并以此作为对真实功率谱P(ejω)的估计。工程上,常使用离散Fourier变换〔DFT,编程上使用其快速算法FFT〕,即,,进行计算。因为由这种方法求出的功率谱是通过自相关函数间接得到的,所以又称为间接法或Blackman-Tuckey(BT)法,该方法是FFT出现之前常用的谱估计方法。3.时域中系统对随机信号响应的统计特性分析及仿真根据系统卷积性质,计算系统输出信号的统计特性。有如下性质:频域中系统对随机信号响应的统计特性分析及仿真根据卷积定理,输入、输出信号功率谱的关系为RY(ejω)=RX(ejω)|H(ejω)|2。在计算系统输出信号功率谱时,如果在时域时计算困难,可以按照上式在频域计算。〔二〕实验内容:直接法估计随机信号功率谱(1)生成1024点数据的随机信号其中f1=30Hz,f2=100Hz,为在[0,2п]内的均匀分布的随机变量,N(n)是数学期望为0,方差为1的高斯白噪声。(2)用周期图法计算的功率谱,并绘图。(3)用MATLAB函数periodogram重新计算的功率谱,并与(2)做比拟。2.间接法估计随机信号功率谱(1)计算以上的自相关函数。(2)通过计算自相关函数的Fourier变换,求的功率谱并绘图。(3)利用MATLAB函数psd、pwelch重新计算的功率谱,并与(2)做比拟。3.系统对随机信号响应的统计特性分析及仿真(1)生成含500点数据的高斯分布白噪声随机信号。(2)设计一个带通系统,其上、下截止频率分别为4KHz和3KHz。(3)计算通过以上带通滤波器的自相关函数和功率谱密度。实验步骤(1)生成1024点随机信号X〔n〕,高斯白噪声N〔n〕。(2)周期法绘制其功率谱。(3)利用MATLAB函数periodogram计算功率谱,与〔2〕比拟。(4)计算其自相关函数。(5)计算自相关函数的傅里叶变换。(6)求其功率谱。(7)利用MATLAB函数psd、pwelch重新计算的功率谱,并与(2)做比拟。(8)生成含500点数据的高斯分布白噪声随机信号。(9)设计一个带通系统,其上、下截止频率分别为4KHz和3KHz。(10)计算通过以上带通滤波器的自相关函数和功率谱密度。大连理工大学实验报告学院〔系〕:信息与通信工程学院专业:电子信息工程班级:1401姓名:******学号:******实验时间:2023/11/04实验室:C221指导教师:郭成安实验II:系统对随机信号响应的统计特性分析、功率谱分析及应用实验实验目的和要求掌握直接法估计随机信号功率谱的原理和实现方法;掌握间接法估计随机信号功率谱的原理和实现方法;掌握系统对随机信号响应的统计特性分析及仿真实现方法。熟悉MATLAB信号处理软件包的使用。实验原理和内容〔一〕实验原理:1.直接法估计随机信号功率谱原理直接法又称为周期图法,它是把随机信号x(n)的N点观察数据xN(n)视为一能量有限信号,直接取xN(n)的傅里叶变换,得到XN(ejω),然后取其模值的平方,并除以N,作为对x(n)真实的功率谱P(ejω)的估计。工程上,常使用离散Fourier变换〔DFT,编程上使用其快速算法FFT〕,即PX(k)=进行计算。2.间接法估计随机信号功率谱间接法的理论根底是Wiener-Khintchine定理,具体的实现方法是先由xN(n)估计出自相关函数,然后对求傅里叶变换得到xN(n)的功率谱,记之为XN(ejω),并以此作为对真实功率谱P(ejω)的估计。工程上,常使用离散Fourier变换〔DFT,编程上使用其快速算法FFT〕,即,,进行计算。因为由这种方法求出的功率谱是通过自相关函数间接得到的,所以又称为间接法或Blackman-Tuckey(BT)法,该方法是FFT出现之前常用的谱估计方法。3.时域中系统对随机信号响应的统计特性分析及仿真根据系统卷积性质,计算系统输出信号的统计特性。有如下性质:频域中系统对随机信号响应的统计特性分析及仿真根据卷积定理,输入、输出信号功率谱的关系为RY(ejω)=RX(ejω)|H(ejω)|2。在计算系统输出信号功率谱时,如果在时域时计算困难,可以按照上式在频域计算。〔二〕实验内容:直接法估计随机信号功率谱(1)生成1024点数据的随机信号其中f1=30Hz,f2=100Hz,为在[0,2п]内的均匀分布的随机变量,N(n)是数学期望为0,方差为1的高斯白噪声。(2)用周期图法计算的功率谱,并绘图。(3)用MATLAB函数periodogram重新计算的功率谱,并与(2)做比拟。2.间接法估计随机信号功率谱(1)计算以上的自相关函数。(2)通过计算自相关函数的Fourier变换,求的功率谱并绘图。(3)利用MATLAB函数psd、pwelch重新计算的功率谱,并与(2)做比拟。3.系统对随机信号响应的统计特性分析及仿真(1)生成含500点数据的高斯分布白噪声随机信号。(2)设计一个带通系统,其上、下截止频率分别为4KHz和3KHz。(3)计算通过以上带通滤波器的自相关函数和功率谱密度。实验步骤(1)生成1024点随机信号X〔n〕,高斯白噪声N〔n〕。(2)周期法绘制其功率谱。(3)利用MATLAB函数periodogram计算功率谱,与〔2〕比拟。(4)计算其自相关函数。(5)计算自相关函数的傅里叶变换。(6)求其功率谱。(7)利用MATLAB函数psd、pwelch重新计算的功率谱,并与(2)做比拟。(8)生成含500点数据的高斯分布白噪声随机信号。(9)设计一个带通系统,其上、下截止频率分别为4KHz和3KHz。(10)计算通过以上带通滤波器的自相关函数和功率谱密度。主要仪器设备微型计算机、Matlab开发环境〔本报告采用MATLAB2023a〕五、编程代码与操作方法实验内容〔1〕:直接法估计随机信号功率谱(1)生成1024点数据的随机信号其中f1=30Hz,f2=100Hz,为在[0,2п]内的均匀分布的随机变量,N(n)是数学期望为0,方差为1的高斯白噪声。(2)用周期图法计算的功率谱,并绘图。(3)用MATLAB函数periodogram重新计算的功率谱,并与(2)做比拟。代码如下:clc;clear;%清屏及去除缓存N=1024;f1=30;f2=100;fs=1000;t=(0:N-1)/fs;%第一问Nn=randn(1,N);%高斯白噪声fai=random('unif',0,1,1,2)*2*pi;%均匀随机变量xn=2*cos(2*pi*f1*t+fai(1))+5*cos(2*pi*f2*t+fai(2))+Nn;Xk=fft(xn);Sk=(abs(Xk).^2)/N;%随机信号分析书第二版上%第二问f=(0:N/2-1)*fs/N;w=10*log10(Sk(1:N/2));figure,plot(f,w);title('周期法估算随机信号功率谱');xlabel('f/Hz'),ylabel('Sx(f)(dB/Hz)');%第三问Sk=periodogram(xn);figure,plot(f,w);title('periodogram函数算随机信号功率谱');xlabel('f/Hz'),ylabel('Sx(f)(dB/Hz)');figure,subplot(2,1,1),plot(f,w);title('直接法估算随机信号功率谱');xlabel('f/Hz'),ylabel('Sx(f)(dB/Hz)');subplot(2,1,2),plot(f,w,'r');%延迟与原随机信号同一张图中比拟title('periodogram函数算随机信号功率谱')xlabel('f/Hz'),ylabel('Sx(f)(dB/Hz)');实验内容(2):间接法估计随机信号功率谱(1)计算以上的自相关函数。(2)通过计算自相关函数的Fourier变换,求的功率谱并绘图。(3)利用MATLAB函数psd、pwelch重新计算的功率谱,并与(2)做比拟。代码如下:clc;clear; %清屏及去除缓存N=1024;f1=30;f2=100;fs=1000; %序列长度和采样频率t=(0:N-1)/fs; %时间序列Nn=randn(1,N); %高斯白噪声fai=random('unif',0,1,1,2)*2*pi;%均匀随机变量xn=2*cos(2*pi*f1*t+fai(1))+5*cos(2*pi*f2*t+fai(2))+Nn;Rxx=xcorr(xn);Sk=abs(fft(Rxx));f=(0:N-1)*fs/N/2;w=10*log10(Sk(1:N));%第一问m=-1023:1023;figure,plot(m,Rxx); %画自相关函数title('自相关函数');xlabel('tao'),ylabel('R(tao)');%第二问figure,plot(f,w); %画功率谱title('自相关函数的功率谱');xlabel('f/Hz'),ylabel('Sx(f)(dB/Hz)');%第三问Nseg=256; %分隔字段为256window=hanning(Nseg); %汉宁窗noverlap=Nseg/2;f1=(0:Nseg/2)*fs/Nseg; %频率坐标轴Sx=psd(xn,Nseg,fs,window,noverlap,'none');figure,plot(f1,10*log10(Sx));gridon; %psd函数估计功率谱title('psd函数估计功率谱');xlabel('f/Hz'),ylabel('Sx(f)(dB/Hz)');Sx1=pwelch(xn,window,128,Nseg,fs,'onesided')*fs/2;%pwelch函数估计功率谱figure,plot(f1,10*log10(Sx1));gridon; title('pwelch函数估计功率谱');xlabel('f/Hz'),ylabel('Sx(f)(dB/Hz)');%三幅图的比照figure,subplot(3,1,1),plot(f,w/2);xlabel('f/Hz'),ylabel('Sx(f)(dB/Hz)');title('自相关函数的功率谱');subplot(3,1,2),plot(f1,10*log10(Sx));xlabel('f/Hz'),ylabel('Sx(f)(dB/Hz)');title('psd函数估计功率谱');subplot(3,1,3),plot(f1,10*log10(Sx1));xlabel('f/Hz'),ylabel('Sx(f)(dB/Hz)');title('pwelch函数估计功率谱');实验内容〔3〕:系统对随机信号响应的统计特性分析及仿真(1)生成含500点数据的高斯分布白噪声随机信号。(2)设计一个带通系统,其上、下截止频率分别为4KHz和3KHz。(3)计算通过以上带通滤波器的自相关函数和功率谱密度。代码如下:clearclc %清屏及去除缓存M=100;%样本数N=500;%样本长度N=500,对应时长25msxt=random('norm',0,1,M,N);%产生M×N个高斯随机数,相当于M个样本ht=fir1(101,[0.30.4]);%ht为带通滤波器的冲激响应HW=fft(ht,2*N);%滤波器频率响应Sxx=abs(fft(xt,2*N,2).^2)/(2*N);%周期图法估计M个白噪声样本的功率谱Sxxav=mean(Sxx);%M个样本的平均,对Sxx的行求平均HW2=abs(HW).^2;%系统的功率传输函数Syy=Sxxav.*HW2;%输出信号的功率谱Ryy=fftshift(ifft(Syy));%用IFFT求输出信号的自相关函数w=(1:N)/N;%功率谱密度横轴坐标t=(-N:N-1)/N*(N/20000);%自相关函数横轴坐标figure,plot(w,abs(Sxx(1:N)));title('输入信号的功率谱密度');xlabel('f/Hz'),ylabel('Sxx(f)');%输入信号的功率谱密度figure,plot(w,abs(HW2(1:N)));%系统的功率传输函数title('系统的功率传输函数');xlabel('f/Hz'),ylabel('H2(f)');figure,plot(w,abs(Syy(1:N)));%输出信号的功率谱密度title('输出信号的功率谱密度');xlabel('f/Hz'),ylabel('Sxy(f)');figure,plot(t,Ryy);%输出信号的自相关函数title('输出函数的自相关函数');xlabel('tao'),ylabel('Ryy〔tao)');实验数据
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 二零二五年度夫妻双方财产分配与子女抚养协议书
- 2025年度汽车修理厂与汽车维修培训学校合作协议
- 二零二五年度拆除承重墙工程风险评估与应急预案合同
- 二零二五年度参股协议书:新能源储能技术参股合作协议
- 二零二五年度影视基地租赁意向金合同模板
- 2025陕西阿康紫阳县农村供水管理有限公司招聘21人笔试参考题库附带答案详解
- 2025贵州榕晟体育文化产业有限责任公司招聘2人笔试参考题库附带答案详解
- 2025湖南高速工程咨询有限公司第一次招聘劳务派遣人员22人笔试参考题库附带答案详解
- 教师语言(吉林联盟)知到智慧树章节测试课后答案2024年秋吉林师范大学
- 2025年甘肃省民航机场集团劳务派遣工招聘58人笔试参考题库附带答案详解
- 《国别和区域研究专题》教学大纲
- 2025年日历表含农历(2025年12个月日历-每月一张A4可打印)
- 福建省部分地市2025届高中毕业班第一次质量检测 英语试卷(含答案)
- 《无创呼吸机护理》课件
- 2025年临床医师定期考核必考复习题库及答案(1060题)
- 腹膜透析管理标准操作规程
- 新人教版四年级数学下册全册教案 完整
- 2025年行政执法人员执法资格考试必考题库及答案(共232题)
- 浙教版小学三年级下册《劳动》全册教学课件
- 午休安全纪律教育
- 四年级全一册《劳动与技术》第一单元活动3《学习使用家用电器》课件
评论
0/150
提交评论