资源描述：
《数字信号处理fft快速傅里叶变换程序仿真》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、实用文档北京信息科技大学硕士研究生数字信号处理课程仿真实验报告内容仿真结果报告文字和分析格式总分分数学科控制工程姓名崔领袖班级研1507学号20150202072015年月日文案大全实用文档一、实验目的1.通过本次实验熟悉用Matlab产生方波以及白噪声。2.学会运用数字信号处理的知识进行采样，并从采的的点数不同来更加深刻的理解信号采样定理。3.学会运用FFT处理分析信号问题，并充分学习理解频谱在信号中的应用。4.增强独立完成实验的能力。二、实验要求1、首先生成幅值为1，-1，各持续半个周期，周期T=16m
2、s的方波信号；然后生成均值为0，方差为0.1，强度为0.1的白噪声信号。将白噪声信号加在方波信号进而得到复合信号。（1）对进行采样，采样周期为1ms，采用FFT算法分别计算128点、256点、1024点采样信号的频谱。（2）从上述结果中如何估算信号的各次谐波分量，给出分析结果，并与没有噪声情况下的理论分析结果对照。2、首先生成幅值为1，周期为16ms的方波信号，然后生成均值为0，方差分别为0,0.05，0.1，0.2的白噪声信号，将得到方波信号进行计算并与白噪声相互叠加得到复合信号，并对进行采样，采样周期为
3、1ms，采用FFT算法分别计算128点、256点、1024点采样信号的频谱，定性分析噪声对信号频谱的影响。三、实验程序及结果分析实验一：标准方波程序：（）clear;t=-0.016:0.00001:0.016;y=square(125*pi*t);plot(t,y,'r');gridon;axis([-0.0160.016-1.11.1])title('方波信号（2周期）');xlabel('周期16ms');ylabel('振幅');文案大全实用文档标准方波仿真图白噪声生成程序：（）cleart=-0.
4、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.0160.016-1.11.1])gridon;title('白噪声信号');xlabel('时间t');ylabel('振幅');文案大全实用文档白噪声仿真图均值方差信号强度如下图：生成带有白噪声的方波复合信号程
5、序：（）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');gridon;axis([-0.0160.016-1.51.5])title('复合信号（2周期）');xlabel('周期16ms');文案大全实用文档ylabel('振幅');复合信号仿真图标准方波采样程序：clear;fs=1000;N=128;%设置每个周期采样点t=
6、-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','.');gridon;axis([-128128-1.51.5])title('复合信号128点采样')
7、;xlabel('采样点n');ylabel('振幅');clear;fs=1000;N=256;%设置每个周期采样点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,2);x=fft(y);n=-N:1:N-1;%m=-0.5:0.5/N:0.5-0.5/N;%f=x/(2*N
8、);stem(n,y1,'r','.');gridon;axis([-256256-1.51.5])title('复合信号256点采样');xlabel('采样点n');ylabel('振幅');clear;fs=1000;N=1024;%设置每个周期采样点t=-N/fs:1/fs:(N-1)/fs;y1=square(125*pi*t);%产生周期16MS的周期信号u1=(rand(1,length(t))-0
