资源描述:
《实验三信 号 频 谱 分 析 》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、南京工程学院通信工程学院实验报告题目:信号频谱分析课程名称信号与系统实验专业电子信息工程班级学生姓名学号设计地点指导教师实验时间:2014年月日实验三:信号频谱分析一、实验目的:1、掌握傅立叶级数(FT),学会分析连续时间周期信号的频谱及MATLAB实现;2、掌握傅立叶变换F(jw),了解傅立叶变换的性质以及MATLAB实现。3、掌握信号抽样与恢复的原理,能够用MATLAB实现一般信号的采样与恢复。二、课内验证与设计实验1、画出书中P121且周期为2的方波图形,对其傅立叶级数(正弦)取不同的N(N可自己随意取)值,画出其近似波
2、形,并注意吉伯斯现象。程序如下:t=0:0.01:4;x=zeros(10,max(size(t)));y=zeros(10,max(size(t)));fork=1:2:9;x1=sin(pi*k*t)/k;x(k,:)=x(k,:)+x1;y((k+1)/2,:)=x(k,:);endsubplot(221)plot(t,y(1:9,:))subplot(222)N=5;ft=0;forn=1:N;ft=ft+(4/pi)*(sin((2*n-1)*pi*t)/(2*n-1));endplot(t,ft);subplot(
3、223)N=10;ft=0;forn=1:N;ft=ft+(4/pi)*(sin((2*n-1)*pi*t)/(2*n-1));endplot(t,ft);subplot(224)N=50;ft=0;forn=1:N;ft=ft+(4/pi)*(sin((2*n-1)*pi*t)/(2*n-1));endplot(t,ft);2、利用数值法定义求门宽为2的门函数的傅氏变换,画出频谱图,与前面3中的(3)进行比较并对此信号进行移时与移频,观察频谱的变化。程序如下:t=-3:0.02:3;f=u(t+1)-u(t-1);w1=4*
4、pi;k=0:500;w=k*w1/500;subplot(441);plot(t,f)axis([-3,3,-0.5,2]);title(‘f(t)=u(t+1)-u(t-1)’);subplot(442);plot(w,F);title(‘F(w)’);subplot(443);F1=abs(F);plot(w,F1);subplot(444);P1=angle(F);plot(w,P1*180/pi);f=u((t-1)+1)-u((t-1)-1);w1=4*pi;k=0:500;w=k*w1/500;F=f*exp(-
5、1i*t’*w)*0.02;F=real(F);w=[-fliplr(w),w(2:501)];F=[fliplr(F),F(2:501)];subplot(445);plot(t,f)axis([-3,3,-0.5,2]);title(‘f(t-1)’);subplot(446);plot(w,F);title(‘f(t-1)-->F(w)’);subplot(447);F1=abs(F);plot(w,F1);subplot(448);P1=angle(F);plot(w,P1*180/pi);f=u(t+1)-u(t-1
6、);f=f1.*exp(li*5*t);w1=4*pi;k=0:500;w=k*w1/500;F=f*exp(-1i*t’*w)*0.02;F=real(F);subplot(449);plot(t,f)F=f*exp(-1i*t’*w)*0.02;F=real(F);w=[-fliplr(w),w(2:501)];F=[fliplr(F),F(2:501)];subplot(445);axis([-3,3,-0.5,2]);title(‘f(t)*exp(5jt)’);subplot(4,4,10)plot(w,F);tit
7、le(‘f(t)*exp(5jt)-->F(w)’);axis([-20,20,-2,2]);subplot(4,4,11);F1=abs(F);plot(w,F1);axis([-20,20,0,2]);subplot(4,4,12);P1=angle(F);plot(w,P1*180/pi);axis([-20,20,0,200]);f1=u(t+1)-u(t-1);f=f1.*exp(-li*5*t);w1=4*pi;k=0:500;w=k*w1/500;F=f*exp(-1i*t’*w)*0.02;F=real(F);
8、subplot(4,4,13);plot(t,f)axis([-3,3,-0.5,2]);title(‘f(t)*exp(-5jt)’)subplot(4,4,14);plot(w,F);title(‘f(t)*exp(5jt)-->F(w)’);axis([-20,20,-