DSB-SC信号的处理

				%DSB-SC信号的处理
				close all
				clear all
				%采样点数的设置
				k=input('采样点数N=2^k，请输入k的值，默认为15:');
				if isempty(k),
				    k=15;
				end;
				N=2^k;%采样点数
				dt=0.01;%时域取点间隔
				df=1.0/(N*dt);%频域取点间隔
				T=N*dt;%时域截断区间
				Bs=N*df/2;%频域截断区间
				t=linspace(-T/2,T/2,N);%生成时域采样点
				f=linspace(-Bs,Bs,N);%生成频域采样点
				fm=1;%基带信号频率
				fc=6*fm;%载波频率
				m=sin(2*pi*fm*t);%基带信号的产生
				c=sin(2*pi*fc*t);%载波信号的产生
				s=m.*c;%DSB-SC信号的产生
				S=t2f(s,dt);%傅立叶变换
				figure(1)%开启一个编号为1的绘图窗口
				%设置窗口1左上角的位置在距屏幕左侧100像素、下侧300像素的地方，长为400像素，宽为300像素
				set(1,'position',[100,300,400,300])
				plot(t,s)%画时域图
				%设置窗口的各种参数
				title('时域图')
				xlabel('t(ms)')
				ylabel('s(t)')
				axis([-4,4,min(s),max(s)])
				grid on
				figure(2)%开启一个编号为2的绘图窗口
				%设置窗口2左上角的位置在距屏幕左侧500像素、下侧300像素的地方，长为400像素，宽为300像素
				set(2,'position',[500,300,400,300])
				plot(f,abs(S))%画频域图
				title('频域图')
				xlabel('f(kHz)')
				ylabel('S(f)')
				axis([-10,10,min(abs(S)),max(abs(S))])
				grid on