问题二：

2、单边带调制与解调：

载波频率：250KHZ；          

要求：

1、  原始信号及频谱

2、  调制信号及频谱

3、  解调信号及频谱

4、  要求在一幅图形中给出所用结果

程序：

dt=0.001;

fs=1/dt;

t=-0.2:dt:0.2;

fc=250;

m=sinc(200*t);

m=m.*m;

subplot(421);

plot(t,m);

axis([-0.2 0.2 -1.2 1.2]);

title('原始信号');

c=cos(2*pi*fc.*t);

d=sin(2*pi*fc.*t);

u=(m.*c-imag(hilbert(m)).*d)/2;

v=(m.*c+imag(hilbert(m)).*d)/2;

subplot(422);

plot(t,u,'r');

axis([-0.1 0.1 -1.2 1.2]);

title('上边带调制信号');

subplot(423);

plot(t,v);

title('下边带调制信号');

axis([-0.1 0.1 -1.2 1.2]);

F_m=dt*fft(m);

F_u=dt*fft(u);

F_v=dt*fft(v);

f=linspace(-fs/2,fs/2,length(m));

subplot(424);

plot(f,abs(fftshift(F_m)));

title('原始信号频谱');

subplot(425);

plot(f,abs(fftshift(F_u)));

title('上边带调制信号频谱');

subplot(426);

plot(f,abs(fftshift(F_v)));

title('下边带调制信号频谱');

y=v.*c;

F_y=dt*fft(y);

subplot(427);

plot(f,abs(F_y));

title('解调信号频谱');

H=zeros(size(f));

f_cutoff=100;

l=fs/length(m);

n=floor(f_cutoff/l);

H(1:n)=4*ones(1,n);

l1=length(m)-n;

H(length(m)-n:length(m)-l)=4*ones(l,n);

Z=F_y.*H;

z=fs*real(ifft(Z));

subplot(428);

plot(t,z);

axis([-0.1 0.1 -1.2 1.2]);

title('解调信号');

图形：

dt=0.001;

fs=1/dt;

t=-0.2:dt:0.2;

fc=250;

m=sinc(200*t);

m=m.*m;

subplot(421);

plot(t,m);

axis([-0.2 0.2 -1.2 1.2]);

title('原始信号');

c=cos(2*pi*fc.*t);

d=sin(2*pi*fc.*t);

u=(m.*c-imag(hilbert(m)).*d)/2;

v=(m.*c+imag(hilbert(m)).*d)/2;

subplot(422);

plot(t,u,'r');

axis([-0.1 0.1 -1.2 1.2]);

title('上边带调制信号');

subplot(423);

plot(t,v);

title('下边带调制信号');

axis([-0.1 0.1 -1.2 1.2]);

F_m=dt*fft(m);

F_u=dt*fft(u);

F_v=dt*fft(v);

f=linspace(-fs/2,fs/2,length(m));

subplot(424);

plot(f,abs(fftshift(F_m)));

title('原始信号频谱');

subplot(425);

plot(f,abs(fftshift(F_u)));

title('上边带调制信号频谱');

subplot(426);

plot(f,abs(fftshift(F_v)));

title('下边带调制信号频谱');

y=v.*c;

F_y=dt*fft(y);

subplot(427);

plot(f,abs(F_y));

title('解调信号频谱');

H=zeros(size(f));

f_cutoff=100;

l=fs/length(m);

n=floor(f_cutoff/l);

H(1:n)=4*ones(1,n);

l1=length(m)-n;

H(length(m)-n:length(m)-l)=4*ones(l,n);

Z=F_y.*H;

z=fs*real(ifft(Z));

subplot(428);

plot(t,z);

axis([-0.1 0.1 -1.2 1.2]);

title('解调信号');

（在最下面及右下角少一个频谱图形）

分析：

调制信号图形基本相同，但是幅度发生了变化，将幅度进行调整，及可得基本相同的调调制信号，根据由傅里叶变换得到的频谱进行分析，结果基本符合调制原理，但是由于各种噪音以及误差的影响，部分地方有失真