MATLAB幅度调制与解调源程序.pdfVIP

源程序

clc

T=0.001;%采样周期

fs=1/T;%采样频率

Tp=1;%观察时间

N=Tp/T;%采样点

n=-N/2:(N-1)/2;

t=n*T;

Yi=sin(40*pi*t);%调制信号

subplot(211);

plot(t,Yi);

xlabel(t);

ylabel(幅度);

title(调制信号);

axis([0,.5,-1,1]);

Yk=fft(Yi,2048);%对调制信号进行傅立叶变换

Yw=2*pi/N*abs(fftshift(Yk));

Fw=[-1024:1023]/2048*fs;

subplot(212);

plot(Fw,Yw);

xlabel(频率/hz);

ylabel(调制信号频谱幅度);

title(调制信号频谱);

axis([-40,40,0,3]);

figure

y1=cos(200*pi*t);%调制信号

subplot(211);

plot(t,y1);

xlabel(t);

ylabel(幅度);

.

title(载波信号);

axis([0,.1,-1,1]);

fm=20;

fs=1000;

N=512;

Y2=fft(y1,N);%对调制信号进行傅里叶变换

yw=abs(fftshift(Y2));

fw=[-255:256]/N*fs;

subplot(212);

plot(fw,yw);

xlabel(频率/hz);

ylabel(载波信号频谱幅度);

title(载波信号频谱);

axis([-200,200,0,100]);

clc

T=0.01;%采样周期

fs=1/T;%采样频率

Tp=10;

N=Tp/T;

n=0:N-1;

t=n*T;

Yi=cos(pi*t).*cos(20*pi*t);%调制后信号

subplot(211);

plot(t,Yi);

grid;

Yk=fft(Yi,2048);%对调制后信号进行傅立叶变换

Yw=2*pi/N*abs(fftshift(Yk));

Fw=[-1024:1023]/2048*fs;

subplot(212);

plot(Fw,Yw);

.

grid;

1

0.5

0

-0.5

-1

012345678910

2

1.5

1

0.5

0

-50-40-30-20-1001020304050

fs=1000;

t=0:1/fs:.8;

fc=100;

y=sin(pi*40*t);

x=modulate(y,fc,fs,amdsb-tc);%双边带调幅

subplot(211);

plot(t,x);

title(双边带调幅波形);

subplot(212);

p=fft(x,1024);

p1=fftshift(p);

w=0:511;

p2=abs(p1);