版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
.7/7信号与系统课设报告课程名称信号与系统课程设计实验实验项目傅里叶变换分析连续时间系统分析离散时间系统分析音频加解扰实验仪器计算机、MATLAB软件系别___光电学院_________专业__________班级/学号___学生姓名_______实验日期成绩_______________________指导教师________一.实验目的:"信号与系统"是一门重要的专业基础课,MATLAB作为信号处理强有力的计算和分析工具是电子信息工程技术人员常用的重要工具之一。本课程设计基于MATLAB完成信号与系统综合设计实验,以提高学生的综合应用知识能力为目标,是"信号与系统"课程在实践教学环节上的必要补充。通过课设综合设计实验,激发学生理论课程学习兴趣,提高分析问题和解决问题的能力。课设内容:傅里叶变换1、实验要求:1.周期信号的谱分析,要求任意给定单频周期信号,能够准确计算出其幅度谱和相位谱,并画出图形,要求正确显示频率。2.非周期信号的频谱分析,要求分析语音信号的幅度谱和相位谱,并画出图形。3.对于方波,设计程序计算其傅里叶级数,仿真吉伯斯现象。2、实验过程1.1周期信号的谱分析1代码clear,clc,closeall%清屏f0=500;%设定正弦的频率N=100;n=0:N;%设置采样点Ts=0.0001;%抽样时间间隔%生成正弦波形y=sin<2*pi*f0.*n*Ts>;subplot<3,2,1>;plot<n*Ts,y>;title<'y=sin<t>'>;xlabel<'s'>;legend<'sin<t>'>;%fftshift的作用正是让正半轴部分和负半轴部分的图像分别关于各自的中心对称。%因为直接用fft得出的数据与频率不是对应的,fftshift可以纠正过来Y=fftshift<fft<y>>;subplot<3,1,2>;stem<<n-N/2>/N/Ts,abs<Y>,'g'>;%画出信号的幅度谱title<'Y的幅度频谱'>;xlabel<'Hz'>subplot<3,1,3>;stem<<n-N/2>/N/Ts,angle<Y>,'g'>;%画出信号的相位谱title<'Y的相位频谱'>;xlabel<'Hz'>2运行结果图像及其分析:图像分析先设置正弦信号频率为500Hz,在生成正弦信号后,对其做傅里叶变换,再将其做集中使之与频率相对应,紧接着将它的频度和相位算出来,画图!1.2非周期信号的频谱分析1代码clear,clc,closeallN=4000;%设定采样点数n=[0:N-1];[x1,fs1,nbits]=wavread<'F:\fmt.wav'>;%进行采样sound<x1,fs1,nbits>;%播放歌曲x=x1<1:N>;%取部分点数X=fft<x>;%傅里叶变换并集中f=n/N*fs1;%把点数转换成频率subplot<3,1,1>;plot<n/fs1,x>;%画出采样的信号ylabel<'原歌曲信号'>;xlabel<'时间/s'>;subplot<3,1,2>;plot<f,abs<X>>;%画出采样信号的幅度谱ylabel<'语音的幅度谱'>;xlabel<'频率/Hz'>;subplot<3,1,3>;plot<f,angle<X>>;%画出采样信号的相位谱ylabel<'语音的相位谱'>;xlabel<'频率/Hz'>;2>运行结果图像及其分析:图形分析用适当的时间间隔对语音信号进行采集,然后将它的幅度和相位计算出来,并画在图上!1.3方波分析,计算其傅里叶级数,仿真吉伯斯现象。1代码closeall;clc;clearall;%清屏t=-pi:0.01:pi;%设置显示的时间范围a=0;T=input<'T='>%由用户指定要T次谐波forn=1:2:Tb=sin<n*t+n*pi/2>/n;%n次谐波a=a+b;endy=4*a/pi;x=square<t+pi/2>;%原先的方波plot<t,y,'b',t,x,'r'>%画出图像1.3运行结果图像N=10N=40N=100N=250程序分析:计算出它的傅里叶系数,然后按指定的要求将它的吉伯斯想象画出来。连续时间系统分析:〔一实验要求:1、自行设计一个有初始条件的微分方程,至少二阶。2、解出其零输入解,并画出图形,与手工计算相比较。3、解出其单位冲激响应,并画出图形,与手工计算相比较。4、设定某一激励信号,用卷积方法解出其零状态响应,并画出图形,与手工计算相比较。5、计算系统的幅频响应和相频响应,并画出图形。〔二程序代码clearall,closeall,clc;A=[0,1;-3,-2];B=[0;1];C=[31];D=[0];sys=ss<A,B,C,D>;%状态方程及输出方程四个矩阵x0=[-1,4];%状态变量初始值,要从初始状态r<0->=1,r'<0->=2计算而得t=[0:0.01:10]';e=1*<t>=0>;ezi=0*<t>=0>;%生成零输入信号的抽样值,e<t>=0rzs=lsim<sys,e,t>;%仿真零状态响应rzi=lsim<sys,ezi,t,x0>;%仿真零输入响应rf=lsim<sys,e,t,x0>;%验证,与手工计算比较rzs_manul=-2*exp<-t>+0.5*exp<-2*t>+1.5;%零状态响应手工计算值error_rzs=rzs-rzs_manul;rzi_manul=4*exp<-t>-3*exp<-2*t>;%零输入响应手工计算值error_rzi=rzi-rzi_manul;rf_manul=2*exp<-t>-2.5*exp<-2*t>+1.5;%全响应手工计算值error_rf=rf-rf_manul;%画图figuresubplot<2,2,1>holdon,boxon;title<'系统的零输入响应、零状态响应、完全响应'>set<gca,'FontSize',10>;plot<t,e,'b.','linewidth',3>;plot<t,rzs,'r-','linewidth',3>;plot<t,rzi,'g-.','linewidth',3>;plot<t,rf,'y','linewidth',3>;legend<'input','rzs','rzi','rf'>;xlabel<'t/s'>;axis<[-0.010-14]>%画手工得到的响应subplot<2,2,2>;plot<t,rzs_manul,'r:','linewidth',3>;holdon;plot<t,rzi_manul,'g-.','linewidth',3>;plot<t,rf_manul,'k','linewidth',3>;legend<'零状态响应','零输入响应','全响应'>;title<'手工得到的响应'>xlabel<'时间/s'>;axis<[-0.010-14]>〔三运行结果图像及其分析系统的零输入响应、零状态相应、完全响应和手工计算的零输入响应、零状态相应、完全响应,以及单位冲激响应和卷积之后的零状态响应图像分析:将系统的零输入响应、零状态响应和完全响应计算出来,并与手工相比较,差距不大。并画出来,再将系统的冲击响应求出来,用将激励响应和冲激响应做卷积,把系统的零状态响应求出来,画在图上!系统的幅频响应和相频响应系统的幅度在500Hz时有一个冲击,系统的相位在500Hz时发生翻转。离散时间系统分析:〔一实验要求:1.自行设计某离散时间系统函数,至少是二阶的系统,画出零极点图,判断系统的稳定性。2.求出单位样值响应,并画出图像。3.求出系统的幅频响应和相频响应,并画出图形。〔二程序代码clearall,closeall,clc;%清屏%建立一个二阶系统%系统的差分方程为y〔n-5y<n-1>+0.6y<n-2>=5x<n>+2x<n-1>+x<n-2>den=[1,5,0.6];%定义左侧系数num=[5,2,1];%定义右侧系数t=[-5;0.01;5]';%设置时间范围sys=tf<num,den>;poles=roots<den>%计算多项式的根figure;subplot<2,1,1>%pzmap<sys>;zplane<num,den,'r'>;%画出sys所描述系统的零极点图h=filter<num,den,t>;%计算样值响应subplot<2,1,2>%plot<h>;impz<num,den>;%画出系统的样值响应figure,freqz<num,den>;%计算并画出系统的幅频响应和相频响应〔三运行结果图像及其分析零极点图,单位样值响应建立一个二阶系统,并将它的零极点求出,极点并没有全部在单位圆内,零点都在单位圆内;求出系统的样值响应;画出图像。在将系统的幅频响应和相频响应做出来。系统的幅频相频这是一个带阻滤波器,因为它在零时刻为高,之后为低,然后再为高!音频加解扰:实验要求:利用AM调制解调的基本原理和滤波算法实现音频扰频和解扰功能。首先对信号频率进行翻转,即将高频信号反转到低频处,将低频信号反转到高频处;然后再通过一个0.2-3.4kHz的带通滤波器进行滤波,其输出信号就是扰频信号。当扰频信号通过另一个运行同样程序的系统时,该信号就被解扰,恢复为原始可以被接听的语音信号。1、编写扰频解扰程序,对于语音输入,画出语音信号的频谱、滤波器的响应曲线、扰频信号的频谱,要求扰频输出时能用耳机听到干扰声音。2、利用扰频解扰程序进行解扰,画出解扰后信号的频谱,要求解扰后能听到恢复的声音。〔二代码clear,clc,closeall%清屏N=40000;%采样点数n=[0:N-1];[x1,fs1,nbits]=wavread<'F:\fmt.wav'>;%采语音信号%sound<x1,fs1,nbits>;%播放歌曲x=x1<1:N>;Ts=1/fs1;%抽样时间间隔X=fft<x>;%傅里叶变换f=n/N*fs1;%把点数转换成频率b=abs<X>;subplot<2,1,1>;plot<f,b>;title<'语音的幅度谱'>;xlabel<'频率/Hz'>;%翻转频谱f1=4000;y=x'.*exp<j*2*pi*f1.*n*Ts>;%进行翻转Y=fft<y>;subplot<2,1,2>;plot<f,abs<Y>>;%画出翻转后的幅度图title<'翻转后语音的幅度谱'>;xlabel<'频率/
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024年承台建设合同书规范文本
- 2024年建筑施工承揽协议模板
- 小学开展的活动方案及策划8篇
- 2024年房地产代理销售协议
- 2024年工程劳务合作合同
- 2024年建筑抹灰工程分包合同
- 2024年招投标管理第六版
- 2024年影视明星经纪服务合同
- 2024年房屋装修解除合同协议书
- 2024年式全屋定制合同样本
- 液化石油气充装操作规程(YSP118液化石油气钢瓶)
- 工程样板过程验收单
- 颅内动脉动脉瘤介入治疗临床路径
- 粮食仓储场建设项目可行性研究报告
- 珠宝销货登记表Excel模板
- 深基坑开挖施工风险源辨识与评价及应对措施
- 唯美手绘风花艺插花基础培训PPT模板课件
- 《现代汉语语法》PPT课件(完整版)
- 5G智慧农业建设方案
- 航海学天文定位第四篇天文航海第1、2章
- 浙江大学学生社团手册(08)
评论
0/150
提交评论