干擾信號(hào)MATLAB仿真

1、干擾信號(hào)干擾信號(hào)MATLABMATLAB仿真仿真實(shí)驗(yàn)?zāi)康膶?shí)驗(yàn)?zāi)康姆抡娈a(chǎn)生4種以上的通信干擾信號(hào)。畫出信號(hào)的時(shí)域、頻域波形。實(shí)驗(yàn)思路與步驟實(shí)驗(yàn)思路與步驟 產(chǎn)生一個(gè)高斯白噪聲， 利用Matlab自帶的fir1函數(shù)產(chǎn)生一個(gè)低通濾波器，限制高斯白噪聲的帶寬，由此產(chǎn)生了視頻噪聲。 利用產(chǎn)生的視頻噪聲，分別代入噪聲調(diào)幅干擾的時(shí)域表達(dá)式，并且進(jìn)行100次的積累后求平均值，對(duì)其進(jìn)行快速傅里葉變換后，由此畫出噪聲調(diào)幅干擾頻域波形 重復(fù)上述步驟，分別代入噪聲調(diào)頻干擾和噪聲調(diào)相干擾的時(shí)域表達(dá)式，分別畫出其時(shí)域波形和頻域波形。實(shí)驗(yàn)結(jié)果實(shí)驗(yàn)結(jié)果 視頻噪聲噪聲調(diào)幅干擾 噪聲調(diào)頻干擾 噪聲調(diào)相干擾程序代碼 clear al

2、l;close all;clc;%清除變量 T=50e-6;%采樣時(shí)間 fs=150e6;%采樣頻率 N=T*fs;%采樣點(diǎn)數(shù) detlf=10e6;%濾波器截止頻率 f1=50e6;%調(diào)制信號(hào)中心頻率 m=0.5;%調(diào)制度 kfm=5e6;%調(diào)頻斜率 kpm=5;%調(diào)相斜率 M=100;%積累次數(shù) p=fft(fir1(N-1,detlf/fs*2);%濾波器頻譜s=0;for i=1:100 xn=ifft(fft(random(Normal,0,1,1,N).*p);%高斯白噪聲通過濾波器j=abs(fft(xn);s=s+j;endfigure(1)t=0:1/fs:T-1/fs;pl

3、ot(t*1e6,xn);xlabel(us); title(視頻噪聲時(shí)域波形); figure(2) f=(0:N-1)*fs/N; y1=fft(xn,N); mag1=abs(y1); plot(f*1e-6,y1);%視頻噪聲頻域波形 axis(-2 15 -150 150); xlabel(MHZ); title(視頻噪聲頻域波形);n=1:N;zn=(1+m*cos(2*pi*xn).*cos(2*pi*f1/fs*n);%噪聲調(diào)幅干擾表達(dá)式figure(3)plot(t*1e6,zn);title(噪聲調(diào)幅干擾時(shí)域波形);xlabel(us); figure(4)f=(0:N-1

4、)*fs/N;y2=fft(zn,N);mag2=abs(y2);plot(f*1e-6,mag2);xlabel(MHZ);title(噪聲調(diào)幅干擾頻域波形); sum(1)=0; for i=1:N-1; sum(i+1)=xn(i)+sum(i); end xn=sum/fs; wn=cos(2*pi*f1*t+2*pi*kfm*xn);%噪聲調(diào)頻干擾表達(dá)式 figure(5) plot(t*1e6,wn); title(噪聲調(diào)頻干擾時(shí)域波形); xlabel(us); figure(6) f=(0:N-1)*fs/N; y3=fft(wn,N); mag3=abs(y3); plot(f*1e-6,mag3); xlabel(MHZ); title(噪聲調(diào)頻干擾頻域波形); sum(1)=0; for i=1:N-1; sum(i+1)=xn(i)+sum(i); end xn=sum/fs; on=cos(2*pi*f1*t+kpm*xn);%噪聲調(diào)相干擾表達(dá)式 figure(7) plot(t*1e6,on); title(噪聲調(diào)相干擾時(shí)域波形); xlabel(us); figure(8) f=(0:N-1)*fs/N; y4=