dsb_sc.rar_DSB载波抑制_Dsb时域频域_dsb_调幅仿真信号

%DSB-SC信号时域、频域波形仿真 close all %关闭所有的窗口 clear all %清除所有的变量 %要求输入k,从而计算总取样点数N=2^k k=input('取样点数=2^k,k=[18] '); if isempty(k), k=18; %默认K=18 end; N=2^k; %初始化时域分辨率dt（单位ms）， %频域分辨率df（单位kHz） %并满足df*dt=1/N dt=0.01; df=1.0/(N*dt); %计算时间宽度T（单位ms） %系统带宽Bs（单位kHz） T=N*dt; Bs=N*df/2; %产生时间取样点距阵t %和频域取样点距阵f t=[-T/2+dt/2:dt:T/2]; f=[-Bs+df/2:df:Bs]; %基带信号fm=1kHz，载波fc=6kHz fm=1; fc=6*fm; %基带信号m，载波c m=sin(2*pi*fm*t); c=sin(2*pi*fc*t); s=m.*c; %做傅立叶变换，产生频域抽样S S=t2f(s,dt); %计算出功率谱P P=S.*conj(S)/T; %画时域图，开窗口1 figure(1) %设置窗口位置和大小 set(1,'position',[10,50,300,200]) plot(t,s,'b',t,m,'r:',t,-m,'r:'); %给图形打上格线 grid on %设置坐标范围 axis([-2^10*dt/10,2^10*dt/10,min(s)*1.1,max(s)*1.1]) %axis([-1,1,min(s)*1.1,max(s)*1.1]) %给X轴、Y轴加标签 xlabel('t(ms)') ylabel('s(t)') title('DSB-SC信号时域波形') %画频谱曲线，开窗口2 figure(2) set(2,'position',[350,50,300,200]) plot(f,abs(S),'b') grid on axis([-20,20,min(abs(S)),max(abs(S))*1.1]) xlabel('f(kHz)'); ylabel('S(f)'); title('DSB-SC信号频域波形')