常用信号的频谱分析及时域采样定理.docxVIP

开课学期 2016-2017 学年第 2 学期实验课程信号与系统仿真实验实验项目常用信号的频谱分析及时域采样定理班级学号学生姓名实验时间实验台号A11操作成绩报告成绩实验目的掌握常用信号的频域分析方法；掌握时域采样定理；掌握时域采样信号恢复为原来连续信号的方法及过程。二、实验性质验证性预习内容时域采样定理的内容及信号时域采样过程；连续信号经时域采样后，信号的频谱发生的变化；时域采样信号恢复为原来连续信号的方法及过程。实验内容(编写程序，绘制实验结果)实现周期信号的频谱f(t)=sin(*80t)程序：fa=sin(2.*pi.*80.*t);%原信号fs0=10000; %采样频率tp=0.1;%时间范围t=[-tp:1/fs0:tp];%信号持续时间范围k1=0:999;k2=-999:-1;m1=length(k1);m2=length(k2);f=[fs0*k2/m2,fs0*k1/m1];%信号频率范围w=[-2*pi*k2/m2,2*pi*k1/m1];fx1=eval(fa);%把文本fa赋值给信号fx1FX1=fx1*exp(-j*[1:length(fx1)]*w);%进行傅立叶变换figuresubplot(2,1,1),plot(t,fx1,r);title(原信号);xlabel(时间t(s));%原信号的时域波形图axis([min(t),max(t),min(fx1),max(fx1)]);subplot(212),plot(f,abs(FX1),r);title(原信号频谱);xlabel (频率f(Hz));%频域波形图axis([-100,100,0,max(abs(FX1))+5]);实现非周期信号的频谱,要求记录结果并对结果进行分析讨论.(1)门函数信号的频谱分析,(2)尺度变换之后门函数的频谱分析.程序：令tao=1syms tx=heaviside(t+0.5)-heaviside(t-0.5);F=fourier(x);subplot(211);ezplot(x,[-2,2]);subplot(212);ezplot(F,[-10,10]);程序：令tao=1，a=4syms tx=heaviside(t+(1/8))-heaviside(t-(1/8));F=fourier(x);subplot(211);ezplot(x,[-2,2]);axis([-2,2,-1,2])subplot(212);ezplot(F); axis([-5,5,-0.5,0.5]);分析：经过尺度变换，门函数的时间常数tao改变了，tao从1变成了1/4，门函数的幅度保持不变，但频谱变化幅度比尺度变换前缓慢，频谱的基波分量降低了时域采样及其恢复运行给定实验程序，绘制运行实验结果，总结实验结果，说明采样过程及恢复原信号的原理。程序：syms t w f; %定义符号变量f=(1-2*abs(t))*exp(-j*w*t); %计算被积函数F=int(f, t, -1/2, 1/2); %计算傅里叶系数F(w)F=simple(F);F %化简subplot(3, 1, 1), %绘制三角波的幅频特性曲线F(w)low=-26*pi;high=-low; %设置w的上界和下界ezplot(abs(F), [low:0.01:high]); axis([low high -0.1 0.5]); xlabel();title(三角波的频谱); subplot(3, 1, 2), %绘制经过截止频率为4*pi低通滤波器后的频谱Y1(w)ezplot(abs(F), [-4*pi:0.01:4*pi]);axis([low high -0.1 0.5]);title(低通滤波后的频谱);%采样信号的频谱是原信号频谱的周期延拓，延拓周期为(2*pi)/Ts%利用频移特性F[f(t)*exp(-j*w0*t)]=F(w+w0)来实现subplot(3, 1, 3); %绘制采样后的频谱Y(w)Ts=0.2; %采样信号的周期w0=(2*pi)/Ts; %延拓周期10*pifor k=-2:2 ft=f*exp(-j*w0*k*t); FT=int(ft, t, -1/2, 1/2); ezplot((1/Ts)*abs(FT), [(-4*p