作业1 噪声时域信号分析。
fft的用途之一就是找出隐藏或淹没在噪声时域信号中信号的频率成分。首。
1)先建立试验数据(比如设数据的采样频率为1000hz,给出时间坐标区间从t=0到t=0.25,步长0.001秒);
2)其次生成一个包含两个频率成分的试验信号(比如50hz和120hz正弦波信号);
3)对这个信号加入随机噪声(可设其标准偏差为2),形成一个加噪信号y;
4)找到含噪声信号y的离散傅立叶变换(取它的fft即可),并求出信号的功率谱密度(它是不同频率所含能量的度量);
5)绘制出功率谱密度,讨论结果。
作业2 太阳黑子活动周期的分析。
太阳黑子的活动是周期的,大约每11年达到一个爆发高峰。试证明这一点。
设计提示:1、首先装载太阳黑子的数据,2、以横坐标表示年份,纵坐标表示黑子出现的数量,绘制wolfer图。
3、应用fft技术分析wolfer数,在复平面上可以直接绘制出由y给出的傅立叶系数的分布图。
4、绘制周期图(y的模的平方被定义为功率,功率与频率的关系曲线则被定义为周期图)。
附参考程序:
装载太阳黑子的数据程序如下:
load 这组数据记录了在过去的288年(1700-2024年)间每年太阳黑子出现的数量和大小的观测数据。
绘制wolfer图的程序如下:
year=sunspot(:,1); wolfer=sunspot(:,2);
plot(year,wolfer);
xlabel('years');ylabel(' sunspot data ')
title('sunspot data')
求wolfer数的fft程序如下:
y = fft(wolfer);
在复平面上绘制由y给出的傅立叶系数的分布图程序如下:
y(1)=[
plot(y,'ro')
title('fourier coefficients in the complex plane');
xlabel('real axis');
ylabel('imaginary axis');
set(gca,'ytick',[xtick',[
wolfer数的周期图程序如下:
n=length(y);
power = abs(y(1:n/2)).2;
nyquist = 1/2;
freq = 1:n/2)/(n/2)*nyquist;
plot(freq,power)
xlabel('cycles/year')
title('periodogram')
作业3 数字音频回声系统的模拟。
1. 应用背景。
**厅内听到的声音一般认为主要是由三部分组成:来自声源的直达波,经过墙壁有限次数反射的前期波和经过墙壁多次反射形成的后期波。由于声波传播路径的不同,这三种声音信号到达听众耳朵的先后顺序就有所不同,并且存在互相混叠的现象。
如果我们把整个**厅视为一个音响系统,比如把它看成一个对声源发出的声波信号进行处理的数字滤波器,则它的单位冲激响应反映的就是其回声特性的好坏。
2.问题描述。
数字回声滤波器可以看做是由一些基本的算法模块组成。这些算法模块包括延时相加模块,该模块模拟声音经墙壁反射,产生一定的延时,并与直接到达的声波相迭加,因此,可以用如下的差分方程描述这个模块的运算功能:
其中a是反射系数,表示声波经过反射后产生的衰减,d为反射造成的延时。上式对应的系统函数为。
它可以用离散滤波器模块实现。
由于已知这种延时相加模块的频谱特性在一些频率点是尖峰值,其结果将导致它对声音信号中的对应频率给予了强调,而其他频率成分则产生了相应的消弱。为补偿这种情况,我们引入另外一种具有全通性质的滤波算法,该算法的差分方程为:
对应的系统函数为。
它同样可以用离散滤波器模块实现。
我们通过适当的组合这两个基本算法模块,可以获得对实际**厅回声性能有较好模拟效果的数字滤波器系统。
3.解决方案。
有了上述两种算法模块,就可以通过对他们的适当组合来生成数字回响滤波器的simulink**模型。具体的组合方式根据实际建筑物的结构特征可以有多种,并随着相对与声源不同的位置而有所不同。因此,实际系统的**方案不是唯一的。
4.**研究。
下面是一个对实际效果模拟比较逼真的组合方案。它的框图加图10-18所示。
该系统中,左面的四个子系统是第一种延迟相加模块,滤波器参数分别为。
1. a = 0.75, b = 1, d = 29
2. a = 0.75, b = 0.9, d = 37
3. a = 0.75, b = 0.8, d = 44
4. a = 0.75, b = 0.7, d = 50
图10-18 数字音频回声系统的simulink**模型(echo_
右面的两个子系统是第二种延迟相加模块,他们的参数分别为。
1. a = 0.75, d = 27
2. a = 0.75, d = 31
输入信号是一个脉冲信号,输出信号直接送到“scone”示波器,结果如图10-19所示。
图10-19 数字音频回声系统的单位冲激响应。
从图10-19可以看出,所设计的数字音频回声模拟系统确实能较好地反映出**厅的回响性能:输入一个单位脉冲,输出信号在零时刻的取值代表直接到达的声音,一段时间间隔后的取值代表前期的反射声音信号,再往后杂乱密集的取值则代表后期的反射声音信号。
matlab设计作业
1.在同一个图形窗口内画出衰减 曲线及其包络线,取值范围是 0,6pi t 0 pi 100 6 pi y1 exp 3 t cos t 2 y2 exp 3 t t3 pi 0 9 9 y3 exp 3 t3 plot t,y1,r t,y2,b t3,y3,bo 2.画出所表示的三维曲面。x,y...
matlab作业设计
计算机 及应用 课程设计报告书。目录。一 设计思想。二 设计步骤。三 调试过程。四 结果分析。五 心得体会。1 编写matlab程序,建立该系统的数学模型。对系统进行阶跃响应 脉冲响应分析 求取给定输入u t sin t cos t 时的响应 绘制系统bode图和根轨迹。一 设计思想因为其为单输入单...
MATLAB期末设计作业
基于matlab的组合逻辑电路设计和 作者 马金龙。摘要 应用matlab语言设计数字电路的技巧和方法。本设计介绍了matlab simulink 软件的功能,叙述了相关的 数字电路 并以8线3线编码器 3线8线译码器 数据选择器等为例对其进行了 和分析。同时通过matlab汇编语言实现各逻辑电路的...