%噪聲調幅干擾?
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