softautumn
晨714
采纳率75%
2021-04-05 11:42

FIR滤波器设计的不调用fir函数的matlab,理想滤波的h序列引入了时延和窗函数时波形发生畸变

10

        fs=8000;%采样频率
Fp=1000;
Fs=1500;
wp=2*pi*Fp/fs;
ws=2*pi*Fs/fs;
wc=1/2*(wp+ws);
          N2=ceil(6.2*pi/(ws-wp));%h序列的长度
          tao=ceil((N2-1)/2);%w序列的时延

          n2=1:N2;%h序列的横坐标点阵
          
W=0.5*(1-cos(2*pi*n2/(N2-1)));%line12,窗函数
h=sin(wc*(n2-tao))./(pi*(n2-tao)).*W;%line13,若引入此行的窗函数,或者去掉窗函数但引入时延tao,波形都会畸变。
%h=sin(wc*(n2))./(pi*(n2));  %line14,若注释掉上两行,只采用理想滤波的序列,即此行,效果还行。
len_h=length(h);


num_max=100;%设置输入信号的点数,即x的横坐标点阵
max=1/fs*num_max;

t=1/fs:1/fs:max;
x=sin(2*pi*600*t);%据此设置频率,观察滤波程度
len_x=length(x);

y=zeros(1,len_x);
%卷积x和h,经测试,下面部分是没有问题的。
for num=1:len_x  
    temp=0;
    for i=1:len_h
            if(num+1-i>0)
                 temp=h(i).*x(num+1-i)+temp;
            end
            if(num+1-i<=0)
                  temp=temp;
            end
    end
    y(num)=temp;
end
plot(t,y)
%如果只用line14,注释掉line12和13,可以达到滤波效果。只是幅值衰减没那么明显而已。
%但是如果引入了窗函数和时延,各种频率就会出现波形畸变,而且幅值也不会衰减,求大佬指教,修改下代码为正确可行
%要求:不调用matlab本身自带的fir函数(因为我必须用c语言设计)

  • 点赞
  • 写回答
  • 关注问题
  • 收藏
  • 复制链接分享
  • 邀请回答

2条回答

  • qq_44791407 fightandstrive 1月前

    你引入时延的方式不对,引入时延h序列第1到tao个原始应该置零。在13行后加入h(1:tao)=0即可

    点赞 1 评论 复制链接分享
  • qq_44791407 fightandstrive 1月前

     

    点赞 评论 复制链接分享

相关推荐