连续离散系统频域分析

课程实验报告学年学期2015-2016学年第二学期课程名称信号与系统实验名称连续和离散系统的频域分析实验室北校区5号楼计算机房专业年级电气141学生姓名宋天绍学生学号2014011595提交时间2016.6.19成绩任课教师吴凤娇水利与建筑工程学院实验二：连续和离散系统的频域分析一：实验目的1：学习傅里叶正变换和逆变换，理解频谱图形的物理含义2：了解连续和离散时间系统的单位脉冲响应3：掌握连续时间系统的频率特性二：实验原理1.傅里叶正变换和逆变换公式正变换：()()jtFftedt逆变换：1()()2jtftFed2.频域分析tjtjeddete)(21)(21)(将激励信号分解为无穷多个正弦分量的和。detrtjzs)()(21)(，R(ω)为)(trzs傅里叶变换;d)(各频率分量的复数振幅激励单位冲激响应时的零状态响应)(t)(th单位阶跃响应时的零状态响应激励)(tu)(tg3各函数说明：(1)impulse冲激响应函数：[Y,X,T]=impulse(num,den);)1()2()1()1()2()1()()()(11nasasambsbsbsAsBsHnnmmnum分子多项式系数；num=[b(1)b(2)…b(n+1)];den分母多项式系数；den=[a(1)a(2)…a(n+1)];Y,X,T分别表示输出响应，中间状态变量和时间变量；如：352)(2ssssH，等价于)(2)()(3)(5)(tetetrtrtr定义den=[153];num=[12];[Y,X,T]=impulse(num,den);(2)step阶跃响应函数：[Y,X,T]=step(num,den);num分子多项式；den分母多项式Y,X,T分别表示输出响应，中间状态变量和时间变量；如：352)(2ssssH，den=[153];num=[12];[Y,X,T]=step(num,den);（3）impz数字滤波器的冲激响应[h,t]=impz(b,a,n)b分子多项式系数；a分母多项式系数；n采样样本h离散系统冲激响应；t冲激时间，其中t=[0:n-1]',n=length(t)时间样本数（4）freqs频域响应[h,w]=freqs(b,a,f)b,a定义同上，f频率点个数h频域响应，w频域变量)1()2()1()1()2()1()()()(11masasanbsbsbsAsBsHmmnn三．实验内容1．周期信号傅里叶级数已知连续时间信号()()2/π8cos3/π4coscos)(321++++=tAtAtAtx，其中321,,AAA取值如下：（X为学号的后两位）]10,1[,5.02321XXAXAXA]20,11[,55321XXAXAXAXAXAXA32151020,X要求画出信号的时域波形和频域波形（幅度谱和相位谱）。分析该信号有几个频率成分，频率分别是多少，振幅为多少，相位为多大。理解并体会连续信号可以分解为无穷多正弦波叠加。(1)Commandwindow程序清单：%%信号的频域成分表示法例子：正弦波的叠加t=0:20/400:20;w1=1;w2=4;w3=8;fai1=0;fai2=pi/3;fai3=pi/2;%在命令窗口分别输入A1，A2，A3振幅值A1=input('InputtheamplitudeA1forw1=1:');A2=input('InputtheamplitudeA2forw2=4:');A3=input('InputtheamplitudeA3forw3=8:');%连续时间信号形x(t)f1=A1*cos(w1*t+fai1);f2=A2*cos(w2*t+fai2);f3=A3*cos(w3*t+fai3);x=A1*cos(w1*t+fai1)+A2*cos(w2*t+fai2)+A3*cos(w3*t+fai3);figure(1);subplot(211),plot(t,f1,'r',t,f2,'g',t,f3,'b','linewidth',4)title('连续时间信号时域图形x(t)')ylabel('x(t)')xlabel('时间（秒）')legend({'f1=A1*cos(w1*t+fai1)','f2=A2*cos(w2*t+fai2)','f3=A3*cos(w3*t+fai3)'})subplot(212),plot(t,x,'linewidth',4)title('连续时间信号时域图形x(t)')ylabel('x(t)')xlabel('时间（秒）')figure(2)subplot(211),stem([w1w2w3],[A1A2A3])v=[01001.5*max([A1,A2,A3])];axis(v);%限定XY轴坐标范围title('幅频特性')ylabel('振幅')xlabel('频率(弧度/秒)')subplot(212),stem([w1w2w3],2*pi*[fai1fai2fai3])fai=[01001.5*max(2*pi*[fai1fai2fai3])];axis(fai);%限定XY轴坐标范围title('相频特性')ylabel('相位(度)')xlabel('频率(弧度/秒)')(2)M文件函数清单functionFS(w1,w2,w3,A1,A2,A3)%调用格式：FS(1,4,8,10,5,30)%信号的频域成分表示法%例子：正弦波的叠加t=0:20/400:20;fai1=0;fai2=pi/3;fai3=pi/2;%连续时间信号形x(t)f1=A1*cos(w1*t+fai1);f2=A2*cos(w2*t+fai2);f3=A3*cos(w3*t+fai3);x=A1*cos(w1*t+fai1)+A2*cos(w2*t+fai2)+A3*cos(w3*t+fai3);figure(1);subplot(211)plot(t,f1,'r',t,f2,'g',t,f3,'b','linewidth',2)title('连续时间信号时域图形x(t)')xlabel('时间（秒）')ylabel('x(t)')legend({'f1=A1*cos(w1*t+fai1)','f2=A2*cos(w2*t+fai2)','f3=A3*cos(w3*t+fai3)'})subplot(212)plot(t,x,'linewidth',4)title('连续时间信号时域图形x(t)')xlabel('时间（秒）')ylabel('x(t)')figure(2)subplot(211)stem([w1w2w3],[A1A2A3])v=[01001.5*max([A1,A2,A3])];axis(v);%限定XY轴坐标范围title('幅频特性')xlabel('频率(弧度/秒)')ylabel('振幅')subplot(212)stem([w1w2w3],2*pi*[fai1fai2fai3])fai=[01001.5*max(2*pi*[fai1fai2fai3])];axis(fai);%限定XY轴坐标范围title('相频特性')xlabel('频率(弧度/秒)')ylabel('相位(度)')实验代码及过程：%%信号的频域成分表示法例子：正弦波的叠加t=0:20/400:20;w1=1;w2=4;w3=8;fai1=0;fai2=pi/3;fai3=pi/2;%在命令窗口分别输入A1，A2，A3振幅值A1=input('InputtheamplitudeA1forw1=1:');A2=input('InputtheamplitudeA2forw2=4:');A3=input('InputtheamplitudeA3forw3=8:');%连续时间信号形x(t)f1=A1*cos(w1*t+fai1);f2=A2*cos(w2*t+fai2);f3=A3*cos(w3*t+fai3);x=A1*cos(w1*t+fai1)+A2*cos(w2*t+fai2)+A3*cos(w3*t+fai3);figure(1);subplot(211),plot(t,f1,'r',t,f2,'g',t,f3,'b','linewidth',4)title('连续时间信号时域图形x(t)')ylabel('x(t)')xlabel('时间（秒）')legend({'f1=A1*cos(w1*t+fai1)','f2=A2*cos(w2*t+fai2)','f3=A3*cos(w3*t+fai3)'})subplot(212),plot(t,x,'linewidth',4)title('连续时间信号时域图形x(t)')ylabel('x(t)')xlabel('时间（秒）')figure(2)subplot(211),stem([w1w2w3],[A1A2A3])v=[01001.5*max([A1,A2,A3])];axis(v);%限定XY轴坐标范围title('幅频特性')ylabel('振幅')xlabel('频率(弧度/秒)')subplot(212),stem([w1w2w3],2*pi*[fai1fai2fai3])fai=[01001.5*max(2*pi*[fai1fai2fai3])];axis(fai);%限定XY轴坐标范围title('相频特性')ylabel('相位(度)')xlabel('频率(弧度/秒)')在弹出的命令行输入数值：InputtheamplitudeA1forw1=1:75——————学号85-10InputtheamplitudeA2forw2=4:80——————学号85-5InputtheamplitudeA3forw3=8:85——————学号85实验结果：2傅里叶的正变换和逆变换dxexfwFxj)()(调用符号工具箱中F=fourier(f)函数返回傅里叶变换F(w)f=ifourier(F)函数返回被积函数f(t)02468101214161820-100-50050100连续时间信号时域图形x(t)x(t)时间（秒）f1=A1*cos(w1*t+fai1)f2=A2*cos(w2*t+fai2)f3=A3*cos(w3*t+fai3)02468101214161820-400-2000200400连续时间信号时域图形x(t)x(t)时间（秒）012345678910050100幅频特性振幅频率(弧度/秒)0123456789100510相频特性相位(度)频率(弧度/秒)(1)分别求)100sin()(ttf,()()()22ftEutut对应的傅里叶变换程序清单：%%矩形脉冲的傅里叶变换symstt0EFwtauff=E*(heaviside(t-tau/2)-heaviside(t+tau/2));%heaviside单位阶跃信号Fw=fourier(f);simplify(Fw)%简化函数计算过程，结果中的dirac是单位冲击信号%%正弦信号的傅里叶正变换symstwfFwf=A1*sin(100*pi*t);Fw1=simplify(fourier(f))%fourier正变化函数，返回值频域F(w)结果是：Fw1=pi*(dirac(-100*pi-w)-dirac(100*pi-w))*i(2)分别求)2()(wsaEwF，)()(0的原函数)(tf%%傅里叶逆变换SymswFtfrealE=1;tau=2;F=E*tau*sinc(w*tau/(2*pi));%定义F(w)f=ifourier(F);%傅里叶逆变换函数f=simple