白化滤波器 matlab,白化滤波器-matlab-程序.doc
随机信号分析实验
PAGE
3 -
白化滤波器
原理
在统计信号处理中,往往会遇到等待处理的随机信号是非白色的,例如云雨、海浪、地物反射的杂乱回波等,它们的功率谱即使在信号通带内也非均匀分布。这样会给问题的解决带来困难。克服这一困难的措施之一是对色噪声进行白化处理。主要内容是设计一个稳定的线性滤波器,将输入的有色噪声变成输出的白噪声。
在这里,我们就对一般的具有功率谱的平稳随机过程X(t)白化处理问题进行讨论。为了具体的进行分析和计算,假设可以表达成有理数的形式,即
其中分子、分母为多项式。这个假设对于通常见到的功率谱是很近似的,而且有可行的方法用有理数去逼近任意的功率谱密度。
由于是功率谱,它的平稳随机过程相关函数的傅里叶变换具有非负的实函数和偶函数的性质。这些性质必然在其有理函数的表示式中体现出来,特别是,的零、极点的分布和数量会具有若干个特点。
由于是实函数,因此有:,是实数,的零、极点是共轭成对的。从而也可以把的表示式写成如下形式:
把开拓到复平面s中去,另。用s代替就可以把函数扩大到整个复平面。的零、极点必将对称于轴,如图13所示:
图13
由于是偶函数,因此不难判断,的零、极点是象限对称的,从而对于轴也是对称的。
由于,因此分子的虚根必然是偶数,否则会出现负值。这就是说轴上的零、极点必将成对的出现。
由于是可积的,因此分子的阶数不能大于分母的阶数,这就是说零点总数不会大于极点总数,而且分母不可能有虚根,这意味着轴上没有极点。
综合上述情况,在s平面的零、极点的可能位置如上图所示:
令:
则有
其中代表零、极点均在s左平面的部分,代表零、极点均在s右平面的部分。若在轴上有零点的话,必是成对的。则将一个放在内,将另一个放在内。实质上,对应的时域函数在负时间域为零,而对应的时域函数在正时间域为零。
根据上述的讨论,可以求得白化滤波器的解析式为:
由于
故得:
若运用傅里叶变换进行分析计算,以s代替,可得白化滤波器公式:
其中
我们知道,的傅里叶反变换是白化滤波器在时域的单位冲击响应,零、极点在s左半平面,因此的零、极点也是在s左半平面。故它对应的时域函数在负时域时为零,也就是说,上述白化滤波器是物理可实现的。
白化滤波器的设计方法是:
首先计算色噪声自相关函数,根据色噪声的自相关函数,计算出色噪声的功率谱(色噪声的自相关函数和功率谱构成一对傅里叶变换对),然后根据公式
程序
%*******色噪声的产生*******************%
Fs=44100;[x1,Fs]=wavread('E:\matlab\work\混合信号\色噪声.wav');
L1=length(x1)/10;
x=x1(1:44100);
l1=0:L1-1;
t=l1/Fs;
figure(1);
plot(t,x,'-r');
title('色噪声');
%******统计色噪声*******************
E=mean(x); %色噪声的均值E = -0.0054
S=var(x); %色噪声的方差S = 0.0324
%**************求色噪声概率密度函数*****************%
eachi=linspace(min(x),max(x),42);
yyi=hist(x,eachi); %计算各个区间的个数
yyi=yyi/length(x); %对各个区间的个数归一化处理
figure(8); %绘制色噪声的概率密度函数
plot(eachi,yyi,'-k')
title('色噪声的概率密度函数')
%***************色噪声自相关函数 *****%
Rx=xcorr(x,x); %色噪声的自相关函数Rx
tau=(-L1+1:L1-1)/Fs;
figure(2);
plot(tau,Rx,'-r')
title('色噪声的自相关函数'); %色噪声的自相关函数波形
xlabel('\tau'),ylabel('R_x(\tau)');
grid on;
hold on;
%*************** 色噪声功率谱密度 *****%
R=fft(Rx); %自相关函数的傅里叶变换即是功率谱密度
cm=abs(R);
fl=(0:length(R)-1)*44100/length(R);
figure(3)
plot(fl(1:length(fl)/2),cm(1:length(fl)/2),'-b')
title('色噪声的功率谱')
白化滤波器 matlab,白化滤波器-matlab-程序.doc相关推荐
- pid算法matlab仿真程序和c程序,pid算法matlab仿真程序和c程序.doc
pid算法matlab仿真程序和c程序.doc 1 增量式 PID 控制算 法 Matlab 仿真程序 设一被控对象 G(s) =50/(0.125s^2+7s),用增量式 PID 控制 算法编写仿真 ...
- dmc matlab程序,matlab编的DMC程序.doc
matlab编的DMC程序 clear all; % close all; %系统模型建立 num=[0.8]; den=[225 1]; [a,b,c,d]=tf2ss(num,den); % st ...
- 微机保护算法matlab编程,微机继电保护matlab算法仿真(有源程序).doc
微机继电保护matlab算法仿真(有源程序) 微 机 继 电 保 护 算 法 仿 真 电控学院 两点乘算法 matlab中编写的程序 N=12; t=(0:0.02/N:0.02)'; m=size( ...
- Matlab语音信号去噪程序,使用低通巴特沃斯滤波器
Matlab语音信号去噪程序,使用低通巴特沃斯滤波器. 1.读取一段歌曲的信号,绘制时域频域图,并播放. 2.添加正弦噪声: 3.设计巴特沃斯低通滤波器: 4.使用滤波器去除噪声,并画出时域频域图,播 ...
- 【 MATLAB 】使用 MATLAB 求由差分方程表示的滤波器的响应的两种方法
例题: 一个3阶低通滤波器由下面差分方程描述: y(n) = 0.0181 x(n) + 0.0543 x(n-1) + 0.0543 x(n-2) + 0.0181 x(n-3) + 1.76 y( ...
- fir滤波器应用matlab,fir滤波器应用.ppt
fir滤波器应用 电子与通信工程 许永全 FIR滤波器应用 内容 数字滤波器概述 FIR滤波器基本介绍 FIR滤波器在matlab及FPGA中的调用程序 FIR滤波器在无线信号处理的一个应用 数字滤波 ...
- 用matlab设计滤波器实验报告,数字信号出来实验报告--matlab滤波器设计
数字信号出来实验报告--matlab滤波器设计 广 西 工 学 院 实 验 报 告 用 纸 实验名称 IIR数字滤波器的设计 实验成绩 指导老师 陈艳 系(院) 计算机工程系 班级 学号 学生姓名 一 ...
- Matlab实现滤波器,进行ASK、FSK、多音信号的滤波
目录 绪论 一 :滤波器设计 二:ASK.FSK.多音信号产生 三:滤波 绪论 (1)以下是matlab破解版本的百度网盘资源: 链接:https://pan.baidu.com/s/1ORFFWE4 ...
- 数字信号处理matlab设计滤波器
目 录 1 设计目的.内容及要求 1 1.1 设计目的 1 1.2 设计内容 1 1.3 设计思考 2 1.4 设计要求 2 2 设计方案与设计原理 3 2.1 设计思路 3 2.2 采样定理 3 ...
- matlab无源滤波器,matlab在一阶无源带阻滤波器的应用
matlab在一阶无源带阻滤波器的应用 目录 1设计相关理论3 1.1串联谐振电路3 1.1.1电路模型分析3 1.1.2电路谐振条件4 1.1.3电路频率响应4 1.1.4电路品质因数5 1.2无源 ...
最新文章
- Datawhale组队学习:初级算法梳理课程任务
- HDU 5289 Assignment(单调队列)
- ADC/DAC设计常见40问
- 你大概走了假敏捷:认真说说敏捷的实现和问题(手绘版)
- Java 实现 HTTP 请求的三种方式
- DL之CNN:利用CNN算法实现对句子分类+进行情感分析(预测句子情感)
- 简单易用的.NET免费开源RabbitMQ操作组件EasyNetQ解析
- K/3 MRP运算数据不准的原因及解决方案
- spring 启动加载数据_12个很棒的Spring数据教程来启动您的数据项目
- 鸿蒙os开发小程序,9岁小学生展示鸿蒙OS开发:这操作太秀了
- 加密技术,给邮件安全加上一把锁
- 锁失效_关于bigtable中chubby锁失效时的一点思考
- html5标签 不可修改,SEO优化不可忽视的5个Html代码标签
- Go基础编程:运算符
- UAV021(六):系统架构优化、SBUS协议、遥控器控制电机转动
- java json转map_Java 把json对象转成map键值对的方法
- 维智WSD-A2系列伺服驱动器用户手册(MECHATROLINK-Ⅱ总线通信型)
- 压缩或解压文件出现循环冗余检查的解决办法
- 2022半入耳式蓝牙耳机怎么选?半入耳式耳机排行榜
- 计算机专业 英语复试专业问题准备(1)