dsb.zip_DSB调制解调_dsb_模拟调制_线性系统_线性调制

%DSB 调制解调 %1）输入信号波形； %2）载波波形； %3）AM信号波形及AM包络； %4）相干解调后的信号波形； %5）AM信号功率谱 clear clc close all fs=100; % 采样率100Hz t=0:1/fs:10; % 设置仿真时间起止以及步长 fc=10; % 载波10Hz fm=1; % 调制信号1Hz 将输入信号幅值默认为1 mt=cos(2*pi*fm*t); % 定义调制信号 ct=cos(2*pi*fc*t); % 定义载波信号 psnt=ammod(mt,fc,fs); % DSB信号psnt envelope=abs(hilbert(psnt)); % 使用函数计算DSB信号的上包络 取负可得下包络 % DSB相干解调 tmp=psnt.*ct; [k,Wn,beta,ftype]=kaiserord([fc-fm fc+fm],[1 0],[0.05 0.01],fs); % 使用kaiser窗设计FIR滤波器 window=kaiser(k+1,beta); b=fir1(k,Wn,ftype,window,'noscale'); % 使用标准频率响应的加窗设计函数 jt=filter(b,1,tmp); % 滤波 % 功率谱 nfft=1024; window=hanning(256); % 构造hanning窗 noverlap=128; range='onesided'; % 求0 - Fs/2范围 [ps,f]=pwelch(psnt,window,noverlap,nfft); % PWELCH方法估计功率谱 f=f./pi*fs/2; % 将f的单位化为Hz %结果输出部分 subplot(5,1,1);%划分画图区间 plot(t,mt,'g');%画出调制信号波形 title('输入信号波形'); xlabel('Variable t'); ylabel('Variable mt'); axis([0 10 -1 1]); subplot(5,1,2); plot(t,ct,'b');%画出载波信号波形 title('输入载波波形'); xlabel('Variable t'); ylabel('Variable ct'); axis([2.5 5 -1 1]);%规定输出的坐标轴范围 subplot(5,1,3); plot(t,psnt,'r'); hold on; plot(t,envelope , 'g'); plot(t,-envelope, 'g'); hold off; title('已调信号波形及包络');%画出已调信号波形及包络 xlabel('Variable t'); ylabel('Variable psnt'); axis([2.5 5 -1 1]); subplot(5,1,4); plot(t,jt); title('相干解调输出');%画出相干解调后的波形图 xlabel('Variable t'); ylabel('Variable jt'); subplot(5,1,5); plot(f,ps); title('已调信号功率谱密度');%已调信号功率谱密度 xlabel('Variable f'); ylabel('Variable ps');