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

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

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);

y1=xcorr(y,'unbiased'); % 计算序列的自相关函数

y2=abs(fft(y1,nfft));

index=0:round(nfft/2-1);

k=index*Fs/nfft;

y3=10*log10(y2(index+1));

c1=xcorr(c,'unbiased');

c2=abs(fft(c1,nfft));

index2=0:round(nfft/2-1);

k2=index2*Fs/nfft;

c3=10*log10(c2(index2+1));

u1=xcorr(u,'unbiased');

u2=abs(fft(u1,nfft));

index3=0:round(nfft/2-1);

k3=index*Fs/nfft;

u3=10*log10(u2(index3+1));

subplot(3,2,1),plot(t,y),grid; % 画图显示结果

subplot(3,2,2),plot(t,c),grid;

subplot(3,2,3),plot(t,u),grid;

subplot(3,2,4),plot(k,y3),grid;

subplot(3,2,5),plot(k2,c3),grid;

subplot(3,2,6),plot(k3,u3),grid;

2.

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

t=0:0.0001:0.5;

f=130;

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

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

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