数字信号处理实验报告

');gridtitle('Phase(36HzSamplingrequency)');xlabel('f');ylabel('???¡§Xk??');fs=100;dt=1/fs;xn=sin(2*pi*f1*n*dt)+2*sin(2*pi*f2*n*dt);Xk=fft(xn,N);mag=2*abs(Xk)/N;pha=angle(Xk);f=fs/N*n;subplot(223)plot(f,mag,'r');gridtitle('Magnitude(100HzSamplingFrequency)');xlabel('f');ylabel('|Xk|');subplot(224)plot(f,pha,'r');gridtitle('Phase(100HzSamplingFrequency)');xlabel('f');ylabel('???¡§Xk??');N=64;n=0:N-1;f1=15;f2=18;fs=50;dt=1/fs;xn=sin(2*pi*f1*n*dt)+2*sin(2*pi*f2*n*dt);Xk=fft(xn,N);mag=2*abs(Xk)/N;pha=angle(Xk);f=fs/N*n;figure(3);subplot(221)plot(f,mag,'m');gridtitle('Magnitude(64SamplingPoints)');xlabel('f');ylabel('|Xk|');subplot(222)plot(f,pha,'m');gridtitle('Phase(64SamplingPoints)');xlabel('f');ylabel('???¡§Xk??');N=512;n=0:N-1;dt=1/fs;xn=sin(2*pi*f1*n*dt)+2*sin(2*pi*f2*n*dt);Xk=fft(xn,N);