数字信号处理-MATLAB程序

例题3.1.1

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

k=0:199;

t=k/Fs;

x=sin(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.实现抑制载波的幅度调制。已调信号()()cos()c y t x t t ω=，式中()x t 为调制信号；cos()c t ω为载波信号。此处可取()cos(),80m c x t t ωπωπ== /rad s ，10m ωπ= /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);

z1=fft(x,256);

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

subplot(2,1,2);plot([-128:127],fftshift(abs(z1)));legend('x(t)的频谱'); 实验所得图形如下：

（2）分析已调信号()

y t的频谱，绘出其时域波形和频谱。Matlab程序如下：

Fs=1000;

Fc=40;

N=1000;

n=0:N-2;

t=n/Fs;

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