matlab实验报告 模拟调制解调

实验报告
姓名：李鹏博实验名称：模拟调制解调
学号：2011300704 课程名称：数字信号处理
班级：03041102 实验室名称：航海西楼303
组号： 1 实验日期：2014.06.20
一、实验目的、要求
1．掌握掌握模拟调制以及对应解调方法的原理。

2．掌握模拟调制解调方法的计算机编程实现方法，即软件实现。

二、实验原理
调制的目的是把要传输的模拟信号或数字信号变换成适合信道传输的信号，这就意味着把基带信号（信源）转变为一个相对基带频率而言频率非常高的代通信号。

三、实验环境
PC机，Windows2000，office2000，Matlab6.5以上版本软件
四、实验过程、数据记录、分析及结论
实验过程
1.根据参数产生消息信号s和载波信号。

2.编程实现调制过程。

调用函数y=ammod(s,Fc,Fs)完成幅度调制，y=fmmod(s,Fc,Fs,
FREQDEV) 完成频率调制，y=pmmod(s,Fc,Fs, PHASEDEV) 完成相位调制。

3.编程实现信号的传输过程。

产生白噪声noise，并将其加到调制信号序列。

4.编程实现信号的解调。

调用函数x=amdemod(y,Fc,Fs)完成幅度调制信号的解调，
x=fmdemod(y,Fc,Fs, FREQDEV) 完成频率调制信号的解调，x=pmdemod(y,Fc,Fs, PHASEDEV) 完成相位调制信号的解调。

数据记录
消息信号s幅度调制和FFT 消息信号s频率调制和FFT
消息信号s相位调制和FFT 幅度调制信号加入白噪声及FFT
频率调制信号加入白噪声及FFT 相位调制信号加入白噪声及FFT
三种调制方式后经过滤波如上图解调后信号及其FFT
解调后信号及其FFT 解调后信号及其FFT
五、讨论
对于实验中的FREQDEV 、PHASEDEV 含义不清楚，具体数值不知道该如何确定。

通过这次的实验，我们对信息和通信系统有了更进一步的认识，尤其是在系统设计方面，尽管是非常基础的调制与解调的传输，也是经过若干设备协同工作，才能保证信号有效传输，而小到仅仅是一个参数，都有可能导致整个系统无法正常运行。

另一方面，我们通过本次的实验，着实领教了MATLAB强大的功能和实力。

实验所得这些经验虽然并不高深，但是对于刚入门的初学者来说，对以后步入专业领域进行设计或研发无疑具有重大的意义。

附：
close all
clear all
f1 = 10;
Fc = 8000;
Fs = 40000;
t = [0:Fs]/Fs;
f = [-Fs/2:Fs/2];
f0 = [0:Fs];
FREQDEV = 10; PHASEDEV = 5;
s = sin(2*pi*f1*t);
ya = ammod(s,Fc,Fs);
fya=fftshift(fft(ya));
yf = fmmod(s,Fc,Fs, FREQDEV);
fyf=fftshift(fft(yf)); yp = pmmod(s,Fc,Fs, PHASEDEV);
fyp=fftshift(fft(yp));
figure(1)
subplot(3,1,1)
plot(t,s)
subplot(3,1,2)
plot(t,ya)
subplot(3,1,3)
plot(f,abs(fya))
figure(2)
subplot(3,1,1)
plot(t,s)
subplot(3,1,2)
plot(t,yf)
subplot(3,1,3)
plot(f,abs(fyf))
figure(3)
subplot(3,1,1)
plot(t,s)subplot(3,1,2)
plot(t,yp)
subplot(3,1,3)
plot(f,abs(fyp))
%加入噪声
Ya =
awgn(ya,1,'measured');
fYa=fftshift(fft(Ya));
Yf =
awgn(yf,1,'measured');
fYf=fftshift(fft(Yf));
Yp =
awgn(yp,1,'measured');
fYp=fftshift(fft(Yp));
figure(4)
subplot(2,1,1)
plot(t,Ya)
subplot(2,1,2)
plot(f,abs(fYa))
figure(5)
subplot(2,1,1)
plot(t,Yf)
subplot(2,1,2)
plot(f,abs(fYf))
figure(6)
subplot(2,1,1)
plot(t,Yp)
subplot(2,1,2)
plot(f,abs(fYp))
%滤波
b = fir1(1024,
[2*7900/Fs,2*9100/Fs]);
h1 = filter(b,1,Ya);
h2 = filter(b,1,Yf);
h3 = filter(b,1,Yp);
figure(10)
subplot(311)
plot(t,h1)
subplot(312)
plot(t,h2)
subplot(313)
plot(t,h3)
%解调
yya = amdemod(h1,Fc,Fs)
fyya =
fftshift(fft(yya));
yyf =
fmdemod(h2,Fc,Fs,FREQDE
V)
fyyf =
fftshift(fft(yyf));
yyp =
pmdemod(h3,Fc,Fs,PHASED
EV)
fyyp =
fftshift(fft(yyp));
figure(7)
subplot(3,1,1)
plot(t,ya)
subplot(3,1,2)
plot(t,yya)
subplot(3,1,3)
plot(f,abs(fyya))
figure(8)
subplot(3,1,1)
plot(t,yf)
subplot(3,1,2)
plot(t,yya)
subplot(3,1,3)
plot(f,abs(fyya))
figure(9)
subplot(3,1,1)
plot(t,yp)
subplot(3,1,2)
plot(t,yya)
subplot(3,1,3)
plot(f,abs(fyya))。