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语言设计)