9
1倍有效宽度结果：
5倍有效宽度结果：
（4）、
x(t)?e?tu(t)的频谱为X，根据公式0.1Xmax?X(Δf) Δf=9.95
当信号长度为2倍的有效长度时，取抽样频率f=10，50，100，200Hz。 f=10Hz时结果：
f=50Hz时结果：
f=100Hz时结果：
f=200Hz时结果：
11
综上所述，当信号抽样频率为100Hz左右时，误差较小。 （5）、 源代码：
[x1,fs,bits] = wavread('man.wav'); sound(x1,fs);
[X,f]=ctft1(x1,100,1024); subplot(211); plot (f,abs(X)); title('幅度'); subplot(212);
plot (f,angle(X)); title('相位');
[x1,fs,bits] = wavread('WOMAN.wav'); sound(x1,fs);
[X,f]=ctft1(x1,100,1024); subplot(211); plot (f,abs(X)); title('幅度'); subplot(212);
plot (f,angle(X)); title('相位');
结果：
女生：
男生：
区别：女生高频成分较多，而男声的高频成分较少。
（6）、
误差原因是用离散频谱近似等于连续的频谱，用和代替了积分。 3．幅度调制和连续信号的Fourier变换
本题研究莫尔斯码的幅度调制与解调。本题中信号的形式为
x(t)?m1(t)cos(2πf1t)?m2(t)sin(2πf2t)?m3(t)sin(2πf1t)
(1)
其中信号x(t)由文件ctftmod.mat定义(该文件在硬盘上)，可用命令Load ctftmod 将文件ctftmod.mat定义的变量装入系统内存。运行命令Load ctftmod后，装入系统的变量有
af bf dash dot f1 f2 t x
其中
bf af： 定义了一个连续系统H(s)的分子多项式和分母多项式。可利用freqs(bf,af,w)求出该系统的频率响应，也可用sys=tf(bf,af)得到系统的模型，从而用lsim求出信号通过该系统的响应。
dash dot： 给出了莫尔斯码中的基本信号dash和dot的波形 f1 f2： 载波频率 t: 信号x(t)的抽样点
x: 信号x(t)的在抽样点上的值
信号x(t)含有一段简单的消息。Agend 007的最后一句话是
The future of technology lies in ···
还未说出最后一个字，Agend 007就昏倒了。你(Agend 008)目前的任务就是要破解Agend 007的最后一个字。该字的信息包含在信号x(t)中。信号x(t)具有式(1)的形式。式中的调制频率分别由变量f1和f2给出，信号m1(t)，m2(t)和m3(t)对应于字母表中的单个字母，这个字母表已用国际莫尔斯码进行编码，如下表所示：
H ···· A ·? O ????? I ·· B ?··· P ·??· C ?·?· D ?·· E · F ··?· G ? ?· J ·????? K ?·? L ·?·· M ??? N ?· Q ? ?·? R ·?· S ··· T ? U ··? V ···? W ·?? X ?··? Y ?·?? Z ??··
(1) 字母B可用莫尔斯码表示为b=[dash dot dot dot]，画出字母B莫尔斯码波形； (2) 用freqs(bf,af,w)画出由bf和af定义的系统的幅度响应；
(3) 利用lsim求出信号dash通过由sys=tf(bf,af)定义的系统响应，解释你所获得的结果； (4)用解析法推导出下列信号的Fourier变换
m(t)cos(2πf1t)cos(2πf2t) m(t)cos(2πf1t)sin(2πf2t)
相关推荐：
loading