高斯噪声调频matlab,MATLAB仿真含高斯噪声条件下解调信号的时域图

clear all

dt=0.005; %设定步长 t=0:dt:3;

am=5; %调制信号幅度

fm=5; %调制信号频率

mt=am*cos(2*pi*fm*t); %生成调制信号

j_mt(1)=0;

for i=1:length(t)-1 %对调制信号求积分

j_mt(i+1)=j_mt(i)+mt(i)*dt;

end

fc=40;

ct=cos(2*pi*fc*t); %生成载波

kf=10; %调频灵敏度

sft=cos(2*pi*fc*t+kf*j_mt); %生成已调信号

figure(1) %绘制图形

subplot(3,1,1);plot(t,ct);

xlabel('t');title('载波时域图');

subplot(3,1,2);plot(t,mt);

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

subplot(3,1,3);plot(t,sft);

xlabel('t');title('已调信号时域图');

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

%*********高斯白噪声****************

sn=30; %设定信噪比

db=am^2/(2*(10^(sn/10))); %计算高斯白噪声方差

n=sqrt(db)*randn(size(t)); %生成高斯白噪声

nsfm=n+sft; %含高斯白噪声已调信