plot(t,v(1:length(t)));
axis([-0.2,0.2,-1.5,1.5]);
xlabel('时间');
title('上边带已调信号'); %上边带已调信号波形figure(3)
subplot(2,1,1);
plot(l,jt,'r');
axis([-1,1,-1000,1000])
xlabel('时间');
title('下边带解调信号'); %下边带解调信号波形subplot(2,1,2);
plot(l1,jt1);
axis([-1,1,-1000,1000])
xlabel('时间');
title('上边带解调信号'); %上边带解调信号波形figure(4)
V=fftshift(fft(v));
V0=abs(V);
V1=V0.^2;
df=0.5;
L=length(V);
f=-L/2*df:df:L/2*df-df;
subplot(2,1,1);
plot(f,V1);
title('下边带功率谱'); %下边带已调信号功率谱xlabel('f/HZ');
ylabel('V1');
U=fftshift(fft(u));
U0=abs(U);
U1=V0.^2;
df=0.5;
L=length(U);
f=-L/2*df:df:L/2*df-df;
subplot(2,1,2);
plot(f,U1);