fs=60000;

%频率为524Hz的声音极其波形 f1=524; t=0:1/fs:0.02; y1=sin(2*pi*f1*t); sound(y1,fs) plot(t,y1)

频率为524Hz的正弦信号声音波形

结果分析：

当正弦信号频率增大时，周期变短，音调变高。

（2）仿真程序：

f1=262;f2=294;f3=330;f4=262;f5=262;f6=294;f7=330;f8=262; f9=330;f10=349;f11=392;f12=392;f13=330;f14=349;f15=392; f16=392; fs=44100;

6

t=0:1/fs:0.5; y1=sin(2*pi*f1*t); y2=sin(2*pi*f2*t); y3=sin(2*pi*f3*t); y4=sin(2*pi*f4*t); y5=sin(2*pi*f5*t); y6=sin(2*pi*f6*t); y7=sin(2*pi*f7*t); y8=sin(2*pi*f8*t); y9=sin(2*pi*f9*t); y10=sin(2*pi*f10*t); y11=sin(2*pi*f11*t); y12=sin(2*pi*f12*t); y13=sin(2*pi*f13*t); y14=sin(2*pi*f14*t); y15=sin(2*pi*f15*t); y16=sin(2*pi*f16*t);

y=[y1 y2 y3 y4 y5 y6 y7 y8 y9 y10 y11 y12 y13 y14 y15 y16]; sound(y,fs)

仿真结果：

儿歌《两只老虎》的曲片段。

结果分析：

不同频率的无意义的信号在时域的叠加，形成了有意义的信号。

（3）仿真程序： time=4; Fs=44100;

[y,Fs]=audioread('st50_M2_7_3.wav'); p=audioplayer(y,fs); play(p);

t=(1:length(y))/Fs; plot(t,y,'s')

仿真结果：

7

播放音频：（要用兼容.wav格式播放器播放）

结果分析：

1．单位时间内时域图密集的部分，信号声音更加尖锐； 2. 单位时间内时域图稀疏的部分，信号声音更加低沉；

二、M3-3

1、问题重述

利用MATLAB提供的conv函数计算序列的卷积和，并验证卷积和的交换律、分配律与结合律。 ↓ ↓

x[k]={0.85,0.53,0.21,0.67,0.84,0.12},h[k]={0.68,0.37,0.83,0.52,0.71}

2、问题分析

conv函数可计算起点为k=0的两个序列的卷积。可利用此函数和卷积的性质，完成此题。

3、仿真程序与仿真结果 （1）仿真程序：

x=[0.85,0.53,0.21,0.67,0.84,0.12]; h=[0.68,0.37,0.83,0.52,0.71]; z=conv(x,h); N=length(z); stem(0:N-1,z);

仿真结果：

8

结果分析：

两个离散序列卷积和所得到的序列，其起点等于两个序列的起点之和，终点等于两个序列的终点之和。

（2）仿真程序：

%定义函数conv_m

function [y] = conv_m(x1,x2) y = conv(x1,x2);

x1 = [8, 5, 2, ];%定义x1序列 x2 = [6, 8, 1, ];%定义x2序列 x3 = [6, 3, 8, ];%定义x3序列

y11=conv_m(x1,x2); %交换律 y12=conv_m(x2,x1); %交换律

y21 = conv_m(conv_m(x1,x2),x3);%结合律 y22 = conv_m(x1,conv_m(x2,x3));%结合律

x4=x2+x3;

y31=conv_m(x1,x4);

y32=conv_m(x1,x2)+conv_m(x1,x3); %分配律

仿真结果：

结果分析：

9