[实验二带通信号的低通等效和随机信号的频带转换.docVIP

实验二 带通信号的低通等效和随机信号的频带转换 一、实验目的 根据内插与抽取原理，用Matlab 编程实现等效低通信号与频带信号的采样过程，比较二者采样率大小。掌握利用MATLAB实现连续信号采样、频谱分析和采样信号恢复的方法。 二、实验原理 三、 实验内容 对带通信号 . (1)画出该信号和它的幅度谱。 (2)求出该信号的解析信号，并画出它的幅度谱;(3)求出并画出该信号的包络。(4)分别假设f0 =100Hz和f0=50Hz，求该信号的低通等效，并画出它的幅度谱。若对该低通等效信号进行抽取，试设计取样频率，并重画抽取之后低通等效信号的频谱。(5)取该带通信号的低通等效部分： 试设计它的取样频率，并画出它的时域波形和它的幅度谱。 四、 实验步骤 五、 实验结果分析结论 六、实验程序 %通信系统建模与仿真 实验二 带通信号的低通等效和随机信号的频带变换 %主要参数：ts采样间隔 fs采样频率 %张希希 2014.6.17 clear all ts=0.002; %载波频率为100Hz，所以采样时间间隔ts=0.002 fs=1/ts; %采样频率 t=-1:ts:1; df=fs/length(t); %DFT的分辨率 f=-250:df:250-df; %生成频率矢量 x=2*sinc(20*t).*cos(2*pi*100*t+sinc(5*t)); X=fft(x)/fs; %求信号的频谱，因为原始信号是模拟信号，根据抽样定理，需要再计算出的fft后除以fs才能得到x的傅里叶变换 xa=hilbert(x); %求x的解析信号xa Xa=fft(xa)/fs; %求xa的频谱 subplot(2,1,1);plot(t,x);grid;title(信号x波形);xlabel(时间t); subplot(2,1,2);plot(f,fftshift(abs(X)));title(信号x幅度谱);xlabel(频率f); figure(1) subplot(2,1,1);plot(t,xa);title(信号xa包络);xlabel(时间t);%xa是复数，故画出其包络 subplot(2,1,2);plot(f,fftshift(abs(Xa)));title(信号xa幅度谱);xlabel(频率f); figure(2) plot(t,abs(x));title(信号x包络) f01=100; xl1=xa.*exp(-j*2*pi*f01*t); Xl1=fft(xl1)/fs; figure(3) subplot(2,1,1);plot(t,real(xl1));title(f01=100Hz时的低通信号同相分量);xlabel(时间t); subplot(2,1,2);plot(f,fftshift(abs(Xl1)));title(f01=100Hz时的低通信号幅度谱);xlabel(频率f); f02=50; xl2=xa.*exp(-j*2*pi*f02*t); Xl2=fft(xl2)/fs; figure(4) subplot(2,1,1);plot(t,real(xl2));title(f02=50Hz时的低通信号同相分量);xlabel(时间t); subplot(2,1,2);plot(f,fftshift(abs(Xl2)));title(f02=50Hz时的低通信号幅度谱);xlabel(频率f); wm=12.5; %低通等效信号带宽 wc=wm; %滤波器截止频率 Ts=pi/wm; %周期 临界采样 fs=1/Ts; %采样频率 ws=2*pi/Ts; %采样角频率 t1=-Ts*(length(t)-1)/2:Ts:Ts*(length(t)-1)/2; % f1=-fs*(length(f)-1)/2:fs:fs*(length(f)-1)/2; % n=-100:100; nTs=n*Ts; %时域采样点 xln=xa.*exp(-j*2*pi*f01*t1);%采样信号 Xln=fft(xln); figure(5) subplot(2,1,1);stem(t1,xln);axis([-30,30,-2,2]);title(低通等效信号采样后的信号);xlabel(n); subplot(2,1,2);plot(f1,fftshift(abs(Xln)));title(低通等效信号采样后的信号幅度谱);xl