白噪聲信道實驗

1、實驗二 白噪聲信道實驗一、代碼如下N=4096;dt=0.01; %采樣間隔t1=-N/2:N/2*0.01; %時間分布 Ts=0.01*N; %采樣時長df=1/Ts; %頻率分辨率fl=1/dt; %頻率范圍f=-fl/2:df:fl/2; %頻率分布%1、產(chǎn)生并分析理想低通高斯白噪聲的時頻特性，繪制其時頻波形。n0=2;fh=10;Rn1=n0*fh*sinc(2*fh*t1);figure(1);subplot(2,1,1);plot(t1,Rn1);axis(-1 1 -5 20);title('低通理想白噪聲的時域'); xlabel('t(s)'

2、)grid on;subplot(2,1,2);Rw1=fftshift(fft(Rn1,N+1);plot(f,abs(Rw1);axis(-30 30 0 100);title('低通理想白噪聲的頻域'); xlabel('w(Hz)')grid on;B=10;fc=20;%2、產(chǎn)生并分析理想帶通高斯白噪聲的時頻特性，繪制其時頻波形。Rn2=n0*B*sinc(B*t1).*cos(2*pi*fc*t1);%ts=2*B*sinc(B*nt).*cos(2*pi*fs*nt);figure(2);subplot(2,1,1);plot(t1,Rn2);ax

3、is(-1 1 min(Rn2) max(Rn2);grid on;Rf2=fftshift(fft(Rn2,N+1);subplot(2,1,2);plot(f,abs(Rf2);axis(-30 30 0 max(abs(Rf2);title('帶通理想白噪聲頻域');xlabel('頻率(Hz)');grid on;figure(3);%3、任意輸入一個信號如，繪制其時頻域波形。st=cos(2*pi*15*t1);subplot(2,1,1);plot(t1,st);%grid on;xlabel('t(s)');title('s

4、(t)的時域');axis(-1 1 -1.1 1.1);sf=fftshift(fft(st,N+1);sfa=abs(sf);subplot(2,1,2);plot(f,sfa);grid on;xlabel('w(Hz)');title('s(t)的頻域')axis(-30 30 min(sfa) max(sfa);a=1;b1=0no=a*randn(1,length(t1)+b1;rt=st+no;%4、產(chǎn)生一個正態(tài)分布的高斯白噪聲，繪制的時頻域波形。figure(4);subplot(2,1,1);plot(t1,rt);xlabel(

5、9;t(s)');title('r(t)的時域');%加噪axis(-1 1 min(rt) max(rt);rf0=fft(rt,N+1)rf=fftshift(fft(rt,N+1);rfa=abs(rf);subplot(2,1,2);plot(f,rfa);grid on;xlabel('w(Hz)');title('r(t)的頻域')axis(-30 30 min(rfa) max(rfa);%5、將通過帶通濾波器，繪制通過帶通濾波器后的時頻域波形。fz=15;B=1;pf=zeros(1,length(f);pf(length

6、(f)/2+round(fz-B)/df):length(f)/2+round(fz+B)/df)=1;pf(length(f)/2-round(fz+B)/df):length(f)/2-round(fz-B)/df)=1;Y=rf0.*fftshift(pf);y=ifft(Y);figure(5)subplot(2,1,1);plot(t1,real(y);title('去噪信號');xlabel('時間(s)'); axis(0 2 min(real(y) max(real(y); grid on; subplot(2,1,2);plot(f,abs(f

7、ftshift(Y);title('去噪頻譜');xlabel('頻率(Hz)');axis(-30 30 min(abs(Y) max(abs(Y); grid on;%6、計算高斯白噪聲信道的信噪比為-10dB到10dB，每2dB遞進(jìn)，發(fā)送端信號為0、1時%的誤碼率情況。figure(6);m=1;n=10000; a=rand(m,n); %產(chǎn)生n×m的均勻分布的隨機(jī)數(shù)for i=1:length(a) if a(i)>0.5 b(i)=1; %產(chǎn)生m×m的雙極性信號 else b(i)=0; end;end;S=0.5; %信號功率snr=-10:2:10; %信噪比(dB)SNR=exp(snr*log(10)/10);for l=1:length(snr) N=S/SNR(l); %噪聲功率 rk=b+sqrt(N)*randn(m,n); %加上高斯白噪聲 for k=1:length(rk) if rk(k)>0.5; dec(k)=1; else dec(k)=0; end; end; %dec=sign(rk); %判決，大于0判為1，小于0判為-1 err(l)=sum(abs(dec-b)/2)/length(b); %計算