资源描述：
《信号与系统课程设计实验-傅里叶变换等x》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、实验报告课程名称信号与系统课程设计实验实验项目傅里叶变换实验解微分方程实验离散时间系统函数实验码分复用实验仪器计算机，MATLAB软件系别___光电学院_________专业____电子信息工程_______班级/学号__电信0905学生姓名________实验日期2011.6成绩_______________________指导教师_____焦老师_______一．实验目的：“信号与系统”是一门重要的专业基础课，MATLAB作为信号处理强有力的计算和分析工具是电子信息工程技术人员常用的重要工具之一。本课程设计基于MATLAB完成信号与系统综合设计实验，以提高学生的综合应用
2、知识能力为目标，是“信号与系统”课程在实践教学环节上的必要补充。通过课设综合设计实验，激发学生理论课程学习兴趣，提高分析问题和解决问题的能力。二、课设内容：1、傅里叶变换分析：1).周期信号的谱分析，要求任意给定单频周期信号，能够准确计算出其幅度谱和相位谱，并画出图形，要求正确显示频率。源程序：clc;clear;closeall;%清屏f=110;%信号频率fs=10*f;%采样频率N=fs/f;n=0:N;subplot(3,1,1);x=cos(2*pi*f*n/fs);%建立余弦波形plot(n,x);x=fft(x);%傅里叶变换x=fftshift(x);%移位
3、subplot(3,1,2);stem((n-N/2)/N*f,abs(x));%离散型输出频谱图ylabel('Cn的幅度');subplot(3,1,3);stem((n-N/2)/N*f,angle(x));%相谱图ylabel('Cn的相位');xlabel('omega/omega0');截图：2).非周期信号的频谱分析，要求分析语音信号的幅度谱和相位谱，并画出图形。n=[0:N-1];[ss,fs,bits]=wavread('E:matlab1worksuyan(1).wav');%加载文件s=ss(1:N);f1=fft(s);%傅里叶变换f=n/
4、N*fs;%把点数转换成频率subplot(3,1,1);plot(ss);title('原信号');subplot(3,1,2)plot(f,abs(f1))title('幅度谱');xlabel('频率/Hz');subplot(3,1,3)plot(f,angle(f1))title('相位谱');xlabel('频率/Hz');sp1=max(abs(s));sound(s/sp1,fs);3)、对于方波，设计程序计算其傅里叶级数系数，仿真吉伯斯现象。clear,clc,closeall;%清屏N=500000;t=-2:0.001:2;T=1;ft1=rectpu
5、ls(t,T);ft2=rectpuls(t-2,T);ft3=rectpuls(t+2,T);ft=ft1+ft2+ft3;N=input('Numberofharmonics=');c0=0.5;xN=c0*ones(1,length(t));forn=1:NxN=xN+cos(pi*n*t)*sinc(n/2);%正弦分量为零endplot(t,ft,t,xN);Numberofharmonics=5Numberofharmonics=50Numberofharmonics=10Numberofharmonics=5002、连续时间系统分析：1）、自行设计一个有初始条
6、件的微分方程，至少二阶。2）、解出其零输入解，并画出图形，与手工计算相比较。3）、解出其单位冲激响应，并画出图形，与手工计算相比较。4）、设定某一激励信号，用卷积方法解出其零状态响应，并画出图形，与手工计算相比较。5）、计算系统的幅频响应和相频响应，并画出图形。基本题2：解微分方程，习题2-6clearall,closeall,clc;%b=[1,3],a=[132],初始状态r(0-)=1,r'(0-)=2%先由系统函数转成状态方程及输出方程源程序：A=[0,1;-2,-3];B=[0;1];C=[31];D=[0];b=[1,3];a=[132];num=[13];de
7、n=[132];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);%单位冲击响应rzs=lsim(sys,e,t);%零状态响应rzi=lsim(sys,ezi,t,x0);%零输入响应rf=lsim(sys,e,t,x0);%完全响应h=impulse(sys,t);%冲击响应[H,w]=freqs(num,
