资源描述:
《北邮通原软件实验》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、实验一实验目的:假设基带信号为m(t)=sin(2000πt)+2cos(1000πt),载波频率为20kHz,请仿真出AM,DSB-SC,SSB信号,观察已调信号的波形和频谱。1.AM信号:(1)信号的表达式(3)流程图(2)源代码%AM信号的产生fs=800;%采样频率KHzT=200;%截短时间msN=T*fs;%采样点数dt=1/fs;t=[-T/2:dt:T/2-dt];df=1/T;f=[-fs/2:df:fs/2-df];fm=1;%kHzfc=20;%kHzm=sin(2*pi*fm*t)+2*cos(1*fm*pi*t);s=(1+0.3*m).*cos(2*p
2、i*fc*t);%AM信号S=t2f(s,fs);figure(1)plot(f,abs(S1))title('AM信号频谱')xlabel('f')ylabel('S(f)')axis([-25,25,0,max(abs(S1))]);%xset('window',2)figure(2)plot(t,s1)title('AM信号波形')xlabel('t')ylabel('s(t)')axis([-3,3,-3,3]);(4)实验结果2.DSB-SC信号(1)信号的产生和表达式(2)流程图(3)源代码fs=800;%KHzT=200;%msN=T*fs;dt=1/fs;t=[-
3、T/2:dt:T/2-dt];df=1/T;f=[-fs/2:df:fs/2-df];fm=1;%kHzfc=20;%kHzm=sin(2*pi*fm*t)+2*cos(1*fm*pi*t);s=m.*cos(2*pi*fc*t);%DSB-SC信号S=t2f(s,fs);figure(1)plot(f,abs(S2))title('DSB-SC信号频谱')xlabel('f')ylabel('S(f)')axis([-25,25,0,max(abs(S2))]);figure(2)plot(t,s2)title('DSB-SC信号波形')xlabel('t')ylabel('s
4、(t)')axis([-1,4,-3,3]);(4)实验结果2DSC-SB频谱3.SSB信号(1)信号的产生和表达式(2)流程图(3)源代码:%SSB信号的产生fs=800;%KHzT=200;%msN=T*fs;dt=1/fs;t=[-T/2:dt:T/2-dt];df=1/T;f=[-fs/2:df:fs/2-df];fm=1;%kHzfc=20;%kHzm=sin(2*pi*fm*t)+2*cos(1*fm*pi*t);M=t2f(m,fs);MH=-j*sign(f).*M;%在频域进行希尔伯特变换mh=real(f2t(MH,fs));%希尔伯特变换后的信号s=m.*c
5、os(2*pi*fc*t)-mh.*sin(2*pi*fc*t);%SSBsignalS=t2f(s,fs);figure(1)plot(f,abs(S3))title('SSB信号频谱')xlabel('f')ylabel('S(f)')axis([-25,25,0,max(abs(S3))])figure(2)plot(t,s3)title('SSB信号波形')xlabel('t')ylabel('s(t)')axis([0,6,-3,3])(4)实验结果实验二实验目的:假设基带信号为m(t)=sin(2000πt)+2cos(1000πt)+4sin(500πt+π/3),
6、载波频率为40kHz,仿真产生FM信号,观察波形与频谱,并与卡松公式作对照。FM的频率偏移常数是5kHz/V。(1)信号表达式(2)流程图(3)源代码fs=800;%kHzT=16;%msN=T*fs;dt=1/fs;t=[-T/2:dt:T/2-dt];df=1/T;f=[-fs/2:df:fs/2-df];fm=1;%kHzKf=5;%kHz/Vfc=40;%kHzm=sin(2*pi*fm*t)+2*cos(1*pi*fm*t)+4*sin(0.5*pi*fm*t+pi/3);phi=2*pi*Kf*cumsum(m)*dt;%求相位s=cos(2*pi*fc*t+phi)
7、;%s(t)S=t2f(s,fs);figure(1)plot(f,abs(S).^2)title('FM信号功率谱')xlabel('f')ylabel('S(f)')axis([-80,80,0,max(abs(S).^2)]);%功率谱密度为
8、S
9、^2figure(2)plot(f,abs(S))title('调制信号频谱')xlabel('f')ylabel('S(f)')axis([-80,80,0,max(abs(S))]);figure(3)plot(t,s)tit