信号与系统 2023(春季) 作业要求 - 第七次作业
01 基础练习
一、傅里叶变换性质
1、傅里叶变换
请使用傅里叶变换的性质求解下面信号的频谱。如果信号中包括有 f ( t ) f\left( t \right) f(t) , 则它对应的傅里叶变换为 F ( ω ) F\left( \omega \right) F(ω) 。
(1)必做题
(1) 双三角形信号
请分别应用傅里叶变换的微分性质、卷积性质 两种方法求取下面梯形信号的傅里叶变换。
▲ 图1.1.1 双三角形脉冲信号
微分和卷积定理求解信号频谱
(2)
提示: 需要应用到 频域微分特性。
(3)
提示:需要应用到时域微分和频域积分特性。
(4) 调制高斯信号
请使用傅里叶变换频移特性 求解调制高斯信号的频谱。
▲ 图1.1.2 调制高斯信号
提示: 关于高斯信号的傅里叶变换可以参照下面给出的公式。
▲ 图1.1.3 高斯信号的傅里叶变换
(5)
(2)选做题
(6)
请使用傅里叶变换的微分特性求取下面锯齿脉冲信号的频谱。
▲ 图1.1.4 锯齿脉冲信号
(7)
请使用傅里叶变换的频移特性求下面单个正弦脉冲信号的频谱。
▲ 图1.1.5 单个正弦脉冲信号
(8)
(9)
应用FT性质求解频谱
2、傅里叶反变换
(1)必做题
(1) 已知信号的频谱如下, 求 f ( t ) f\left( t \right) f(t) 。
提示: 利用频域微分定理, 并结合
(2)选做题
(2)
提示:参见前面第一题的 第 (9) 小题的答案。
(3)
求解下面“信号平方频谱” 对于的时域信号 f ( t ) f\left( t \right) f(t) 的表达式。
\
求解傅里叶反变换
3、信号平方的频谱
已知信号
f
(
t
)
f\left( t \right)
f(t) 的频谱为:
如下图所示:
▲ 图1.1.6 信号的频谱
请绘制出 f 2 ( t ) f^2 \left( t \right) f2(t) 的 频谱图。
提示: 这个题目需要应用到 FT 的频域卷积定理, f 2 ( t ) f^2 \left( t \right) f2(t) 的频谱实际上是它自己的频谱与自己卷积, 再考虑前面还有 2 π 2\pi 2π 的因子。
信号平方的频谱
二、周期信号的频谱
已知单个梯形脉冲和单个余弦脉冲的傅里叶变换,求下图所示的周期梯形信号和周期全波信号的傅里叶变化。
1、必做题
(1) 全波整流周期信号
单个半波余弦脉冲信号的傅里叶变换为:
▲ 图1.2.1 全波整流信号波形
提示: 将单个周期信号的傅里叶变换进行“离散化”。
(2) 周期信号与三角波信号乘积
求下面信号的频谱。 该信号是周期为 τ \tau τ , 幅度为 ±1 的三角波信号, 与宽度为 τ 1 \tau _1 τ1 , 高度为 1 的对称三角脉冲信号相乘。
▲ 图1.2.2 调制在三角波上的三角波信号
提示: 关于等腰三角形脉冲信号的频谱可以查询课件中给出的答案。 使用傅里叶变换时域乘积对应频域卷积定理。
2、选做题
单个梯形信号的傅里叶变换:
▲ 图1.2.3 周期梯形信号
提示: 将单个周期信号的傅里叶变换进行“离散化”。
周期信号傅里叶变换
三、证明题
1、必做题
已知信号
f
(
t
)
f\left( t \right)
f(t) 的“希尔伯特变换” 的表达式为:
定义信号
f
(
t
)
f\left( t \right)
f(t) 的解析信号 为:
试证明:
其中:
F
(
ω
)
,
Z
(
ω
)
F\left( \omega \right),Z\left( \omega \right)
F(ω),Z(ω) 分别是
f
(
t
)
,
z
(
t
)
f\left( t \right),z\left( t \right)
f(t),z(t) 的傅里叶变换。
u
(
ω
)
u\left( \omega \right)
u(ω) 是关于
ω
\omega
ω 的单位阶跃函数。
解析信号的频谱
02 实验作业
一、观察信号频谱变化
1、作业要求
在网络学堂下载一段音乐及其经过处理后的音乐数据文件,聆听相应的音频效果。
使用MATLAB中使用spectrogram 命令观察音乐及其变换后的数据是时频联合分布, 总结节奏变化一是尺度变化两种操作在听觉和信号时频联合分布之间的差别。
2、作业相关文件说明
音乐数据文件与作业要求文档合并在一起,HMW7-DATA.ZIP文件。
HMW7-DATA中的文件说明:_
1.mp3 :原始音乐,MP3格式
1.wav :原始音乐,WAV格式
1fast.wav :加快节奏50%文件
1fast2.wav :加快节奏100%文件
1fastrate.wav:尺度减少50%文件;
1fastrate2.wav:尺度减小75%文件
1slow.wav :减慢节奏50%文件;
1slow2.wav :减慢节奏100%文件
1slowrate.wav:尺度拉伸50%文件
1slowrate2.wav:尺度拉伸100%文件
soundstretch.exe:处理程序
3、MATLAB相关命令
>> [mdata, fs] = audioread('d:\'temp\'1.wav');'
>> spectrogram(mdata(1:fs*2,1), 2048, 1024, 2048, fs, 'yaxis')
附录:
本次作业中部分图形绘制程序:
#!/usr/local/bin/python
# -*- coding: gbk -*-
#============================================================
# TEST1.PY -- by Dr. ZhuoQing 2023-04-06
#
# Note:
#============================================================
from headm import *
t = linspace(-5, 5, 10000)
def G(t, startn, endn):
return heaviside(t-startn,0.5)-heaviside(t-endn,0.5)
def Gt(t, center, width):
startn = center-width/2
endn = startn + width
return heaviside(t-startn,0.5)-heaviside(t-endn,0.5)
#ft = (t+2)*G(t,-2,-1) - t*G(t,-1,0) + t*G(t,0,1) - (t-2)*G(t,1,2)
#ft = (-t) * G(t, -2, 2) / 3
#ft = -sin(t*pi/2)*G(t,-2,2)
#ft = G(t,-4,-2) + G(t, 2,4)
#ft = abs(cos(t*1.25*2*pi/4))
#ft = pi-arccos(cos(t*3.8/2))
#ft[ft>1]=1
ft = 1-arccos(cos(t*10))/pi*2
f2 = (1-arccos(cos(t))/pi) * G(t,-pi, pi)
f3 = ft*f2
plt.plot(t, f3, lw=3)
plt.plot(t, f2, 'r--', lw=0.5)
plt.plot(t, ft, 'b--', lw=0.5)
#plt.arrow(0, 0.0, 0.0, 1.5, shape="full", lw=2, length_includes_head=True, head_width=0.1)
plt.xlabel("t")
plt.ylabel("f(t)")
plt.grid(True)
plt.grid(False)
plt.axis([min(t)-0.5, max(t)+0.5, -1.5, 1.5])
plt.tight_layout()
plt.show()
#------------------------------------------------------------
# END OF FILE : TEST1.PY
#============================================================
■ 相关文献链接:
● 相关图表链接:
相关文章
- 自然语言交流系统 phxnet团队 创新实训 个人博客 (十)
- 机房收费系统之思想性总结
- vcruntime140.dll 丢失64位系统
- Linux系统运维之路(运维纯属爱好,开发才是本职)
- 秒杀系统设计优化
- 用Python构建你自己的推荐系统
- Ubuntu 14.04 LTS中怎样解决系统设置残缺的问题
- 支付系统整体设计:整体架构设计以及注意要点(一)
- Linux系统针对网卡中断的优化处理
- windows 操作系统使用 pscp 拷贝文件到 Linux 系统遇到的错误消息
- 终于有阿里p9架构师分享出困扰我多年的Android系统开发进阶实战文档
- 如何在同一系统中启动多个 TOMCAT
- 内存管理 | 伙伴系统
- Linux系统root密码重置教程
- 2023年中职网络安全—Linux系统渗透提权2204解析(详细)
- 史上最详细教你制作“U盘启动盘”重装Windows10系统
- Win10注册表损坏进不了系统怎么U盘重装系统?
- 转:斯坦福大学心理学家:思维模式和信念系统,影响人生的底层代码
- Liquibase----Windows10系统安装Liquibase
- 【2023年第十一届泰迪杯数据挖掘挑战赛】C题:泰迪内推平台招聘与求职双向推荐系统构建 建模及python代码详解 问题二
- OSI的物理层 数据通信 数据通信系统的模型 信道的基本概念 编码方式 调制方法 香农公式 信噪比 硬核解析 图解刨析 宫廷解释
- 综合案例:利用ansible常用七种模块实现备份系统的部署及告警(十一)
- 2023年中职网络安全—Linux系统渗透提权2204解析(详细)