噪声调幅干扰

雷达通信领域噪声调幅干扰的Matlab小例子

%噪声调幅干扰 
function y=noiseAM（u0Nwpp）;


fj=35e6;
Tr=520e-6;    
fs=4*fj; 
t1=0:1/fs:Tr; 
if nargin==0
 N=length（t1）; wpp=0;u0=1; 
end
u=wgn（1Nwpp）;  %产生高斯白噪声6dBW% 
df1=fs/N;n=0:N/2;f=n*df1;  
%加滤波器 生成带限高斯白噪声
wp=10e6; %调制噪声带宽10M
ws=14e6; 
rp=1; rs=60; 
[n1wn1]=buttord（wp/（fs/2）ws/（fs/2）rprs）; 
[ba]=butter（n1wn1）; 
u1=filter（bau）; %带限
j2=fft（u1）;
%figuresubplot（211）plot（t1u1）title（‘高斯带限噪声‘）;
%subplot（212） plot（f10*log10（abs（j2（n+1）*2/N）））title（‘功率谱‘）;% good!

rand（‘state‘ 0）;
%phi=2*pi*rand（1N）;               %phi2=filter（baphi）; % 0－2pi均匀分布
y=（u0+u1）.*cos（2*pi*fj*t1+2）;         %
%J=fft（y）;                        %  or  J=10*log10（abs（fft（y512）））; fftshift（

%figuresubplot（211） plot（t1y）title（‘噪声调幅干扰时域波形‘）axis（[0  Tr  -5  5]）;;
%subplot（212） plot（f10*log10（abs（J（n+1）*2/N）））title（‘功率谱‘）;%???
%clear all; clc;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        948  2009-03-26 16:19  noiseAM.m

----------- ---------  ---------- -----  ----

                  948                    1