function?y=ask_2（snr_in）;
snr=10^（snr_in/10）;
x=rand（1001）;%產生隨機序列
h=length（x）;%計算序列長度
fc=10000;fs=100000;%設定采樣頻率以及高頻和低頻（調制頻率）?
ts=10/fc;%設定碼元寬度;每碼元100個點?
t=0:1/fs:ts*h-1/fs;%產生全部碼元總的時間序列起始時刻為0步長為采樣頻率的倒數即采樣時間結束時刻為序列結束時刻
l=length（t）/h;%計算出這個序列的長度
t=reshape（t‘lh）;%將序列處理成矩陣形式，每行元素為一個碼元內的所有采樣點l，行數為碼元數n??：reshape（lh）
t=t‘;%因矩陣是按列存儲，故須作轉置處理
y=zeros（hl）;
for?i=1:h
????if?x（i）>=0.5
?????????y（i:）=3.5*sin（2*pi*fc.*t（i:））;
????else?y（i:）=sin（2*pi*fc.*t（i:））;%賦值
????end
end
y0=reshape（y‘l*h1）;
power=y0‘*y0;
powerzao=power/snr;%以下是產生高斯白噪聲序列
a=sqrt（powerzao/（l*h））;
zao=a*randn（l*h1）;
powerzao2=zao‘*zao;
y=y0+zao;%加噪以后的信號
t=reshape（t‘l*h1）;
subplot（311）;plot（ty）;
grid?on
title（‘2ASK信號的時域圖形‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????904??2008-12-20?16:02??ask_2.m

?????文件????????862??2008-11-05?16:13??ask_4.m

?????文件????????957??2008-12-20?16:03??fsk_2.m

?????文件????????937??2008-12-20?16:03??fsk_4.m

?????文件????????607??2008-12-22?09:56??main.m

?????文件????????683??2008-11-13?14:00??psk_2.m

?????文件???????2154??2008-11-03?16:57??QAM_16.m

?????文件????????867??2008-12-03?15:53??qpsk.m

?????文件????????779??2008-12-20?15:54??TX.m

?????文件????????865??2008-12-04?10:05??TZSB.m

-----------?---------??----------?-----??----

?????????????????9615????????????????????10