fir1.m

呵呵,第一次来很多东西不懂 还请多多指教啊! 几个编的简单的FFT和滤波器设计 希望能有点用

clear all
[x,y] = solve('abs(2*x*cos(0.1)+y)','abs(2*x*cos(0.4)+y-1)') 
subplot(311);
n=-10:10;
syms w;
a=sum((cos(0.1*n)+cos(0.4*n)).*exp(-j*w*n))
w=0:0.1:3;
a=-1020887807958581/9007199254740992*exp(10*i*w)-2478316896340541/9007199254740992*exp(9*i*w)-2716463807345489/9007199254740992*exp(8*i*w)-798849192683869/4503599627370496*exp(7*i*w)[phone]/2251799813685248*exp(6*i*w)+8312487042844827/18014398509481984*exp(5*i*w)+4016586991621411/4503599627370496*exp(4*i*w)+5934367304475473/4503599627370496*exp(3*i*w)[phone]/281474976710656*exp(2*i*w)+2157297584338775/1125899906842624*exp(i*w)+2+2157297584338775/1125899906842624*exp(-i*w)[phone]/281474976710656*exp(-2*i*w)+5934367304475473/4503599627370496*exp(-3*i*w)+4016586991621411/4503599627370496*exp(-4*i*w)+8312487042844827/18014398509481984*exp(-5*i*w)[phone]/2251799813685248*exp(-6*i*w)-798849192683869/4503599627370496*exp(-7*i*w)-2716463807345489/9007199254740992*exp(-8*i*w)-2478316896340541/9007199254740992*exp(-9*i*w)-1020887807958581/9007199254740992*exp(-10*i*w);
plot(w,a);
w=0:0.1:1;
b=(2.*-6.761949634801.*cos(w)+13.4563361040).*exp(-j*w);
subplot(312)
plot(w,b)
w=0:0.1:2;
c=-1020887807958581/9007199254740992*exp(10*i*w)-2478316896340541/9007199254740992*exp(9*i*w)-2716463807345489/9007199254740992*exp(8*i*w)-798849192683869/4503599627370496*exp(7*i*w)[phone]/2251799813685248*exp(6*i*w)+8312487042844827/18014398509481984*exp(5*i*w)+4016586991621411/4503599627370496*exp(4*i*w)+5934367304475473/4503599627370496*exp(3*i*w)[phone]/281474976710656*exp(2*i*w)+2157297584338775/1125899906842624*exp(i*w)+2+2157297584338775/1125899906842624*exp(-i*w)[phone]/281474976710656*exp(-2*i*w)+5934367304475473/4503599627370496*exp(-3*i*w)+4016586991621411/4503599627370496*exp(-4*i*w)+8312487042844827/18014398509481984*exp(-5*i*w)[phone]/2251799813685248*exp(-6*i*w)-798849192683869/4503599627370496*exp(-7*i*w)-2716463807345489/9007199254740992*exp(-8*i*w)-2478316896340541/9007199254740992*exp(-9*i*w)-1020887807958581/9007199254740992*exp(-10*i*w).*(2.*-6.761949634801.*cos(w)+13.4563361040).*exp(-j*w);
subplot(313)
plot(w,c)