通信原理程设计全.docVIP

第一题 clear echo off t0=0.15; ts=0.001; fc=250; snr=20; fs=1/ts; t=(0:ts:t0); dt=0.01; df=0.3; snr_lin=10^(snr/10); m=[ones(1,t0/(3*ts)),-2*ones(1,t0/(3*ts)),zeros(1,t0/(3*ts)+1)]; c=cos(2.*pi.*fc.*t); u=[2+0.85*m].*c; clf figure(1) %（1）m(t)和u(t)的波形 subplot(2,1,1) plot(t,m(1:length(t))) xlabel(Time) title(基带信号m(t)) subplot(2,1,2) plot(t,u(1:length(t))) xlabel(Time) title(调制信号u(t)) [M,m,df1]=fftseq(m,ts,df); %（2）m(t)和u(t)的频谱图 M=M/fs; [U,u,df1]=fftseq(u,ts,df); U=U/fs; [C,c,df1]=fftseq(c,ts,df); f=[0:df1:df1*(length(m)-1)]-fs/2; figure(2) subplot(2,1,1) plot(f,abs(fftshift(M))) xlabel(frequency) title(基带信号频谱M(f)) subplot(2,1,2) plot(f,abs(fftshift(U))) xlabel(frequency) title(调制信号频谱U（f)) signal_power=sum(u.*u)*dt/(length(t).*dt); %（3）周期T0=t0，求已调信号的功率 signal_power noise_power=signal_power/snr_lin; %（4）噪声功率 noise_power noise_std=sqrt(noise_power); n=noise_std*randn(1,length(u)); r=u+n; figure (5) %（5） 噪声n(t)和已调信号加噪声（u(t)+ n(t)）的波形； subplot(2,1,1) %噪声时域显示 plot(t,n(1:length(t))) title(噪声信号) subplot(2,1,2) %加噪调制信号时域显示 plot(t,r(1:length(t))) title(信号加噪声信号) [N,n,df1]=fftseq(n,ts,df); %对噪声进行fft N = N/fs; %频率采样 [R,r,df1]=fftseq(r,ts,df); %对最后的信号进行fft R = R/fs; figure (6) %（6）噪声n(t)和已调信号加噪声（u(t)+ n(t)）的频谱图； subplot(2,1,1) %噪声频域显示 plot(f,abs(fftshift(N))) axis([-200 200 0 0.15]) title(噪声信号频谱) axis([0 0.15 -3 3]) subplot(2,1,2) %加噪调制信号频域显示 plot(f,abs(fftshift(R))) title(信号加噪声信号频谱) axis([-500 500 0 0.15]) M文件 function [M,m,df]=fftseq(m,ts,df) fs=1/ts; n1=fs/df; n2=length(m); n=2^(max(nextpow2(n1),nextpow2(n2))); M=fft(m,n); m=[m,zeros(1,n-n2)]; df=fs/n; 第二题 clear echo off b=3000; sn0=-20; sn1=30; ns=0.01; N=(sn0:ns:sn1); snrlin=10.^(N/10); ct=b.*log2(1+snrlin.*(