ww.m

算法计算n=1到N的Fourier系数偶次谐波抽样点上的直流分量

%program4_2
N=8;
%计算n=-N到－1的Fourier系数
n1=-N:-1;
c1=-4*j*sin(n1*pi/2)/pi^2./n1.^2;
%计算n=0时的Fourier系数
c0=0;
%计算n=1到N的Fourier系数
n2=1:N;
c2=-4*j*sin(n2*pi/2)/pi^2./n2.^2;
cn=[c1 c0 c2];
n=-N:N;
subplot(2,2,1);
stem(n,abs (cn));
ylabel('Cn的幅度');
xlabel('\omega/\omega0');

subplot(2,2,2);
stem(n,angle (cn));
ylabel('Cn的相位');
xlabel('\omega/\omega0');
%program4_3
t=-2:0.001:2;%信号的抽样点
N=7;
c0=0.5;
fN=c0*ones(1,length(t));%计算抽样点上的直流分量
for n=1:2:N %偶次谐波为零
fN=fN+cos(pi*n*t)*sinc(n/2);
end
subplot(2,2,3)
plot(t,fN);
title (['N=',num2str(N)])
axis([-2 2 -0.2 1.2]);

title('集成05-3 孙振宇 12');