randb.m

十几个数字信号处理实验课的程序

n=0:63;
xb=exp(-0.1*n).*sin(2*pi*0.0625*n);
w=randn(1,64);
y=xb+w;
figure
subplot(2,1,1);stem(xb);title('xb(n)');
subplot(2,1,2);stem(y);title('y(n)');
figure
subplot(2,1,1);stem(abs(fft(xb,64)));title('xb的幅度谱');
subplot(2,1,2);stem(abs(fft(y,64)));title('y的幅度谱');

figure
subplot(2,1,1);stem(angle(fft(xb,64)));title('xb的相位谱');
subplot(2,1,2);stem(angle(fft(y,64)));title('y的相位谱');