数字信号处理fft快速傅里叶变换 程序仿真.docVIP

北京信息科技大学硕士研究生 数字信号处理 课程仿真实验报告 内容 仿真结果 报告文字和分析 格式 总分 分数 学科 控制工程 姓名 崔 领 袖 班级 研 1507 学号 2015020207 2015 年 月 日 实验目的 1. 通过本次实验熟悉用Matlab产生方波以及白噪声。 2. 学会运用数字信号处理的知识进行采样，并从采的的点数不同来更加深刻的理解信号采样定理。 3. 学会运用FFT处理分析信号问题，并充分学习理解频谱在信号中的应用。 4. 增强独立完成实验的能力。 实验要求 首先生成幅值为1，-1，各持续半个周期，周期T=16ms的方波信号；然后生成均值为0，方差为0.1，强度为0.1的白噪声信号。将白噪声信号加在方波信号进而得到复合信号。 （1）对进行采样，采样周期为1ms，采用FFT算法分别计算128点、256点、1024点采样信号的频谱。 （2）从上述结果中如何估算信号的各次谐波分量，给出分析结果，并与没有噪声情况下的理论分析结果对照。 2、首先生成幅值为1，周期为16ms的方波信号，然后生成均值为0，方差分别为0,0.05，0.1，0.2的白噪声信号，将得到方波信号进行计算并与白噪声相互叠加得到复合信号 ，并对进行采样，采样周期为1ms，采用FFT算法分别计算128点、256点、1024点采样信号的频谱，定性分析噪声对信号频谱的影响。 实验程序及结果分析 实验一： 标准方波程序：（） clear; t=-0.016:0.00001:0.016; y=square(125*pi*t); plot(t,y,r); grid on; axis([-0.016 0.016 -1.1 1.1]) title(方波信号（2周期）); xlabel(周期16ms); ylabel(振幅); 标准方波仿真图 白噪声生成程序：（） clear t=-0.016:0.00001:0.016; u1 = (rand(1,length(t)) - 0.5 ) * sqrt( 12 * 0.1 ) + 0; plot(t,u1,r); mean_u1=mean(u1) var_u1=var(u1) sigPower = sum(abs(u1).^2)/length(u1) axis([-0.016 0.016 -1.1 1.1]) grid on; title(白噪声信号); xlabel(时间t); ylabel(振幅); 白噪声仿真图 均值方差信号强度如下图： 生成带有白噪声的方波复合信号程序：（） clear; t=-0.016:0.00001:0.016; y1=square(125*pi*t); u1 = (rand(1,length(t)) - 0.5 ) * sqrt( 12 * 0.05 ) + 0; y=y1+u1; plot(t,y,r); grid on; axis([-0.016 0.016 -1.5 1.5]) title(复合信号（2周期）); xlabel(周期16ms); ylabel(振幅); 复合信号仿真图 标准方波采样程序： clear; fs=1000; N=128; %设置每个周期采样点 t=-N/fs:1/fs:(N-1)/fs; y1=square(125*pi*t); %产生周期16MS的周期信号 u1 = (rand(1,length(t)) - 0.5 ) * sqrt( 12 * 0.1 ) + 0; u1=u1-mean(u1); y=y1+u1; subplot(3,1,1); x=fft(y); n=-N:1:N-1; %m=-0.5:0.5/N:0.5-0.5/N; %f=x/(2*N); stem(n, y1,r,.); grid on; axis([-128 128 -1.5 1.5]) title(复合信号128点采样); xlabel(采样点n); ylabel(振幅); clear; fs=1000; N=256; %设置每个周期采样点 t=-N/fs:1/fs:(N-1)/fs; y1=square(125*pi*t); %产生周期