信号与系统课程设计报告材料
信号与系统课程设计-傅里叶变换及matlab仿真
则 （1-11）
（10）时域卷积定理：
则 （1-12）
傅立叶变换及逆变换的MATLAB实现
MATLAB的Symbolic Math Toolbox提供了能直接求解傅里叶变换及逆变换的函数fourier()及ifourier()。
三、设计的方法及步骤：
(1) F=fourier(f)
(2) F=fourier(f,v)
(3) F=fourier(f,u,v)
说明：(1) F=fourier(f)是符号函数f的傅立叶变换，缺省返回是关于ω的函数。如果f=f(ω),则fourier函数返回关于t的函数。
(2)F=fourier(f,v)返回函数F是关于符号对象v的函数，而不是默认的ω，即
傅立叶逆变换定义是： （1-2）
称为 的频谱密度函数。
傅立叶变换的性质
（1）线性性质：
（1-3）
（2）频移性质： （1-4）
（3）时移性质： （1-5）
（4）尺度变换性质： （1-6）
（5）对称性质： （1-7）
（6）时域微分性质： （1-8）
（7）频域微分性质： （1-9）
（8）时域积分性质： （1-10）
axis([-1,3,-0.2,1.2])
j=sqrt(-1);
F=1./(j*t);
y=pi*imp(t);
subplot(1,2,2)
plot(t,abs(F));
axis([-1,1,0,20]);
ylabel('F(jw)');
xlabel('w');
d on,
plot(t,y);
5、冲激函数：
傅里叶分析方法不仅应用于电力工程、通信和控制领域之中，而且在力学、光学、量子物理和各种线性系统分析等许多有关数学、物理和工程技术领域中得到广泛而普遍的应用。
信号与系统课程设计(信号调制与解调)(采样定理)(LTI系统分析)
课题一信号调制与解调题目说明：从语音，图像的原始信息变过来的原始信号频谱分量频率较低，不适宜在信道中长距离传输。
因此，在通信系统的发送通端常需要有调制过程将其转换为适合传输的信号，在接收端则需要有调节过程，将信号还原成原来的信息，以便更准确的利用信息。
原理分析：调制就是按调制信号的变化规律去改变某些参数。
解调是调制的逆过程，即从已调制信号中恢复或提取调制信号的过程。
幅度调制是正弦型载波的幅度随调制信号变化的过程。
采用模拟调制利用正旋波载波的幅度调制，频率调制和相位调制的方式进行信号的处理。
同步解调端本振信号频率必须与发射端调制的载波信号的频率和相位相同才能实现同步解调。
脉冲调制信号只有在脉冲出现才需要存在，在其他时间内等于零，这样就有可能在这空余的时间间隔中去传输其他路德信号，发送端和接受端的转换开关按照同样的顺序和周期轮流接通各个通道，在信道中传送的是各个脉冲幅度调制信号的和，各个脉冲出现在不同的时间段。
而通过接收端的开关以后各路接受端接收到的相当于某一路信号脉冲幅度的结果，可以用低通滤波器进行解调。
实验内容：1.将一正旋信号x(n)=sin(2πn/256)分别以100000Hz的载波和1000000Hz的取样频率进行调制,写出MATLAB脚本实现抑制载波幅度调制，实现同步解调，滤波输出的波形。
2.分别作出cos(10t)cos(w c t)和[1+0.5sin(10t)]cos(w c t)的波形图和频谱图，并对上面调制信号进行解调，观察与源图的区别。
模块设计1：1.产生一个输入信号 2.产生一个载波信号3.构造用于解调的低通滤波器4.低通滤波解调5.画图MATLAB程序1：>> clear; %清除已存在变量n=0:0.0001:256; %自变量e=sin(2*pi*n/256); %调治信号s=cos(100000*n); % 载波信号a=e.*s; % 调制b=a.*s; % 解调[nb,na]=butter(4,100,'s'); % 低通滤波sys=tf(nb,na); % 构建sys对象c=lsim(sys,b,n); %低通滤波subplot(2,2,1) % 图形输出语句plot(n,e);title('调制信号'); %图形标题>> xlabel('n'),ylabel('e(n)'); %横纵坐标变量>> grid on %坐标网格>> subplot(2,2,2) % 图形输出语句>> plot(n,a);>> title('调幅信号'); %图形标题>> xlabel('n'),ylabel('a(n)'); %横纵坐标变量>> grid on %坐标网格>> subplot(2,2,3) % 图形输出语句>> plot(n,b);>>title('解调波形'); %图形标题>> xlabel('n'),ylabel('b(n)'); %横纵坐标变量>> grid on %坐标网格>> subplot(2,2,4) % 图形输出语句>> plot(n,c);>> title('滤波后的波形')；%图形标题>>xlabel('n'),ylabel('e(n)'); %横纵坐标变量>> grid on %坐标网格模块设计2：1.产生两个输入信号 2.用克诺内科内积产生两个周期行序列脉冲3.调制并向加4.构造用于解调的低通滤波器5.低通滤波解调 6画图MATLAB程序2：>> clear; % 清除变量t=0:0.001:9.999; % 定义自变量取值范围和间隔e1=cos(10*t).*cos(600*t); % 输入信号e2=(1+0.5*sin(10*t)).*cos(600*t); %输入信号p0=ones(1,2500);p1=kron(p0,[1,0,0,0]); %第一个序列脉冲p2=kron(p0,[0,0,1,0]); % 第二个序列脉冲a=p1.*e1+p2.*e2; 调制并向加[nb,na]=butter(4,20,'s'); % 用于解调的低通滤波器sys=tf(nb,na); %构建sys对象b1=a.*p1; % 取得第一路信号的脉冲调制信号c1=lsim(sys,b1,t);%通过低通滤波解调输出b2=a.*p2; %取得第二路信号的脉冲调制信号c2=lsim(sys,b2,t); % 通过低通滤波解调输出subplot(4,2,1) % 图形输出语句plot(t,e1);title('第一路输出信号'),xlabel('t'),ylabel('e(t)');grid on%图形横纵坐标，标题，坐标网格subplot(4,2,2) % 图形输出语句plot(t,e2);title('第二路输出信号'),xlabel('t'),ylabel('e(t)');grid on%图形横纵坐标，标题，坐标网格subplot(4,2,3) % 图形输出语句plot(t,e1.*p1);title('第一路脉冲调制信号'),xlabel('t'),ylabel('e(t)');grid on %图形横纵坐标，标题，坐标网格subplot(4,2,4) % 图形输出语句plot(t,e2.*p2);title('第二路脉冲调制信号'),xlabel('t'),ylabel('e(t)');grid on %图形横纵坐标，标题，坐标网格subplot(4,2,5) % 图形输出语句plot(t,a);title('合成的传输信号'),xlabel('t'),ylabel('e(t)');grid on%图形横纵坐标，标题，坐标网格subplot(4,2,6) % 图形输出语句plot(t(5001:5250),a(5001:5250));title('局部放大后的合成信号'),xlabel('t'),ylabel('e(t)');grid on%图形横纵坐标，标题，坐标网格实验总结：通过对理论知识的学习，使自己对信号的调制与解调具有一定的认知水平，然后开始做实验，此时要理论结合实践，作出波形图后要考虑与理论波形进行比较，比较的方法是，首先判断所测波形是否正确，若不正确找出错误原因，若正确则分析实测波形与理论波形不完全相同的原因。
信号与系统课程设计报告
二○一一～二○一二学年第一学期电子信息工程系信号与系统课程设计报告班级：电子信息工程2009级3班学号：200904135104姓名：徐奎课程名称：数字信号处理课程设计学时学分：1周1学分指导教师：陈华丽二○一一年十二月三十日1、课程设计目的：数字信号处理”课程是信息和通信工程专业必修的专业技术基础课程，课程以信号与系统作为研究对象，研究对信号进行各种处理和利用的技术。
通过该课程的学习，学生应牢固掌握确定性信号和系统的分析方法、相关算法、系统实现等的相关知识的，借助于数字滤波器的设计及实现，学生可掌握数字系统的分析以及设计方法。
数字信号处理是理论性和工程性都很强的学科，本课程设计的目的就是使该课程的理论与工程应用的紧密结合, 使学生深入理解信号处理的内涵和实质。
本课程设计要求学生在理解信号处理的数学原理的基础上，应用计算机编程手段，实现一种信号分析或处理的设计，达到对所学内容融会贯通，综合各部分知识，按照题目要求独立设计完成。
2、课程设计内容：滤波器设计产生一个连续信号，包含低频，中频，高频分量，对其进行采样，进行频谱分析，分别设计低通，带通，高通滤波器对信号进行滤波处理，观察滤波前后信号的频谱。
独立完成以上设计，有能力的同学设计一个友好的人机交互界面，不限编程语言。
3、设计内容和步骤：①设定的连续信号为：s=sin(2*pi*t*5)+sin(2*pi*t*15)+sin(2*pi*t*30)可知：信号中包含了5Hz、15Hz、30Hz频率分量，对其采样的频率取100Hz。
用plot函数画出其时域波形，代码如下，结果如下图所示：% 程序功能：产生一个连续信号，包含低频，中频，高频分量，对其进行采样，进行频谱分析，分别设计低通，带通，高通滤波器对信号进行滤波处理，观察滤波前后信号的频谱。
Fs=100;t=(1:100)/Fs;s1=sin(2*pi*t*5);s2=sin(2*pi*t*15);s3=sin(2*pi*t*30);s=s1+s2+s3;figure(1);plot(t,s); % 画出信号的时域波形xlabel('t');ylabel('s');title('原始信号的时域波形');% 程序功能：画出信号的频谱图。
《信号与系统》课程思政教学设计
《信号与系统》课程思政教学设计一、教学目标1. 知识与技能掌握信号与系统的基础理论和分析方法。
能够应用所学知识解决实际工程问题。
2. 思政目标培养学生的爱国情怀和科学精神。
增强学生的职业道德和社会责任感。
提升学生的创新思维和团队协作能力。
二、教学内容与方法1. 教学内容信号与系统的基本概念、分类及性质。
信号的时域和频域分析。
系统的稳定性、因果性和线性时不变性。
2. 思政元素融入引入我国科学家在信号与系统领域的研究成果，激发学生的民族自豪感和科学探索精神。
讨论信号与系统在国家安全、通信、医疗等领域的应用，培养学生的社会责任感和职业道德。
3. 教学方法理论讲授：系统介绍信号与系统的基本理论和方法。
案例分析：结合实际应用案例，分析信号与系统的实际应用。
小组讨论：组织学生围绕思政主题进行小组讨论，促进思想交流和团队协作。
课程设计：安排与课程内容相关的设计任务，提升学生的实践能力和创新思维。
三、思政教学重点1. 科学精神培养通过介绍信号与系统领域的发展历程和科学家事迹，培养学生的科学探索精神和创新意识。
鼓励学生勇于挑战传统观念，追求科学真理。
2. 职业道德教育强调工程师的职业道德和社会责任，引导学生在未来职业生涯中坚守诚信、公正和负责任的原则。
通过案例分析，讨论工程实践中的道德困境和解决方案。
3. 团队协作与沟通能力提升通过小组讨论和课程设计等环节，锻炼学生的团队协作和沟通能力。
培养学生学会倾听他人意见、尊重他人观点并有效表达自己的思想。
四、教学评价与反馈机制1. 知识掌握评价通过作业、测验和考试等方式评价学生对信号与系统知识的掌握情况。
2. 思政表现评价观察并记录学生在课堂讨论、小组活动和课程设计中的思政表现。
将思政表现纳入课程考核体系，激励学生积极参与思政教育活动。
3. 教学反馈定期收集学生对课程内容和教学方法的反馈意见，及时调整教学策略以满足学生需求。
与学生保持良好沟通，及时解答学生在学习和思政方面的困惑和问题。
《信号与系统》课程设计
《信号与系统》课程设计
（3）
• • • • • • • • • • • • • • • • • • •
t=0:0.01:40; w1=exp(-3*t); subplot(3,2,1);plot(t,w1);axis([0,3,-0.2,2]);grid; title('f1复数模随时间变化的绘图') w2=exp(-3*t); subplot(3,2,2);plot(t,-w2);axis([0,3,-2,0.2]);grid; title('-f1（t）'); w3=exp((-3)*(-t))axis([0,100,-0.2,50]);grid; title('f1（-t）'); w4=exp(-3*2*t); subplot(3,2,4);plot(t,w4);axis([0,2,-0.2,2]);grid; title('f1（2t）'); w5=exp(-3*(t+2)); subplot(3,2,5);plot(t,w5);axis([0,3,-0.2,5]);grid; title('f1（t+2）'); w6=exp(-3*(2-2*t)); subplot(3,2,6);plot(t,w6);axis([0,2,-0.2,5]);grid; title('f1（2-2t）');
一.设计目的
1.加深对信号与系统的课本知识的理解和应用。
2.加深和巩固对典型信号：复指数信号的学习和理解，分析实部、虚 部、模及相角随时间变化的曲线并了解其时域特性。 3.应用MATLAB对实际问题进行仿真，通过对课程实践的制作，加深 对信号的时移、翻转、放缩的理解和掌握。
《信号与系统》课程设计
信号与系统课程设计
沈阳大学沈阳大学3．3系统与连续时间信号系统是连续事物或各个部分的一个复杂的整体，有形或无形事物的组成体。
系统可以分为即时系统与动态系统；连续系统与离散系统；线性系统与非线形系统；样时变系统和非时变系统等等。
在连续时间系统中，如一个连续时间系统接收，根据定义在连续时间（-∞<t<∞）有定义的信号称为连续时间信号，在范围内输入信号x(t)，并产生输出信号y(t)。
连续时间信号是在连续时间范围内定义的信号值，信号的幅值可以是连续数值，也可以是离散数值。
当信号幅值连续是，则称之为模拟信号。
3．4采样定理取样定理论述了在一定条件下，一个连续时间信号完全可以用该信号在等时间间隔上的瞬时值（或称样本值）表示，这些样本值包含了连续时间信号的全部信息，利用这些样本值可以恢复原信号。
可以说取样定理在连续时间信号与离散时间信号中架起了一座桥梁。
其具体内容如下：取样定理：设为带限信号，带宽为0F ，则当取样频率02F F s ≥时，可从取样序列)()(s a nT x n x =中重构，否则将导致)(n x 的混叠现象。
带限信号的最低取样频率称为Nyquist （奈奎斯特）速率。
图1给出信号采样原理图图1 信号采样原理图由图1可见，)()()(t t f t f Ts s δ⋅=，其中，冲激采样信号)(t Ts δ的表达式为：∑∞-∞=-=n sT nT t t s)()(δδ （1）其傅立叶变换为∑∞-∞=-n s s n )(ωωδω，其中ss T πω2=。
设)(ωj F ，)(ωj F s 分别为)(t f ，)(t f s 的傅立叶变换，由傅立叶变换的频域卷积定理，可得：沈阳 大 学∑∑∞-∞=∞-∞=-=-=n ssn s s s n j F T n j F j F )]([1)(*)(21)(ωωωωδωωπω （2）若设)(t f 是带限信号，带宽为m ω如图（2），由式（2）可见，)(t f 经过采样后的频谱)(ωj F s 就是将)(ωj F 在频率轴上搬移至 ,,,,,02ns s s ωωω±±±处（幅度为原频谱的s T 1倍）。
信号与与系统课程设计
信号与与系统课程设计一、教学目标本节课的教学目标是让学生掌握信号与系统的基本概念、原理和分析方法。
具体包括：1.知识目标：–了解信号与系统的定义、特点和分类；–掌握信号的时域、频域分析方法；–理解系统的基本特性，如线性、时不变性等。
2.技能目标：–能够运用信号与系统的分析方法解决实际问题；–熟练使用相关软件工具进行信号处理和系统分析；–具备一定的科研能力和创新精神。
3.情感态度价值观目标：–培养对信号与系统学科的兴趣和热情；–树立正确的科学观，注重实践与理论相结合；–增强团队协作意识，提高沟通与表达能力。
二、教学内容本节课的教学内容主要包括以下几个部分：1.信号与系统的定义、特点和分类；2.信号的时域、频域分析方法；3.系统的基本特性，如线性、时不变性等；4.实际应用案例分析。
5.引言：介绍信号与系统课程的背景、意义和目标；6.信号与系统的定义、特点和分类：讲解信号与系统的概念，分析各种信号与系统的特点和分类；7.信号的时域、频域分析方法：讲解信号的时域、频域分析方法，并通过实例进行分析；8.系统的基本特性：讲解系统的基本特性，如线性、时不变性等，并通过实例进行分析；9.实际应用案例分析：分析信号与系统在实际应用中的案例，如通信系统、控制系统等。
三、教学方法为了提高教学效果，本节课将采用以下教学方法：1.讲授法：讲解信号与系统的基本概念、原理和分析方法；2.讨论法：学生进行课堂讨论，培养学生的思考能力和团队协作精神；3.案例分析法：分析实际应用案例，让学生更好地理解信号与系统的应用价值；4.实验法：安排课后实验，让学生动手实践，提高实际操作能力。
四、教学资源为了支持教学内容和教学方法的实施，本节课将准备以下教学资源：1.教材：选用权威、实用的教材，如《信号与系统》、《信号处理与系统分析》等；2.参考书：提供相关领域的参考书籍，如《线性系统理论》、《数字信号处理》等；3.多媒体资料：制作精美的PPT课件，提供动画、视频等多媒体资料；4.实验设备：准备相应的实验设备，如信号发生器、示波器、滤波器等，以便进行课后实验。
高校青教赛 信号与系统教学设计范例
高校青教赛信号与系统教学设计范例信号与系统课程设计教案一、matlab工作空间介绍。
二、信号处理部分：1）信号的产生，matlab工具箱，自己编程函数仿真，导入实际数据。
2）信号的卷积，奇偶分解，各种性质的验证。
3）信号分解的基本原理。
4）信号分解的算法实现，自己编程验证。
5）结合实验给出实验分析和结论。
三、离散信号处理部分：1）信号分解算法的离散化。
2）信号分解的基本原理。
3）信号分解的算法实现，自己编程验证。
4）结合实验给出实验分析和结论。
四、信号滤波处理部分：1）将信号进行傅里叶分解。
2）在频率域进行理想滤波。
3）将信号变换到时间域。
4）结合实验结果给出实验分析和结论。
五、连续系统分析部分：1）电路系统建模或者已有微分系统方程。
2）根据输入求解系统的响应。
3）求解系统的单位冲激响应。
4）编程实现，验证系统的因果性，稳定性。
六、离散系统分析部分：1）电路系统建模或者已有差分系统方程。
2）根据输入求解系统的响应。
3）求解系统的单位脉冲响应。
4）编程实现，验证系统的因果性，稳定性。
实验报告组成：1、实验基本原理2、理论分析求解3、实验编程验证4、实验结果分析。
一、基本函数：1、函数变量的定义。
syms是定义符号变量sym是将字符或者数字转换为字符比如syms x y %就是定了符号变量x y以后x y就可以直接使用sys('a+b')%就是将a+b转化为符号表达式。
2、单位阶跃信号。
Heaviside()。
syms t;f=heaviside(t-4)；或者f=@(t)heaviside(t-4); ezplot(f,[0 5])3、单位冲激信号f=@(x)dirac(x-2);二、示例演示分析示例1：1设f(t) e 2tu(t)，画出该信号的及其幅频图。
21、概述：掌握信号傅立叶变换的计算方法。
2、设计任务，即要设计的主要内容和要求等掌握信号傅立叶变换的计算方法以及程序求解方法。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
课程设计报告课程名称信号与系统课程设计指导教师设计起止日期学院信息与通信工程专业电子信息工程学生班级/学号成绩指导老师签字目录1、课程设计目的 (1)2、课程设计要求 (1)3、课程设计任务 (1)4、课程设计容 (1)5、总结 (11)参考文献 (12)附录 (12)1、课程设计目的“信号与系统”是一门重要的专业基础课，MATLAB作为信号处理强有力的计算和分析工具是电子信息工程技术人员常用的重要工具之一。
本课程设计基于MATLAB完成信号与系统综合设计实验，以提高学生的综合应用知识能力为目标，是“信号与系统”课程在实践教学环节上的必要补充。
通过课设综合设计实验，激发学生理论课程学习兴趣，提高分析问题和解决问题的能力。
2、课程设计要求（1）运用MATLAB编程得到简单信号、简单信号运算、复杂信号的频域响应图；（2）通过对线性时不变系统的输入、输出信号的时域和频域的分析，了解线性时不变系统的特性，同时加深对信号频谱的理解。
3、课程设计任务（1）根据设计题目的要求，熟悉相关容的理论基础，理清程序设计的措施和步骤；（2）根据设计题目的要求，提出各目标的实施思路、方法和步骤；（3）根据相关步骤完成MATLAB程序设计，所编程序应能完整实现设计题目的要求；（4）调试程序，分析相关理论；（5）编写设计报告。
4、课程设计容（一）基本部分（1）信号的时频分析任意给定单频周期信号的振幅、频率和初相，要求准确计算出其幅度谱，并准确画出时域和频域波形，正确显示时间和频率。
设计思路：首先给出横坐标，即时间，根据设定的信号的振幅、频率和初相，写出时域波形的表达式；然后对时域波形信号进行傅里叶变化，得到频域波形；最后使用plot函数绘制各个响应图。
源程序：clc; clear;close all;Fs =128; % 采样频率T = 1/Fs; % 采样周期N = 600; % 采样点数t = (0:N-1)*T; % 时间，单位：Sx=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);endXN=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或者16000t2=0:1/f:1; % 2拍，时间长短不同t4=0:1/f:0.5; % 1拍t8=0:1/f:0.25; % 1/2拍misc_note5=523.25; % 不同音符频率misc_note6=587.33;misc_note2=392;misc_note1=349.23;misc_note6_down=293.66;m1=sin(2*pi*misc_note5*t4); % 波形m2=sin(2*pi*misc_note5*t8);m3=sin(2*pi*misc_note6*t8);m4=sin(2*pi*misc_note2*t2);m5=sin(2*pi*misc_note1*t4);m6=sin(2*pi*misc_note1*t8);m7=sin(2*pi*misc_note6_down*t8);m8=sin(2*pi*misc_note2*t2);m=[m1 m2 m3 m4 m5 m6 m7 m8];% save('DongFangHong','m');sound(m);N = 1024; % 傅里叶变换点数flabel = 0:f/N:f/2-f/N;figure; hold on;H_m = fft(m1,N); plot(flabel,10*log10(abs(H_m(1:N/2))));H_m = fft(m2,N); plot(flabel,10*log10(abs(H_m(1:N/2))));H_m = fft(m3,N); plot(flabel,10*log10(abs(H_m(1:N/2))));H_m = fft(m4,N); plot(flabel,10*log10(abs(H_m(1:N/2))));H_m = fft(m5,N); plot(flabel,10*log10(abs(H_m(1:N/2))));H_m = fft(m6,N); plot(flabel,10*log10(abs(H_m(1:N/2))));H_m = fft(m7,N); plot(flabel,10*log10(abs(H_m(1:N/2))));H_m = fft(m8,N); plot(flabel,10*log10(abs(H_m(1:N/2))));hold off;xlabel('频率');ylabel('幅度响应');grid on;axis([0 1000 -15 30]);结果图；结果分析：根据乐谱编写程序，能够听到预期的歌曲；分析其频谱，包含几个不同频率的脉冲，同样符合预期。
（5）调制分析单位冲激响应为πππ=sin(50)sin(100)h(t)t tt的系统的滤波特性，画出其幅频响应曲线。
设计思路：根据傅里叶变换特性，时域相乘对应频域卷积，将题目中的时域波形进行分解，并将两个函数分别对应的傅里叶变换进行卷积。
参考常用的傅里叶变换对，所求的频域特性为方波和冲击函数的卷积。
源程序：clc; clear;close all;f = 800; % 采样频率t = -0.2:1/f:0.2; % 时域长度x_t = zeros(1,length(t));for k=1:length(t)if t(k)==0x_t(k) = 0;elsex_t(k) = sin(50*pi*t(k))*sin(100*pi*t(k))/(pi*t(k));endendfigure; subplot(2,1,1);plot(t,x_t); % 时域波形xlabel('时间');ylabel('时域波形');grid on;N = 1024; % 傅里叶变换点数H_x = fft(x_t,N);flabel = 0:f/N:f/2-f/N;subplot(2,1,2);plot(flabel,10*log10(abs(H_x(1:N/2))));xlabel('频率');ylabel('频率响应');grid on;结果图；结果分析：中心在0频的方波搬移到了50Hz，符合预期。
（二）提高部分（7）频分复用a、自行给出二路语音信号，分别显示其频谱，并播放语音。
b、对二路语音信号进行频分复用，显示复用后的频谱，播放语音。
c、设计程序对频分复用的信号进行解调，显示解调结果，并回放语音。
设计思路：对信号的处理主要包括以下几个部分：a、根据语音信号文件，使用matlab的函数wavread进行读取，并结合fft观察信号的时域波形和频域波形；b、两路信号分别于不同频率的载波相乘，再相加，得到复用之后的信号；c、根据所发送信号的特性和其对应的载波，设计带通滤波器，从复用信号中分别滤出两路信号；d、对高频的调制信号进行解调至0频；e、设计低通滤波器，并将解调后的信号通过该滤波器，得到预期的发送信号，播放该信号，观察其与发送信号的异同。