随机信号通过带通滤波器

题目描述:
信号为两个正弦信号之和,一个振幅为1,频率为60Hz,另一个振幅为0.5,频率为40Hz;加性噪声是均值为0,方差为4的高斯白噪声,一起调制到9.5kHz的载波上,经过带宽为200Hz、中心频率为9.5kHz的带通系统,求输出信号的功率谱估计。采样频率取20kHz,滤波器阶数选200。
题目分析:仿真思想非常容易理解,废话不多说,直接上代码,不懂的部分请看代码注释:

fs=20000;               %采样频率
f0=9500;                %中心频率
N=50000;                %仿真点数
deltf=200;              %带宽
M=200;                  %滤波器阶数
f1=60;                  %双音信号的频率
f2=40;
t=0:1/fs:(N-1)/fs;      %时间序列
x0=sin(2*pi*f1*t)+0.5*sin(2*pi*f2*t);%原始信号
x=x0+random('norm',0,2,1,N);%原始信号叠加高斯噪声
xt=x.*cos(2*pi*f0*t);   %调制
f1=f0*2/fs;
df1=deltf/fs;
ht=fir1(M,[f1-df1,f1+df1]);%加带通滤波器
X=conv(xt,ht);          %经过滤波器之后的输出
t=0:1/fs:(N+M-1)/fs;
Nt=X.*cos(2*pi*f0*t);   %解调
h1=fir1(M,deltf/fs);    %加低通滤波器
N1=N+M;
Hw=fft(h1,N1);
Ntw=fft(Nt,N1);
RXw=(abs(fft(X,N1))).^2/N1;  %解调前功率谱密度
RYw=(abs(Ntw.*Hw)).^2/N1;    %解调后功率谱密度
maxw=max(max(RXw),max(RYw));
RXw=10*log10(RXw/maxw+0.001);%归一化解调前功率谱密度
RYw=10*log10(RYw/maxw+0.001);%归一化解调后功率谱密度
f=fs/N1:fs/N1:fs/2;
subplot(121);                %对比功率谱密度变化
plot(f,RXw(1:N1/2));
axis([9100  9900  -30  0]);
title('spectrum density');
xlabel('frequency/Hz');ylabel('Amplitude/dB');
subplot(122);
plot(f,RYw(1:N1/2));
axis([0 200  -30  0]);
title('spectrum density after demodulation')
xlabel('frequency/Hz');ylabel('Amplitude/dB');

仿真结果如下:

随机信号通过带通滤波器相关推荐

  1. 【数字信号处理】相关函数应用 ( 相关函数应用场景 | 噪声中信号检测 | 隐含周期性检测 | 时差估计 | 描述随机信号 )

    文章目录 一.相关函数应用场景 一.相关函数应用场景 相关函数应用场景 : " 噪声 " 中 信号检测 : 发射信号 通过 信道 后 , 假如该信道 是 无线信道 , 如手机的 4 ...

  2. matlab生成常用信号(方波、三角波、随机信号、单位冲激)

    使用Octave在线: 1 方波 y=square(t,DUTY) 产生一个时长为t.幅值为±1的周期性方波信号,其中的DUTY表示占空比,即在信号的一个周期中正值所占的百分比: t=-2*pi/10 ...

  3. 通信原理2习题课汇总(随机信号、信道、最佳接收)

    写在前面:本文仅用于个人学习,授课教师为郭心悦,下面的ppt为郭老师上课课件. 文章目录 第二章:随机信号习题 第三章:信道习题课 第五章:数字信号最佳接收习题课 第二章:随机信号习题 自相关函数:离 ...

  4. 通信原理随机信号分析

    通信原理第二章 随机信号分析 一 随机过程 定义 测试n台性能相同的接收机,在同样条件下,不加信号测试其输出噪声,波形如图 (1)每一条曲线 ξi(t)\xi_i(t)ξi​(t) 都是一个随机起伏的 ...

  5. Matlab生成各种伪随机信号_idinput指令

    1.为何不是冲激信号? 只要是在信号系统这条道上混过的童鞋都知道,冲激信号是测试系统的最理想信号.因为时域的卷积相当于频域的乘积,而冲击信号的频域拥有最广泛的频谱,可以"暴露"系统 ...

  6. 随机信号的傅里叶分析

    随机过程可以作为离散时间信号的模型. 通常,一个随机过程是一族带有序号的随机变量: ..., x[-2], x[-1], x[0], x[1], x[2], ... 上面的每一个x[i]都是一个随机变 ...

  7. n阶自相关matlab代码,随机信号及其自相关函数和功率谱密度的MATLAB实现.doc

    随机信号及其自相关函数和功率谱密度的MATLAB实现 引言: 现代信号分析中,对于常见的具有各态历经的平稳随机信号,不可能用清楚的数学关系式来描述,但可以利用给定的N个样本数据估计一个平稳随机信号的功 ...

  8. 建模股票价格数据并进行预测(统计信号模型):随机信号AR模型+Yule-Walker方程_Python...

    1.背景: 针对股票市场中AR 模型的识别.建立和估计问题,利用AR 模型算法对股票价格进行预测. 2.模型选取: 股票的价格可视为随机信号,将此随机信号建模为:一个白噪声通过LTI系统的输出,通过原 ...

  9. 20210701:随机信号的功率谱估计相关算法的实现

    随机信号的功率谱估计相关算法的实现 今天最重要的事!!! 问题描述 前置数学原理 BT法与周期图法估算功率谱 估算功率谱的示例 写在最后 今天最重要的事!!! 1. 我党100周年生日,中国共产党nb ...

最新文章

  1. 2003網域升級到2008網域以及遷移DNS
  2. .NET Winform也能画出类似QQ、飞信这样的窗口风格和控件效果
  3. 用Docker创建Nexus
  4. SDOI2020游记
  5. Orchard商城模块(Commerce)设计与后台部分
  6. 計算機二級-java-03
  7. 东南大学成贤c语言试卷,成贤高数(下)期中试卷参考答案
  8. 深入理解JavaScript中的属性和特性
  9. sqlprompt安装注册
  10. 如何下载斗鱼回放视频(纯手工)
  11. 医学图像分类 神经网络,神经网络图像识别技术
  12. 【软件工程】数据库设计说明书
  13. chrome谷歌浏览器script标签引入CLODOP CLodopfuncs.js 跨域报错
  14. 国家对五险一金的交纳说明(已更新)
  15. Android切词工具——BreakIterator(1)
  16. Typo: In word xxxx
  17. 17万用来助学,能做什么?
  18. 【测试开发】一个5年测试开发的成长经验,大学毕业就开启他的职业生涯......
  19. python win32gui+pynput 自动回复qq消息
  20. iOS开发笔记之二十四——Xcode下类中供外部调用的方法添加注释说明技巧

热门文章

  1. 哥伦比亚大学计算机工程面试题
  2. rhel5.5_Apache配置openssl支持https服务
  3. ASP.NET MVC导出excel(数据量大,非常耗时的,异步导出)
  4. Win10光驱不见了
  5. JS获取浏览器高度 并赋值给类
  6. [读书笔记]C#学习笔记二: 委托和事件的用法及不同.
  7. Kotlin教程学习-数据类型
  8. JVM参数详解和调优
  9. Centos6.5_x86_64 编译安装Subversion-1.8.5
  10. NetBeans 时事通讯(刊号 # 125 - Nov 17, 2010)