数学建模社区-数学中国

标题: 我的困惑，强人请帮忙 [打印本页]

作者: 紫藤恋晴天 时间: 2012-8-15 22:47
标题: 我的困惑，强人请帮忙
这俩天的程序我有点抓狂，刷的题我也很无耐，谁有matlab软件的看到后，帮忙运行一下，最后的一个
“3、  解调信号及频谱”中的频谱没弄出来，希望能帮忙的友人帮忙看看哪有错误吧

问题二：

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('解调信号');