《MATLAB与通信仿真》实验1(3)

MATLAB仿真,常规幅度调制方式,画出波形和功率谱密度,非相干解调,高斯噪声

y(m)=3; % 计算分段函数值

else if t(m)<=1/3

y(m)=-6;

else y(m)=0;

end

end

end

c=cos(2*pi*f*t);

u=c.*(2*1/0.8+y);

bl=abs(hilbert(u)); % 利用函数计算包络检波后的信号

% 画图显示各信号

subplot(3,1,1),plot(t,u),title('已调信号');

subplot(3,1,2),plot(t,bl),title('非相干解调后的信号');

subplot(3,1,3),plot(t,y),title('原始信号');

3.

function y=gszs(t) % M函数定义

t=0:0.0001:0.5;

f=130;

nfft=1024;

SNR=10;

Fs=10000;

leng=length(t); % 计算向量x的长度

for m=1:leng % 计算函数值

if t(m)<=1/6 % 判断x取值所在范围

y(m)=3; % 计算分段函数值

else if t(m)<=1/3

y(m)=-6;

else y(m)=0;

end

end

end

c=cos(2*pi*f*t);

u=c.*(2*1/0.8+y);

sn=awgn(u,SNR);

noise=sn-u;

A=xcorr(noise,'unbiased'); % 计算噪声的自相关函数

A1=abs(fft(A,nfft)); % 计算噪声的功率谱密度