samp3_1_2.m

来自「数字信号处理的Matlab实现的随书光盘」· M 代码 · 共 16 行

%Samp3_1_2
if(mod(N,2)~=1)a(m)=a(m)/2; end   %此时b(m)为零，a(m)减半
for ii=0:N-1
    xx(ii+1)=a(1)/2;
    for k=1:m-1
        xx(ii+1)=xx(ii+1)+a(k+1)*cos(2*pi*k*ii/N)+b(k+1)*sin(2*pi*k*ii/N);
% 
    end
 end
subplot(2,1,1)
plot((0:N-1)*dt,x)  %绘制原信号
title('原始信号')
subplot(2,1,2)
plot((0:N-1)*dt,xx)   %绘制合成信号
title('合成信号'),xlabel('时间/s')