课程设计报告_郑一雄2.docVIP

课 程 设 计 报 告 课程名称 信号与系统课程设计 指导教师 赵老师 设计起止日期 202012-6-22 学 院 信息与通信工程专 业 电子信息工程学生姓名 班级/学号 成 绩 目 录 1、课程设计目的 1 2、课程设计要求 1 3、课程设计任务 1 4、课程设计内容 1 5、总结 11 参考文献 12 附录 12 1、课程设计目的 clc; clear; close all; Fs =128; % 采样频率 T = 1/Fs; % 采样周期 N = 600; % 采样点数 t = (0:N-1)*T; % 时间，单位：S x=2*cos(5*2*pi*t); n=0:N-1; figure; subplot(3,1,1) plot(t,x); xlabel(时间/S); ylabel(振幅);title(时域波形);grid on; y=fft(x,N); y1=fftshift(y); n1=-(N-1)/2:(N-1)/2; f=n1*Fs/N; subplot(3,1,2) plot(f,10*log10(abs(y1))); % 绘出Nyquist频率之前随频率变化的振幅 xlabel(频率/Hz); ylabel(幅度);title(幅值谱);grid on; subplot(3,1,3) plot(f,angle(y1)); % 绘出Nyquist频率之前随频率变化的相位 xlabel(频率/Hz); ylabel(相位);title(相位谱);grid on; 结果图： 结果分析： cos函数波形为周期信号，其频域响应为两个冲击函数，并且符合对偶性。 （2）傅里叶级数分析 分析周期三角波的傅里叶级数系数，用正弦信号的线性组合构成三角波，要求谐波次数可以任意输入，分析不同谐波次数所构成的三角波，解释是否存在吉伯斯现象。 设计思路： 根据原理——任意周期信号都可以表达成傅里叶级数的形式，对周期三角波进行傅里叶级数分解，利用for循环完成级数求和运算。当N值较小时，傅里叶级数的逼近效果不是很理想，随着N值变大，傅里叶级数越来越接近理想值。 源程序： clc; clear; close all; t=-6:6/1000:6; N = input(please input the nunber:); w0=pi; XN=zeros(1,length(t)); tao=0; % 0.0000001; for n=1:N; fn=-4*(sin(n*pi/2))^2/((n+tao)^2*pi^2); XN=XN+fn*cos(n*w0*t); end XN=XN+0.5; figure; plot(t,XN); title([Gibbs,N=,num2str(N)]); xlabel(Time(sec)); ylabel([X,num2str(N),(t)]); 结果图； 结果分析： 随着N值的增大，傅里叶级数与理想的周期三角波越来越接近，符合预期效果。 （3）系统分析 任意给定微分方程或差分方程描述的系统，画出系统的幅频响应和相频响应。 设计思路： 根据微分方程、差分方程与系统函数的对应关系，结合matlab自带的freqs和freqz两个函数，分析系统的幅频特性和相频特性。 源程序： clc; clear; close all; % --- 微分方程 b = [1 0]; % y的系数 a = [1 3 2]; % x的系数 figure; freqs(b,a); % --- 差分方程 b = [0 1]; % y的系数 a = [1 -2 2]; % x的系数 figure; freqz(b,a); 结果图； 结果分析： 微分方程所代表的系统具有高通特性，差分方程所代表的系统具有低通特性。 （4）音乐合成程序设计 对于任意一小段音乐，利用“十二平均律”计算该音乐中各个乐音的频率，产生并播放这些乐音。分析音乐的频谱，从中识别出不同的乐音。 设计思路： 根据时间长短来区别各个音符拍长短，根据频率高低来区别各个音符的音调，具体细节参考“十二平均律”，将每段乐音连接起来，使用sound函数播放乐音。对每段乐音进行傅里叶变换，分析其频谱。 源程序： clc; clear; close all; f =8000; % 音乐采样频率，可改为4000或者16000 t2=0:1/f:1; % 2拍，时间长短不同 t4=0