xinhaoyuan.asv

给出了数字信号处理课程设计的基本程序

clear all
y1=zhengxian(1,2,32);subplot(4,2,1);stem(y1,'.');title('正弦信号');
             subplot(4,2,2);stem(abs(FFT(y1,32)),'.');%正弦信号
y2=maichong(16,8,2);subplot(4,2,3);stem(y2,'.');title('脉冲信号');
subplot(4,2,4);stem(abs(FFT(y2,32)),'.');%脉冲信号
y3=chongji(16,5,32);subplot(4,2,5);stem(y3,'.');title('冲激信号');
subplot(4,2,6);stem(abs(FFT(y3,32)),'.');%冲激信号
y4=gaoshi(32);subplot(4,2,7);stem(y4,'.');title('高斯噪声');
subplot(4,2,8);stem(abs(FFT(y4,32)),'.');  %高斯噪声
figure();
y=zhengxian(1,2,32);subplot(3,2,1);stem(abs(FFT(y,8)),'.');title('F= 2, N= 32,2M=8'); %正弦信号
 y=zhengxian(1,2,32);subplot(3,2,2);stem(abs(FFT(y,16)),'.');title('F= 2, N= 32,2M=16'); %正弦信号
 y=zhengxian(1,2,32);subplot(3,2,3);stem(abs(FFT(y,32)),'.');title('F= 2, N= 32,2M=32'); %正弦信号
 y=zhengxian(1,2,32);subplot(3,2,4);stem(abs(FFT(y,64)),'.');title('F= 2, N= 32,2M=64'); %正弦信号
 y=zhengxian(1,2,32);subplot(3,2,5);stem(abs(FFT(y,128)),'.');title('F= 2, N= 32,2M=128'); %正弦信号
 y=zhengxian(1,2,32);subplot(3,2,6);stem(abs(FFT(y,256)),'.');title('F= 2, N= 32,2M=256'); %正弦信号
 figure();%整周期抽样和非整周期抽样的频谱
  y=zhengxian(1,1,32);subplot(3,2,1);stem(abs(FFT(y,32)),'.');title('F= 1, N= 32,2M=32'); %正弦信号
  y=zhengxian(1,1.3,32);subplot(3,2,2);stem(abs(FFT(y,32)),'.');title('F= 1.3, N= 32,2M=32'); %正弦信号
  y=zhengxian(1,2,32);subplot(3,2,3);stem(abs(FFT(y,32)),'.');title('F= 2, N= 32,2M=32'); %正弦信号
  y=zhengxian(1,2.4,32);subplot(3,2,4);stem(abs(FFT(y,32)),'.');title('F= 2.4, N= 32,2M=32'); %正弦信号
  y=zhengxian(1,4,64);subplot(3,2,5);stem(abs(FFT(y,64)),'.');title('F= 4, N= 64,2M=64'); %正弦信号
  y=zhengxian(1,3.5,64);subplot(3,2,6);stem(abs(FFT(y,64)),'.');title('F= 3.5, N= 64,2M=64'); %正弦信号
  figure();%高斯噪声功率谱分析
  y=gaoshi(32);subplot(3,1,1);stem(y,'.');title('高斯噪声');%高斯噪声
  u=R(y);