clear?all
global?dt?t?f?df?N?T?%全局變量
close?all?
N=2^20;?
dt=0.001;?%ms?
df=1/（N*dt）;?%KHz?
T=N*dt;?%截短時間?
Bs=N*df/2;?%系統(tǒng)帶寬?
t=linspace（-T/2T/2N）;?%時域橫坐標?
f=linspace（-BsBsN）+eps;?%頻域橫坐標?
fm=1;?%模擬基帶信號頻率為1kHz?
fc=6;?%正弦載波信號頻率為6kHz?
mt=cos（2*pi*fm*t）;?%模擬基帶信號?
m（mt>0）=1;??%發(fā)送比特
m（mt<0）=0;
c=cos（2*pi*fc*t）;?%正弦載波信號?
s=m.*c;?%DSB-SC?AM已調(diào)信號?
S=T2F（s）;?
n=awgn（s20）;?%AWGN信道噪聲
r=s+n;?%接受信號?
y=r.*c;?
Lpfil=LPF（fc）;?%低通濾波器?
Y=T2F（y）;?
yrt=real（F2T（Lpfil.*Y））;?
yr=yrt-mean（yrt）;
yr（yr>0）=1;??%發(fā)送比特
yr（yr<0）=0;
subplot（311）?
plot（tm‘LineWidth‘1.2）
title（‘輸入信號‘）;
axis（[-1+1-0.21.2*max（m）]）?
xlabel（‘t?（ms）‘）
ylabel（‘s（t）?（V）‘）?
subplot（312）?
plot（ts‘LineWidth‘1.2）?
title（‘2ASK調(diào)制后波形‘）;
axis（[-1+11.2*min（c）1.2*max（c）]）?
xlabel（‘t?（ms）‘）
ylabel（‘s（t）?（V）‘）?
subplot（313）
plot（tyr‘LineWidth‘1.2）?
axis（[-1+1-0.21.2*max（yr）]）
title（[‘2ASK解調(diào)后波形‘]）?
xlabel（‘t?（ms）‘）?
ylabel（‘s（t）?（V）‘）

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

?????文件???????1040??2007-10-08?19:25??ASKmod.m

?????文件????????151??2007-10-08?19:26??F2T.m

?????文件????????163??2007-10-08?19:26??LPF.m

?????文件????????829??2007-10-08?19:29??MSKmod.m

?????文件????????783??2007-10-08?19:36??OFDMmod.m

?????文件????????248??2007-10-08?19:31??QAMmod.m

?????文件????????199??2007-10-08?19:27??QPSKmod.m

?????文件????????154??2007-10-08?19:26??T2F.m

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

?????????????????3567????????????????????8