号与系统课设报告.docVIP

课程设计目的 “信号与系统”是一门重要的专业基础课，MATLAB作为信号处理强有力的计算和分析工具是电子信息工程技术人员常用的重要工具之一。本课程设计基于MATLAB完成信号与系统综合设计实验，以提高学生的综合应用知识能力为目标，是“信号与系统”课程在实践教学环节上的必要补充。通过课设综合设计实验，激发学生理论课程学习兴趣，提高分析问题和解决问题的能力。 课程设计时间 第十五、十六周。上机时间安排见附件一。第十六周周五提交课程设计报告并答辩。 参考书目 谷源涛、应启珩、郑君里著，信号与系统——MATLAB综合实验，北京：高等教育出版社，2008年1月。 郑君里、应启珩 、杨为理，信号与系统引论，北京：高等教育出版社，2009年3月。 梁虹等，信号与系统分析及Matlab实现，北京：电子工业出版社，2002年2月。 注意事项 基本部分，共三道题，每人都需要全部完成，要求十五周周五做完。 提高部分，共八道题，每人按照学号分配（见附件二）只做其中的一题。 第十六周周五所提交的课程设计报告如有雷同，一律退回重写。 一、基础部分 （一） 傅里叶变换分析 1、周期信号的谱分析，要求任意给定单频周期信号，能够准确计算出其幅度谱和相位谱，并画出图形，要求正确显示频率。 %周期信号频谱分析 t=0:0.2:8; x=cos(t); %产生单频周期信号 subplot(311); plot(x); xlabel(t); ylabel(cos(t)); title(波形); X=fft(x); %计算DFS系数 a=abs(X); b=angle(X); subplot(312); stem(a); %画出幅度频谱 title(幅度谱); xlabel(w); ylabel(a); subplot(313); stem(b); %画出相位频谱 title(相位谱); xlabel(w); ylabel(b); 图1 图1中是cos函数在一个周期内的波形、幅度谱和相位谱，由图可知余弦函数傅氏变换的幅度谱是两个冲击，相位谱为一条直线。 非周期信号的频谱分析，要求分析语音信号的幅度谱和相位谱，并画出图形。 %非周期信号频谱分析 y=wavread(u);%读取音乐文件 Fn=fft(y);%作傅里叶变换 f=44000; Fnf=abs(Fn); %幅度 Fnx=angle(Fn);%相位 N=length(y); M=0:N-1; figure; subplot(3,1,1); plot(M/7005000,y);%语音信号波形 title(语音信号的波形); xlabel(t); ylabel(y); subplot(3,1,2); stem(M/f,Fnf);%幅度谱 axis([0 0.0005 0 1500]); title(幅度谱); xlabel(Hz); ylabel(Fnf); subplot(3,1,3); stem(M/f,Fnx);%相位谱 axis([0 0.0005 0 10]); title(相位谱); xlabel(Hz); ylabel(Fnx); 图2 由图2可以看出，语音信号的幅度谱变化不大，相位谱是断续的。 3、对于方波，设计程序计算其傅里叶级数系数，仿真吉伯斯现象。 %仿真吉伯斯现象 t=0:0.01:20*pi;%定义时间长度 y=0; num=input(叠加的正弦波个数=) for k=1:num %用for循环实现波形的叠加 y=y+(4/pi)*(1/((k*2)-1))*sin(((k*2)-1)*t); r(k,:)=y; end figure; plot(t,r); title(吉伯斯现象); xlabel(t); ylabel(f(t)); grid on; 图3 由图3可以发现，N个正弦波的叠加后的波形近似于方波。 （二）、连续时间系统分析： 1、自行设计一个有初始条件的微分方程，至少二阶。 2、解出其零输入解，并画出图形，与手工计算相比较。 3、解出其单位冲激响应，并画出图形，与手工计算相比较。 4、设定某一激励信号，用卷积方法解出其零状态响应，并画出图形，与手工计算相比较。 5、计算系统的幅频响应和相频响应，并画出图形。 a=[1,3,2]; %微分方程式左边的系数 b=[1,3]; %微分方程式右边系数 sys=tf(b,a); %系统的传递函数 rsys=ss(sys); %讲传递函数转换为状态空间 t=0:0.01:50; %定义时间常数 x=0*t; e=[1,2]; rzi=lsim(rsys,x,t,e); %零输入解 subplot