信号与线性系统实验报告(6)

Ak=int(y*exp(-2*i*pi*k*t/TT),t,-TT/2,TT/2)/TT;

% 利用symmul命令得到Ak*exp(2*j*k*pi*t/T) 的符号表达式 %fk=symmul(Ak,sym(exp(2*i*k*pi*t/TT))); %fk=symmul(Ak,exp(2*i*k*pi*t/TT)); fk=sym(Ak)*sym(exp(2*i*k*pi*t/TT));

% 计算 k 在区间 [-N,N] 内的 Fourier 系数，存入向量 a（长度为 2*N+1） for m=-N:-1

% a(m+N+1)=numeric(subs(Ak,k,m)); a(m+N+1)=double(subs(Ak,k,m)); end

%a(N+1)=numeric(A0); a(N+1)=double(sym(A0)); for m=1:N

a(m+N+1)=double(subs(Ak,k,m)); end % for

% 利用symmul命令进行符号求和，得到综合信号 f f=symsum(fk,k,-N,-1)+A0+symsum(fk,k,1,N); % 图形显示结果? figure n=-N:N;

as=abs(a)*2; % 由 Fourier 系数得到频谱幅度 subplot(3,1,1) ezplot(y,[-TT,TT]) ylabel('原函数'); subplot(3,1,2) ezplot(f,[-TT,TT]) ylabel('合成函数'); subplot(3,1,3) stem(n,as)

ylabel('幅度频谱图');

仿真结果为：