提高题音乐合成信号与系统课程设计报告学位论文文档格式.docx
《提高题音乐合成信号与系统课程设计报告学位论文文档格式.docx》由会员分享，可在线阅读，更多相关《提高题音乐合成信号与系统课程设计报告学位论文文档格式.docx（19页珍藏版）》请在冰豆网上搜索。
“信号与系统”是一门重要的专业基础课，MATLAB作为信号处理强有力的计算和分析工具是电子信息工程技术人员常用的重要工具之一。
本课程设计基于MATLAB完成信号与系统综合设计实验，以提高学生的综合应用知识能力为目标，是“信号与系统”课程在实践教学环节上的必要补充。
通过课设综合设计实验，激发学生理论课程学习兴趣，提高分析问题和解决问题的能力。
要求：
（1）运用MATLAB编程得到简单信号、简单信号运算、复杂信号的频域响应图；
（2）通过对线性时不变系统的输入、输出信号的时域和频域的分析，了解线性时不变系统的特性，同时加深对信号频谱的理解。
2、课程设计原理
1离散时间信号卷积使用MATLAB中的conv函数
2.在MATLAB中，filter函数的作用是进行数字滤波，其实也就是系统的响应序列求解；
y=filter（b,a,x）;
y=filter（b,a,x）表示向量b和a组成的系统对输入x进行滤波，系统的输出为y。
3.在MATLAB中，zplane函数——绘制系统零极点图
4.在MATLAB中，有专门的frepz（）函数来求解。
调用形式为
[h
w]=freqz（b,a,,n,fs）其中向量b和a为离散系统的系数，fs为采样频率，n为在区间[0
fs/2]频率范围内选取的频率点数，f记录频率点数。
由于freqz函数是采用基2的FFT的算法，n常取2的幂次方，以便提高计算速度。
5.MATLAB提供了如abs（）,angle（）等基本函数用来计算DTFT幅度，相位。
3、课程设计所用设备
（1）个人电脑一台
（2）MATLAB2012设计程序
4、课程设计内容和步骤
（1）根据设计题目的要求，熟悉相关内容的理论基础，理清程序设计的措施和步骤；
（2）根据设计题目的要求，提出各目标的实施思路、方法和步骤；
（3）根据相关步骤完成MATLAB程序设计，所编程序应能完整实现设计题目的要求；
（4）调试程序，分析相关理论；
（5）编写设计报告。
5、设计过程、实现过程、实现结果
（一）基本部分
（1）信号的时频分析
任意给定单频周期信号的振幅、频率和初相，要求准确计算出其幅度谱，并准确画出时域和频域波形，正确显示时间和频率。
设计思路：
首先给出横坐标，即时间，根据设定的信号的振幅、频率和初相，写出时域波形的表达式；
然后对时域波形信号进行傅里叶变化，得到频域波形；
最后使用plot函数绘制各个响应图。
源程序：
%基本题1,时频分析，2011010512
clc,clear,closeall
A=4;
fc=200;
a=0*pi/3;
T=1/fc;
fs=15*fc;
Ts=1/fs;
t=0:
Ts:
3*T-Ts;
N=length（t）;
x=A*sin（2*pi*fc*t+a）;
subplot（211）,plot（t,x）
xlabel（'
t（s）'
）,ylabel（'
x（t）'
）,title（'
时域波形图'
）
X=fft（x）;
XX=fftshift（X）;
magX=abs（XX）;
phaseX=angle（XX）;
f=（0:
N-1）*（fs/N）-fs/2;
subplot（223）,plot（f,magX）
f（Hz）'
|X（jw）|'
信号幅度谱'
subplot（224）,plot（f,phaseX）
angle[X（jw）]'
信号相位谱'
结果图：
结果分析：
cos函数波形为周期信号，其频域响应为两个冲击函数，并且符合对偶性。
（2）傅里叶级数分析
分析周期锯齿波的傅里叶级数系数，用正弦信号的线性组合构成锯齿波，要求谐波次数可以任意输入，分析不同谐波次数所构成的锯齿波，解释是否存在吉伯斯现象。
根据原理——任意周期信号都可以表达成傅里叶级数的形式，对周期举行锯齿波进行傅里叶级数分解，利用for循环完成级数求和运算。
当N值较小时，傅里叶级数的逼近效果不是很理想，随着N值变大，傅里叶级数越来越接近理想值。
%基本题2，傅里叶级数分析，2011010512
clc;
clear;
closeall;
%初始化工作空间
t=-6:
6/1000:
6;
%时间
N=input（'
pleaseinputthenunber:
'
）;
w0=pi;
%角频率
XN=0;
forn=1:
N;
fn=（-1）^（n+1）/（n*pi）;
XN=XN+fn*sin（n*w0*t）;
end
XN=XN+0.5;
%抬升0.5
figure;
plot（t,XN）;
title（['
Gibbs,N='
num2str（N）]）;
Time（sec）'
ylabel（['
X'
num2str（N）,（'
t'
）]）;
结果图；
随着N值的增大，傅里叶级数与理想的周期锯齿波越来越接近，符合预期效果。
（3）系统分析
任意给定微分方程或差分方程描述的系统，画出系统的幅频响应和相频响应。
根据微分方程、差分方程与系统函数的对应关系，结合matlab自带的freqs和freqz两个函数，分析系统的幅频特性和相频特性。
%基本题3系统分析2011010512
clc,clear,closeall,
%微分方程为r'
（t）+6r'
（t）+8r（t）=2*e（t）
b=[0,0,2];
%系统函数H（s）分子多项式系数
a=[1,6,8];
%系统函数H（s）分母多项式系数
[Hw1,w1]=freqs（b,a）;
%求连续时间系统频率响应
figure
（1）
subplot（2,1,1）;
zplane（b,a）;
title（'
连续时间系统零极点图'
subplot（2,2,3）;
plot（w1/pi,abs（Hw1））
rad/s'
ylabel（'
|H（jw）|'
连续时间系统幅频响应'
subplot（2,2,4）;
plot（w1/pi,angle（Hw1））
Phase'
连续时间系统相频响应'
%差分方程为y（n）-0.5*y（n-1）+4/5*y（n-2）=3*x（n）-0.6*x（n-1）;
b=[3,-0.6];
%系统函数H（z）分子多项式系数
a=[1,-0.5,4/5];
%系统函数H（z）分母多项式系数
[Hw,w]=freqz（b,a）;
%求离散时间系统频率响应
figure
（2）
离散时间系统零极点图'
plot（w/pi,abs（Hw））
w/pi）'
|H（e^j^w）|'
离散时间系统幅频响应'
plot（w/pi,angle（Hw））
w/pi'
离散时间系统相频响应'
所显示零极点等图符合验算结果
（4）音乐合成程序设计
对于任意一小段音乐，利用“十二平均律”计算该音乐中各个乐音的频率，产生并播放这些乐音。
分析音乐的频谱，从中识别出不同的乐音。
根据时间长短来区别各个音符拍长短，根据频率高低来区别各个音符的音调，具体细节参考“十二平均律”，将每段乐音连接起来，使用sound函数播放乐音。
对每段乐音进行傅里叶变换，分析其频谱。
%基本题4，音乐合成程序设计，一分钱，2011010512
fs=8e3;
freq=[384,512,426.66,512,426.66,320,384,288,320,384,320,384,426.66,512,384,426.66,384,320,384,256,288,320,288,256,288,320,436.66,384,320,384,436.66384,512,426.66,384,320,384,288,384,288,320,256];
N=[1/2,1/2,1/4,1/4,1/2,1/4,1/4,1/4,1/4,1,1/4,1/4,1/4,1/4,1/4,1/4,1/4,1/4,1/4,1/2,1,1/4,1/4,1/4,1/4,1,1/4,1/4,1/4,1/4,1,1/4,1/4,1/4,1/4,1/4,1/4,1/2,1/4,1/4,1/2,1/2]*fs;
D=zeros（1,sum（N））;
X=1;
fork=1:
42
D（X:
X+N（k）-1）=（sin（2*pi*freq（k）*（1/fs:
1/fs:
N（k）/fs）））;
X=X+N（k）;
end;
sound（D,fs）;
结果分心：
所听音乐符合预期。
（5）调制
分析单位冲激响应为
的系统的滤波特性，画出其幅频响应曲线。
根据傅里叶变换特性，时域相乘对应频域卷积，将题目中的时域波形进行分解，并将两个函数分别对应的傅里叶变换进行卷积。
参考常用的傅里叶变换对，所求的频域特性为方波和冲击函数的卷积。
%基本题5，调制，2011010512
fs=300;
t=-10:
10;
n=[-（N-1）/2:
（N-1）/2];
h=sin（50*pi*t）.*sin（100*pi*t）./（pi*t+0.00001）;
f1=fft（h）;
Fn=fftshift（f1）;
f=n/N*fs;
plot（t,h）,
Amplitude'
H（t）'
plot（f,abs（Fn））,
w'
Modulation'
所得结果符合演算过程，与预期一致。
（二）提高部分
题目1、音乐合成2011010512
设计思路：
通过定义每个音来写乐谱，在加包络谐波，在播放音乐，后画出波形。
源程序
%提高题音乐合成
clear;
closeall;
TUNE_BASED=220*[2^（3/12）2^（5/12）2^（7/12）2^（8/12）2^（10/12）2^（12/12）2^（14/12）];
%C调
H_TUNE=2*TUNE_BASED;
L_TUNE=0.5*TUNE_BA
升级会员 