数字信号处理_MATLAB程序文件

例题3.1.1

Fm=10;Fc=100;Fs=500;

k=0:199;

t=k/Fs;

x=si n(2*pi*Fm*t);

y=x.*cos(2*pi*Fc*t);

Y=fft(y,256);

subplot(2,1,1);plot(y);

subplot(2,1,2);plot([-128:127],fftshift(abs(Y)));

实验所得的图像如下图所示：

实验内容1.实现抑制载波的幅度调制。已调信号y(t)二x(t)cosC c t),式中x(t)为调制信号;

cos(衣)为载波信号。此处可取x(t) =cos( m二t), c =80二rad /s，^^10_:

rad / s。

(1)分析调制信号x(t)的频谱，绘出其时域波形和频谱。

Matlab程序如下：

Fs=1000;

Fc=40;

N=1000;

wc=80*pi;

n=0:N-2;

t=n/Fs;

x=cos(10*pi*pi*t);

y=x.*cos(wc*t);

z仁fft(x,256);

subplot(2,1,1);plot(t,x);legend('x(t) 的时域波形');

subplot(2,1,2);plot([-128:127],fftshift(abs(z1)));lege nd('x(t) ');

的频谱实验所得图形如下：

Matlab 程序如下：

Fs=1000;

Fc=40;

N=1000;

n=0:N-2;

t=n/Fs;

x=cos(10*pi*pi*t);

y=modulate(x,Fc,Fs,'am');

z=fft(y,256);

subplot(2,1,1);plot(t,y);legend('y(t) 的时域波形');

subplot(2,1,2);plot([-128:127],fftshift(abs(z)));lege nd('y(t)

形')

(2)分析已调信号y(t)的频谱，绘出其时域波形和频谱。

的频域波

丫们的频域波形

(3)设计低通滤波器并应用filter 函数，实现信号同步解调

设计的低通滤波器为chebyl 型低通滤波器，技术参数为：

wp=0.2*pi,ws=0.6*pi,Ap=1Db,As=25dB

调制：

[N,wc]=cheb1ord(0.2*pi,0.6*pi,1,25,

' s ');

[num,den]=cheby1(N,1,wc, ' s ');

[b,a]=bili near(nu m,de n, 0.5);

x=cos(10*pi*pi*t);

y=x.*cos(80*pi*t);

w1=filter(b,a,y);

plot(w1);

Y=fft(w1);

plot(w1);axis([0 500 -1 1]);

5 o o

5 o -2 o.

o.

解调：

[N,wc]=cheb1ord(0.2*pi,0.6*pi,1,25,'s'); [nu m,de n]=cheby1(N,1,wc,'s'); [b,a]=bili near(nu m,de n, 0.5); x=cos(10*pi*pi*t);

y=x.*cos(80*pi*t);

w1=filter(b,a,y);

plot(w1);

Y=fft(w1);

Q=demod(Y,Fc,Fs,'am'); plot(t,Q);

axis([0.3 0.5 -0.1 0.1]);

2.实现含有载波的幅度调制。已调信号y(t) =[A • mx(t)]cos(「c t),选取A=1,m=0.5，观察调制结果，与抑制载波的幅度调制有什么不同？Matlab 程序：

Fm=10;Fc=100;Fs=1000;N=1000;

k=0:N-1;

t=k/Fs;

x=si n( 2*pi*Fm*t);

subplot(3,1,1);

plot(t,x);

xlabel('t(s)'); title('Modulate sig nal'); y=modulate(x,Fc,Fs,'am'); subplot(3,1,2);

plot(t,y);

xlabel('t(s)');

title('Modulateed sig nal(AM)'); y1=(1+0.5*x).*si n(1OO*pi*t); subplot(3,1,3);

plot(t,y1);

xlabel('t(s)');

title('Modulateed sig nal');

两种调制结果如下：

Modulateed signal (AM)

(3)实现含有载波的幅度调制y(t)二[A • mx(t)]cos(「c t)，调制信号为

(1)分析调制信号x(t)的频谱，绘出其时域波形和频谱。

(2)分析已调信号y(t)的频谱，绘出其时域波形和频谱。

①Matlab 程序为:

Fs=1000;

Fc=40;

N=1000;

wc=80*pi;

n=0:N-2;

t=0.5/4:0.001:1.5/4;

x=t;

z 仁 fft(x,256);subplot(4,1,1);plot(t,x);

grid on;legend('x(t)= t 的时域波形'); subplot(4,1,2);plot([-

128:127],fftshift(abs(z1)));

grid on;legend('x(t)= t 时的频谱');

y=(1+0.8*x).*cos(100*pi*t);

z=fft(y,256);

subplot(4,1,3);plot(t,y);

grid on;legend('y(t) 的时域波形');

3t 0 c x(t)二 t -t o 设 t o = 5s , 汁4 4 % E Yt 7色_ 4 X 载波为 cos( c t) = cos(100二t)，取 A=1,m=0.8。 3t o 4 4