通信原理课程设计1new.docVIP

第一题 clear echo off t0=0.15; ts=0.001; fc=250; snr=20; fs=1/ts; df=0.3; t=(0:ts:t0); 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; [M,m,df1]=fftseq(m,ts,df); 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; dt=0.01; signal_power=sum(u.*u)*dt/(length(t).*dt); noise_power=signal_power/snr_lin; noise_std=sqrt(noise_power); n=noise_std*randn(1,length(u)); r=u+n; [N,n,df1]=fftseq(n,ts,df); %对噪声进行fft N=N/fs; %频率采样 [R,r,df1]=fftseq(r,ts,df); %对最后的信号进行fft R = R/fs; signal_power noise_power clf figure(1) plot(t,m(1:length(t))) xlabel(Time) title(The message signal) figure(2) subplot(2,1,1) plot(t,c(1:length(t))) xlabel(Time) title(The carrier) subplot(2,1,2) plot(t,u(1:length(t))) xlabel(Time) title(The modulated signal) figure(3) subplot(2,1,1) plot(f,abs(fftshift(M))) xlabel(frequency) title(spectrum of the message signal) subplot(2,1,2) plot(f,abs(fftshift(U))) title(spectrum of the modulated signal) xlabel(frequency) figure (4) subplot(2,1,1) %噪声时域显示 plot(t,n(1:length(t))) title(采样的噪声) subplot(2,1,2) %噪声频域显示 plot(f,abs(fftshift(N))) axis([-200 200 0 0.15]) title(噪声频谱) figure (5) subplot(2,1,1) %加噪调制信号时域显示 plot(t,r(1:length(t))) title(信号加噪声) axis([0 0.15 -3 3]) subplot(2,1,2) %加噪调制信号频域显示 plot(f,abs(fftshift(R))) title(信号加噪声频谱) axis([-500 500 0 0.15]) M文件1 fftseq: 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; M文件2 spower： function p=spower(x) p=(norm(x)^2/length(x)); 第二题 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.*(1/b)); clf plot(snrlin,ct) xlabel(s/n0) title(高斯信道