声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1037|回复: 2

[滤波] 求助一个滤波器生成函数

[复制链接]
发表于 2008-5-16 20:24 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
本帖最后由 wdhd 于 2016-9-10 14:09 编辑

  在下面附件1中有待求的公式。
  附件2是我写的程序,但是在有些数据时会出错
  比如:
  K=16;
  M=190;
  B=10;
  t=0:1/(M*K):(K-1); %序列横座标
  N=length(t); %序列总长
  fs=1/(M*K);
  n1=randn(1,N);
  t1=t;
  b=[0.2,0.4,0.6,0.8,1.0,1.0,1.0,0.7,0.5,0.3,0.1];
  [n2,tn2]=xgzs(n1,t1,b,fs);
  n2=B*n2;

待求的式子

待求的式子

xgzs.m

327 Bytes, 下载次数: 0

我的程序

回复
分享到:

使用道具 举报

 楼主| 发表于 2008-5-16 20:35 | 显示全部楼层

回复 楼主 的帖子

因为可能存在看不到附件的情况,我把程序在下面写出来,请大家帮忙看下!

function [n2,tn2]=xgzs(n1,tn1,b,fs)
% 知道噪声n1(t)的情况下求其相关噪声n2(t)
zc=length(b);                           %序列b(n)的长度
n2=0;
t2=0;
for i=0:1:(zc-1)
        n0=n2;
    t0=t2;
    [n3,t3]=Delay_sig0(n1,tn1,i,fs);
        n2=b(i+1)*n3;
    t2=t3;   
        [n2,t2]=Sig_Add(n0,t0,n2,t2,1,1,fs);
    i
end
n2=n2;
tn2=t2;
 楼主| 发表于 2008-5-16 20:36 | 显示全部楼层

回复 楼主 的帖子

function [y,n]=Delay_sig(x,m,n0,fs)                %  n0 为延迟的样点数
n=m-n0*fs;
y=x;

function [y,n]=Sig_Add(x1,n1,x2,n2,a,b,fs)
n=min(min(n1),min(n2)):fs:max(max(n1),max(n2));
y1=zeros(1,length(n));
y2=y1;
y1(find((n>=min(n1))&(n<=max(n1))==1))=x1;
y2(find((n>=min(n2))&(n<=max(n2))==1))=x2;
y=a*y1+b*y2;
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-19 04:50 , Processed in 0.058390 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表