晨714 2021-04-05 11:42 采纳率: 100%
浏览 73
已采纳

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

        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条回答 默认 最新

    报告相同问题?

    相关推荐 更多相似问题

    悬赏问题

    • ¥20 有没有认识硅云的研发 求介绍
    • ¥15 unittest __import__加载名为‘1’的model无法找到
    • ¥30 魔术师读心术猜数字实验
    • ¥15 关于#Python#的问题,如何解决?
    • ¥15 CMD 生成TXT内容输出成一行
    • ¥20 关于#pandas#的问题,如何解决?(语言-python)
    • ¥30 c语言设计井字棋小游戏
    • ¥15 关于属性关联的双极容度多属性决策vikor方法的公式
    • ¥15 10吨190马力爆炸品运输车,以80km/h的速度行驶1个小时输出的动能或者做的功是多少?
    • ¥15 《phthon编程从入门到实践》第十八章,从django入手,learning_log,ModuleNotFoundError: No module named 'learning_logs