;%最小时刻Tmax=max(tau);%最大时刻T=ones(size(carry1));%一个数字信号‘1’dsig=kron(m,T);%数字信号波形subplot(7,1,1)%子图分割plot(tau,dsig)%画出数字信号的波形图title('数字信号');9gridon%添加网axis([TminTmax-0.21.2])%设置坐标范围subplot(7,1,2)%子图分割plot(tau,carry2)%画出载波波形title('载波波形');gridon%添加网格axis([TminTmax-1.2*A1.2*A])%设置坐标范围subplot(7,1,3)%子图分割plot(tau,ask)%画出调制后的波形title('调制后的波形');gridon%添加网格axis([TminTmax-1.2*A1.2*A])%设置坐标范围3.2数字信号的ASK相干解调%PartII数字信号的ASK相干解调sig_mul=ask.*carry2;%已调信号与载波信号相乘subplot(7,1,4)%子图分割plot(tau,sig_mul)%画出信号相乘后的波形Title('信号相乘后的波形');gridon%添加网格axis([TminTmax-0.21.2])[Ord,omega_c]=buttord(4*pi*f*0.6,4*pi*f*0.8,2,30,'s');%获 得Butterworth模拟低通原型滤波器的阶数及3dB截止频率[num,den]=butter(Ord,omega_c,'s');%由原型滤波器向实际滤波器转换,获得滤波器的分子,分母系数h=tf(num,den);%获得滤波器传递函数%滤波x=lsim(h,sig_mul,tau);%运用模拟滤波器对信号进行滤波subplot(7,1,5)%子图分割plot(tau,x)%画出滤波后的滤形