首页 范文大全 古典文学 职场知识 中国文学 公文书信 外国名著 寓言童话 百家讲坛 散文/诗歌 美文欣赏 礼仪知识 民俗风情
  • 工作总结
  • 工作计划
  • 心得体会
  • 竞聘演讲
  • 会议发言
  • 爱国演讲
  • 就职演说
  • 开业开幕
  • 思想学习
  • 征文演讲
  • 经验材料
  • 述职报告
  • 调研报告
  • 工作汇报
  • 年终总结
  • 申报材料
  • 学习体会
  • 企划方案
  • 活动方案
  • 技巧经验
  • 模板范例
  • 思想宣传
  • 经济工作
  • 工作报告
  • 组织人事
  • 反腐倡廉
  • 慰问贺电
  • 先进事迹
  • 思想汇报
  • 入党申请书
  • 党会发言
  • 先进性教育
  • 入团申请书
  • 个人简历
  • 演讲稿
  • 调查报告
  • 实习报告
  • 和谐社会
  • 观后感
  • 读后感
  • 作文范文
  • 自我鉴定
  • 讲话稿
  • 自查报告
  • 维纳滤波器设计随机信号处理

    时间:2021-03-17 18:05:57 来源:蒲公英阅读网 本文已影响 蒲公英阅读网手机站

    相关热词搜索:滤波器 信号处理 随机

    维纳滤波器设计随机信号处理 本文关键词:滤波器,信号处理,随机,维纳,设计

    维纳滤波器设计随机信号处理 本文简介:随机信号处理实验报告维纳滤波器设计随机信号处理学号:姓名:实验一维纳滤波器设计1实验内容设计一个维纳滤波器:(1)产生三组观测数据,首先根据产生信号,将其加噪,(信噪比分别为),得到观测数据。(2)估计,的AR模型参数。假设信号长度为L,AR模型阶数为N,分析实验结果,并讨论改变L,N对实验结果的影

    维纳滤波器设计随机信号处理 本文内容:

    随机信号处理实验报告

    维纳滤波器设计

    随机信号处理

    学号:

    姓名:

    实验一

    维纳滤波器设计

    1实验内容

    设计一个维纳滤波器:

    (1)

    产生三组观测数据,首先根据产生信号,将其加噪,(信噪比分别为),得到观测数据。

    (2)

    估计,的AR模型参数。假设信号长度为L,AR模型阶数为N,分析实验结果,并讨论改变L,N对实验结果的影响。

    2实验原理

    滤波目的是从被噪声污染的信号中分离出有用的信号来,最大限度地抑制噪声。对信号进行滤波的实质就是对信号进行估计。滤波问题就是设计一个线性滤波器,使得滤波器的输出信号是期望响应的一个估计值。下图就是观测信号的组成和信号滤波的一般模型。

    观测信号

    信号滤波的一般模型

    维纳滤波解决从噪声中提取信号的滤波问题,并以估计的结果与真值之间的误差均方值最小作为最佳准则。它根据估计信号的当前值,它的解以系统的系统函数或单位脉冲形式给出,这种系统常称为最佳线性滤波器。

    维纳滤波器设计的任务就是选择,使其输出信号与期望信号误差的均方值最小。

    假设滤波系统是一个线性时不变系统,它的和输入信号都是复函数,设

    考虑系统的因果性,可得到滤波器的输出

    设期望信号,误差信号及其均方误差分别为

    要使均方误差为最小,需满足:

    整理得,等价于

    上式说明,均方误差达到最小值的充要条件使误差信号与任一进入估计的输入信号正交,这就是正交性原理。

    将展开,得

    整理得

    等价于

    此式称为维纳-霍夫(Wiener-Holf)方程。解此方程可得到最优权系数,此式是Wiener滤波器的一般方程。

    定义

    则维纳-霍夫方程可写成矩阵形式

    求逆,得

    此式表明,已知期望信号与观测数据的互相关函数及观测数据的自相关函数时,可以通过矩阵求逆运算,得到维纳滤波器的最佳解。

    3实验结果及分析

    (1)当L=200,N=6

    信噪比为20dB的滤波效果

    信噪比为10dB的滤波效果

    信噪比为6dB的滤波效果

    (2)当L=200,N=60

    信噪比为20dB的滤波效果

    信噪比为10dB的滤波效果

    信噪比为6dB的滤波效果

    (3)当L=600,N=6

    信噪比为20dB的滤波效果

    信噪比为10dB的滤波效果

    信噪比为6dB的滤波效果

    实验分析:别取信号长度为200、600,滤波器长度为6、60,加噪信噪比为20dB、10dB、6dB,组合进行实验。每组实验得到的最小均方误差统计如下表。由此表可以看出,信号长度越长,最小均方误差(绝对值)越大,精度越差;在信噪比较大(误差影响较小)的滤波过程中,滤波器长度约长,最小均方误差(绝对值)越小,精度越好。

    表1

    最小均方误差统计表

    加噪

    信号

    信噪比

    L=200

    N=6

    L=200

    N=60

    L=600

    N=6

    20dB

    -1785.3971

    -1574.1107

    -6888.4103

    10dB

    -845.5714

    -1727.5260

    -10219.3397

    6dB

    -1213.6367

    -2135.4596

    -7743.9358

    对于相同信号和滤波器(这里取L=200,N=6),信噪比越大,最小距离误差约小;而当信噪比较小时,信号与噪声值接近,导致滤波效果受到影响,最小距离误差变大。

    4源代码

    clear;

    clc;

    %初始化变量

    L=200;

    %信号长度

    N=6;

    %滤波器的阶次

    a=0.96;

    wn=randn(L,1);

    %wn为用于生成信号的噪声信号,随机生成一个L*1矩阵,生成矩阵的元素值在%区间(0.0,1.0)之间

    sn=zeros(L,1);

    %sn为信号,生成一个L*1的零矩阵

    hn=zeros(N,1);

    %hn为系统单位脉冲响应

    生成一个N*1的零矩阵

    rxx=zeros(N,1);

    %rxx为自相关函数,生成一个N*1的零矩阵

    rxd=zeros(N,1);

    %rxd为互相关函数,生成一个N*1的零矩阵

    yn=zeros(L,1);

    %yn为输出信号,生成一个L*1的零矩阵

    xt=zeros(L+N,1);

    %生成一个(L+N)*1的零矩阵

    gn=zeros(L,1);

    %gn为yn与sn最小距离误差信号,生成一个L*1的零矩阵

    %根据给定公式s(n)=as(n-1)+w(n),生成理想信号

    for

    i=2:L

    sn(i,1)=a*sn(i-1,1)+wn(i,1);

    end

    sn(1,1)=wn(1,1);

    subplot(2,2,1);

    plot(sn,r

    ),axis([0

    200

    -10

    10]),xlabel(

    时间

    ),ylabel(

    幅度

    ),title(

    sn

    )

    %生成期望信号方差cd

    cd=(var(wn))/(1-a^2);

    %对信号加噪

    x1=awgn(sn,20);

    x2=awgn(sn,10);

    x3=awgn(sn,6);

    subplot(2,2,2)

    plot(x3,g

    ),axis([0

    200

    -10

    10]),xlabel(

    时间

    ),ylabel(

    幅度

    ),title(

    x3

    );

    %生成输入信号与理想信号的互相关函数,此处x1为输入信号,sn为期望信号

    for

    i=1:N,for

    m=i:1:L,rxd(i,1)=rxd(i,1)+x3(m,1)*sn(m-i+1,1);

    end

    end

    %生成输入信号的自相关函数

    for

    i=1:N,for

    m=i:1:L,rxx(i,1)=rxx(i,1)+x3(m,1)*x3(m-i+1,1);

    end

    end

    %将自相关函数生成托普勒斯矩阵

    rxx1=toeplitz(rxx);

    %生成逆矩阵

    irxx=inv(rxx1);

    %生成滤波器系数h(n)

    hn=irxx*rxd;

    for

    i=1:L

    xt(i+N,1)=x3(i,1);

    end

    %实现滤波

    for

    i=1:L,for

    m=1:N,yn(i,1)=yn(i,1)+xt(i+N+1-m,1)*hn(m,1);

    end

    end

    %计算最小均方误差信号en

    en=0;

    en=cd-(rxd

    )*hn;

    %生成最小距离误差信号gn

    gn=yn-sn;

    %画出滤波后的信号时域图

    subplot(2,2,3);

    plot(yn),axis([0

    200

    -10

    10]),xlabel(

    时间

    ),ylabel(

    幅度

    ),title(

    yn

    );

    %画出理想信号与输出信号对比图

    subplot(2,2,4);

    plot(sn,r

    ),axis([0

    200

    -10

    10]),xlabel(

    时间

    ),ylabel(

    幅度

    ),title(

    sn与yn对比

    );

    hold

    on;

    plot(yn,b

    ),axis([0

    200

    -10

    10]);

    hold

    off;

    %画出最小距离误差信号图

    figure;

    plot(gn),axis([0

    200

    -2

    2]),xlabel(

    时间

    ),ylabel(

    幅度

    ),title(

    gn

    );

    10

    • 范文大全
    • 职场知识
    • 精美散文
    • 名著
    • 讲坛
    • 诗歌
    • 礼仪知识