《MATLAB与通信仿真》实验1.doc

MATLAB与通信仿真实验1一、目的和要求：要求学生在熟练掌握MATLAB仿真使用的基础上，学会通信仿真系统的基本设计与调试。并结合通信原理的知识，对通信仿真系统进行性能分析。二、实验环境PC机、Matlab三、具体内容及要求信号定义为：，其中A的取值为学生学号的末两位乘以0.1。的取值为学生学号的末位乘以0.05（s）。该信号以常规幅度调制方式调制，其中调制指数为0.8，载波，其中取值为学生学号的末三位（Hz），所得已调信号为。l 试用Matlab编程实现：画出、和的波形和功率谱密度l 并给出常规幅度调制信号的非相干解调。l 若将一个高斯噪声加到已调信号上，所得的SNR是10dB，求噪声功率。分别画出该高斯噪声的波形和功率谱密度，高斯噪声加已调信号的波形和功率谱密度。(计算一个周期的波形即可)1.function y=cmccc(t) % M函数定义t=0:0.0001:0.5;f=130;nfft=1024;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 endend 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取值所在范围 y(m)=3; % 计算分段函数值 else if t(m)=1/3 y(m)=-6; else y(m)=0; end endend 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 endend 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); % 计算噪声的功率谱密度index=0:round(nfft/2-1); k=index*Fs/nfft; A2=10*log10(A1(index+1); B=xcorr(sn,unbiased); % 计算加噪声的已调信号的自相关函数 B1=abs(fft(B,nfft); % 计算加噪声的已调信号的功率谱密度index2=0:round(nfft/2-1); k2=index2*Fs/nfft; B2=10*log10(B1(index2+1); subplot(2,2,1),plot(t,noise),title(高斯白噪声); % 画图显示结果subplot(2,2,2),plot(k,A