周期信号的傅里叶级数分析
例一:
代码如下:
t=-10:0.01:10;
tao=4;T=8;//(1)情形一
w0=2*pi/8;
k=(-10:0.25:10)*(w0/pi);
f=square((pi/4)*(t+2),50)/2+0.5;
subplot(2,2,1),plot(t,f),grid on
Ak=tao/T*sinc(k*tao/2);//周期信号矩形脉冲信号的傅里叶级数系数表达式
subplot(2,2,2),stem(k,Ak),axis([-5 5 -0.5 0.5]),grid on
tao=2;T=8;//(2)情形二
Ak=tao/T*sinc(k*tao/2);
subplot(2,2,3),stem(k,Ak),axis([-5 5 -0.25 0.2]),grid on
tao=1;T=8;//(3)情形三
Ak=tao/T*sinc(k*tao/2);
subplot(2,2,4),stem(k,Ak),axis([-5 5 -0.125 0.1]),grid on
例二:
代码如下:
t=-5:0.001:5;//时间变量
xt=5*cos(t)+2*cos(10*t);
subplot(3,1,1),plot(t,xt),ylabel('x(t)'),xlabel('Time(s)'),grid on
w1=1;w2=10;
H1=1/((i*w1)^2+3*i*w1+2);
H2=1/((i*w2)^2+3*i*w2+2);
y=5*abs(H1)*cos(w1*t+angle(H1))+ 2*abs(H2)*cos(w2*t+angle(H2));
subplot(3,1,2),plot(t,y),ylabel('y(t)'), xlabel('Time(s)'),grid on