频率Hzylabel幅度title幅值谱grido
subplot313plotfa
gley1绘出Nyquist频率之前随频率变化的相位xlabel频率Hzylabel相位title相位谱grido
结果图：
结果分析：cos函数波形为周期信号，其频域响应为两个冲击函数，并且符合对偶性。
（2）傅里叶级数分析分析周期三角波的傅里叶级数系数，用正弦信号的线性组合构成三角波，要求谐波次数可以任意输入，分析不同谐波次数所构成的三角波，解释是否存在吉伯斯现象。
设计思路：根据原理任意周期信号都可以表达成傅里叶级数的形式，对周期三角波进行傅里叶级数分解，利用for循环完成级数求和运算。当N值较小时，傅里叶级数的逼近效果不是很理想，随着N值变大，傅里叶级数越来越接近理想值。
源程序：
clcclearcloseall
13
ft6610006Ni
putpleasei
putthe
u
berw0piXNzeros1le
gthttao000000001for
1N
f
4si
pi22
tao2pi2XNXNf
cos
w0te
dXNXN05figureplottXNtitleGibbsN
um2strNxlabelTimesecylabelX
um2strNt结果图；
结果分析：随着N值的增大，傅里叶级数与理想的周期三角波越来越接近，符合预期效果。（3）系统分析任意给定微分方程或差分方程描述的系统，画出系统的幅频响应和相频响应。设计思路：根据微分方程、差分方程与系统函数的对应关系，结合matlab自带的freqs和freqz两个函数，分析系统的幅频特性和相频特性。源程序：clcclearcloseall微分方程b10y的系数a132x的系数
13
ffigurefreqsba差分方程b01y的系数a122x的系数figurefreqzba结果图；
结果分析：微分方程所代表的系统具有高通特性，差分方程所代表的系统具有低通特性。（4）音乐合成程序设计对于任意一小段音乐，利用“十二平均律”计算该音乐中各个乐音的频率，产生并播放这些乐音。分析音乐的频谱，从中识别出不同的乐音。设计思路：根据时间长短来区别各个音符拍长短，根据频率高低来区别各个音符的音调，具体细节参考“十二平均律”，将每段乐音连接起来，使用sou
d函数播放乐音。对每段乐音进行傅
13
f里叶变换，分析其频谱。
源程序：
clcclearcloseall
f8000音乐采样频率，可改为4000或者16000t201f12拍，时间长短不同t401f051拍t801f02512拍
misc_
ote552325不同音符频率misc_
ote658733misc_
ote2392misc_
ote134923misc_
ote6_dow
29366
m1si
2pimisc_
ote5t4波形m2si
2pimisc_
ote5t8m3si
2pimisc_
ote6t8m4si
2pimisc_
ote2t2m5si
2pimisc_
r