《信号与系统》matlab 作业 莫尔斯码 解调密码

(5)利用(4)中的结果，设计一个从x(t)中提取信号m1(t)的方案，画出m1(t)的波 形并确定其所代表的字母； Untiled2.m m1=lsim(bf,af,x.*cos(2*pi*f1*t),t); plot(t,m1);
d（- 。 。）
(6)对信号m2(t)和m3(t)重复(5)。请问Agent 008 The future of technology lies in ··· Untiled3.m m1=lsim(bf,af,x.*cos(2*pi*f1*t),t); subplot(3,1,1); plot(t,m1); m1=lsim(bf,af,x.*sin(2*pi*f2*t),t); subplot(3,1,2); plot(t,m1); m1=lsim(bf,af,x.*sin(2*pi*f1*t),t); subplot(3,1,3); plot(t,m1);
G − −·
N −·
U ··−
先调用函数得到ctftmod.m中的内容。 >>clear >> load ctftmod >> whos
Name af bf dash dot f1 f2 t x
Size 1x6 1x6 1x2000 1x2000 1x1 1x1 1x8000 1x8000
Bytes Class Attributes 48 double 48 double 16000 double 16000 double 8 double 8 double 64000 double 64000 double
(4)用解析法推导出下列信号的Fourier变换
m ( t ) cos( 2 π f 1t ) cos( 2 π f 2 t ) m ( t ) cos( 2 π f 1t ) sin( 2 π f 2 t ) m ( t ) sin( 2 π f 1t ) sin( 2 π f 2 t )
假设m1(t)对应M1（w）， m2（t）对应M2（w），m3（t）对应M3（w） F(1）=0.5M1(w)+0.25[M1(w+4pif1)+M1(w-4pif1)] F(2)=0.25j[M2(w+2pi(f2-f1))-M2(w-2pi(f2-f1))+M2(w+2pi(f2+f1))-M2(w-2pi(f1+f2))] F(3)=0.25j[M1(w+4pif1)+M1(w-4pif1)]
d（- 。 。） s（。 。 。） p（。 - - 。）
数字信号处理
来自百度文库
t: 信号x(t)的抽样点 x: 信号x(t)的在抽样点上的值 信号x(t)具有式(1)的形式。式中的调制频率分别由变量f1和f2给出，信号m1(t)， m2(t)和m3(t)对应于字母表中的单个字母
字母表已用国际莫尔斯码进行编码，如下表所示：
A ·− B −··· C −·−· D −·· E · F ··−· H ··· I ·· J ·− − − K −·− L ·−·· M −− O −−− P ·−−· Q − −·− R ·−· S ··· T − V ···− W ·−− X −··− Y −·−− Z −−··
af和bf是一个低通滤波器的频率响应。 f1和f2是载波频率： f1=200hz f2=400hz
>> plot(dash)
>> plot(dot)
（1）字母B可用莫尔斯码表示为b=[dash dot dot dot]，画出字母B莫尔斯码波形；
>> b=[dash dot dot dot] >> plot(b)
(2) 用freqs(bf,af,w)画出系统的幅度响应；
>> freqs(bf,af)
从图中看出：这是一个低通滤波器
(3) 利用lsim求出信号dash通过由sys=tf(bf,af)定义的系统响应，解释你所获得的结果； Untiled1.m ydash=lsim(bf,af,dash,t(1:length(dash))); ydot=lsim(bf,af,dot,t(1:length(dot))); subplot(2,1,1); plot(t(1:length(dash)),dash,t(1:length(dash)),ydash,'--'); legend('dash','ydash'); subplot(2,1,2); plot(t(1:length(dot)),dot,t(1:length(dot)),ydot,'--'); legend('dot','ydash'); 从信号的波形上看，原信号和输出信号 的波形基本一致，输出信号在幅度上有 些而改变，在时间上有些延迟。这个结 果说明他们都是低频的，fft变换的频率 范围在低通滤波器通带之内。
本题研究莫尔斯码的幅度调制与解调。本题中信号的形式为
x(t) = m1(t) cos(2πf1t) + m2(t) sin(2πf2t) + m3(t) sin(2πf1t)
其中信号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： 载波频率