郑州航院

nsfm=sfm;

for i=1:length(t)-1 %接受信号通过微分器处理 diff_nsfm(i)=(nsfm(i+1)-nsfm(i))./dt; end

diff_nsfmn = abs(hilbert(diff_nsfm)); %hilbert变换，求绝对值得到瞬时幅度（包络检波）

zero=(max(diff_nsfmn)-min(diff_nsfmn))/2; diff_nsfmn1=diff_nsfmn-zero;

subplot(3,1,3); %绘制无噪声条件下解调信号的时域图

plot((1:length(diff_nsfmn1))./1000,diff_nsfmn1./400,'r'); xlabel('时间t'); title('无噪声条件下解调信号的时域图');

%*****************************************

disp('按任意键可以看到原调制信号、小信噪比高斯白噪声条件下已调信号和解调信号已调信号的曲线') pause

%**************figure(4)****************** figure(4)

subplot(3,1,1);plot(t,mt); %绘制调制信号的时域图 xlabel('时间t');

title('调制信号的时域图');

db1=ac^2/(2*(10^(sn1/10))); %计算对应的小信噪比高斯白躁声的方差

n1=sqrt(db1)*randn(size(t)); %生成高斯白躁声

nsfm1=n1+sfm; %生成含高斯白躁声的已调信号（信号通 %过信道传输）

for i=1:length(t)-1 %接受信号通过微分器处理 diff_nsfm1(i)=(nsfm1(i+1)-nsfm1(i))./dt; end

diff_nsfmn1 = abs(hilbert(diff_nsfm1)); %hilbert变换，求绝对值得到瞬时幅度（包络检波）

zero=(max(diff_nsfmn)-min(diff_nsfmn))/2; diff_nsfmn1=diff_nsfmn1-zero; subplot(3,1,2);

plot(1:length(diff_nsfm),diff_nsfm); %绘制含小信噪比高斯白噪声已调信号的时域图

xlabel('时间t');

title('含小信噪比高斯白噪声已调信号的时域图');

subplot(3,1,3); %绘制含小信噪比高斯白噪声解调信号的时域图

plot((1:length(diff_nsfmn1))./1000,diff_nsfmn1./400,'r'); xlabel('时间t');