数字信号处理_MATLAB程序文件.docxVIP

.word格式 .word格式, 专业 专业.专注 例题3.1.1 Fm=10;Fc=100;Fs=500; k=0:199; t=k/Fs; x=si n(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.实现抑制载波的幅度调制。已调信号y(t)二x(t)cosC ct),式中x(t)为调制信号; cos(衣)为载波信号。此处可取 x(t) =cos( m二t), c =80二 rad /s， ^^10_: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); z仁 fft(x,256); 的频谱subplot(2,1,1);plot(t,x);legend(x(t) 的时域波形); 的频谱 subplot(2,1,2);plot([-128:127],fftshift(abs(z1)));lege nd(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); y=modulate(x,Fc,Fs,am); z=fft(y,256); subplot(2,1,1);plot(t,y);legend(y(t) 的时域波形); 的频域波subplot(2,1,2);plot([-128:127],fftshift(abs(z)));lege nd(y(t) 的频域波 形) 5 o o 5 o -2 o. o. 丫们的频域波形 (3)设计低通滤波器并应用filter 函数，实现信号同步解调 设计的低通滤波器为chebyl型低通滤波器，技术参数为： wp=0.2*pi,ws=0.6*pi,Ap=1Db,As=25dB 调制： [N,wc]=cheb1ord(0.2*pi,0.6*pi,1,25, s); [num,den]=cheby1(N,1,wc, s); [b,a]=bili near(nu m,de n, 0.5); x=cos(10*pi*pi*t); y=x.*cos(80*pi*t); w1=filter(b,a,y); plot(w1); Y=fft(w1); plot(w1);axis([0 500 -1 1]); 调制后的波形如下图: 解调： [N,wc]=cheb1ord(0.2*pi,0.6*pi,1,25,s); [nu m,de n]=cheby1(N,1,wc,s); [b,a]=bili near(nu m,de n, 0.5); x=cos(10*pi*pi*t); y=x.*cos(80*pi*t); w1=filter(b,a,y); plot(w1); Y=fft(w1); Q=demod(Y,Fc,Fs,am); plot(t,Q); axis([0.3 0.5 -0.1 0.1]); 解调后的波形如下图: 2.实现含有载波的幅度调制。已调信号y(t) =[A ? mx(t)]cos(「ct),选取 A=1,m=0.5，观察调制结果，与抑制载波的幅度调制有什么不同？ Matlab 程序： Fm=10;Fc=100;Fs=1000;N=1000; k=0:N-1; t=k/Fs; x=si n( 2*pi*Fm*t); subplot(3,1,1); plot(t,x); xlabel(t(s)); title(Modulate sig nal); y=modulate(x,Fc,Fs,am); subplot(3,1,2); plot(t,y); xlabel(t(s)); title(Modulateed sig nal(AM)); y1=(1+0.5*x).*si n(1OO*pi*t); subplot(3,1,3); plot(t,y1); xlabel(t(s)); title(Modulateed sig nal); 两种调制结果如下： Modulateed signal (AM) (3)实现含有载波的幅度调制y(t)二[A ? mx(t)]cos(「ct)，调制信号为 3t0 c x(t)二 t -to 设 to = 5s , 汁4 4 % E Yt