x=randint（151）?
Fc=7000;?%fft要符合奈奎斯特低通采樣定理
Fd=128;?
Fs=65536;?
[yT]=?dmod（xFcFdFs‘ASK‘4）;?%dmod數字通帶調制amod模擬調制ask幅度鍵控調制[yt]=amod（xFcFs‘amdsb-tc‘）;%
len=length（y）;?%y=7680
?
%%pause???%?Press?any?key?to?see?plot??after?sampling?
subplot（211）;?
plot（y）??
xlabel（‘points‘）;?
title（‘sampling?time?domain?signal‘）;?
subplot（212）;?
fft_y=fft（y4069）;?%?

am_y1=abs（fft_y）;??
f?=?Fs*（0:2048）/4096;?%?
plot（fam_y1（1:2049））?
title（‘sampling?frequency?domain?signal‘）;?
xlabel（‘Frequency‘）;?
figure（7）;
plot（abs（fft_y））
?
len=length（y）;?
std_val=0.2;?
n=std_val*randn（len1）;?
y=y+n;?
z=y;?
?
%%pause?%see?plot?after?adding?nosie
figure（2）;
subplot（211）;?
plot（z）?
xlabel（‘points‘）;?
title（‘adding?noise?time?domain?signal‘）;?
subplot（212）;?
fft_z=fft（z4096）;?
am_z=abs（fft_z）;?
f?=?Fs*（0:2048）/4096;?
plot（fam_z（1:2049））?
title（‘adding?noise?frequency?domain?signal‘）;?
xlabel（‘Frequency‘）;?
?
int_z=zeros（len1）;?
for?i=1:len?
????int_z（i）=int16（y（i）*32768）;?%沒有添加噪聲直接變換
end;??
figure（8）;
plot（int_z）?
nco_freq=Fc-Fs;?
ad_freq=Fs;?
len=length（y）;?
nco=zeros（len1）;?
nco_out=zeros（len1）;?
for?i=1:len?
????nco（i）=cos（nco_freq*2*pi*i/ad_freq）;?%軟件無線電書本p75-p76頁
????nco_out（i）=nco（i）*int_z（i）;?
end;?
?
%%pause?%see?plot?after?mixing?
figure（3）;
subplot（211）;?
plot（nco_out）?
xlabel（‘points‘）;?
title（‘mixing?time?domain?signal‘）;?
subplot（212）;?
fft_nco=fft（nco_out4096）;?
am_nco=abs（fft_nco）;?
f?=?Fs*（0:2048）/4096;?
plot（fam_nco（1:2049））?
title（‘mixing?frequency?domain?signal‘）;?
xlabel（‘Frequency‘）;?%ok
?
m1=16;?
D1=4;?
cic_mid1=zeros（len1）;?%len=7680
cic_out1=zeros（len/D11）;?
cic_mid2=zeros（len/D1*D11）;?
cic_mid3=zeros（len/D1*D1*D11）;?
cic_mid4=zeros（len/D1*D1*D1*D11）;?
cic_mid5=zeros（len/D1*D1*D1*D1*D11）;?
?
for?i=m1:len?
????for?j=1:m1?
????????cic_mid1（i）=cic_mid1（i）+nco_out（i-j+1）;?
????????cic_mid2（i）=cic_mid2（i）+cic_mid1（i-j+1）;?
????????cic_mid3（i）=cic_mid3（i）+cic_mid2（i-j+1）;?
????????cic_mid4（i）=cic_mid4（i）+cic_mid3（i-j+1）;?
????????cic_mid5（i）=cic_mid5（i）+cic_mid4（i-j+1）;?
??