clear?all;?
x=randn（101）>0;?
fl=400;fh=800;fs=10000;?
ts=1/200;??%200
tt=（0:1/fs:ts）;?
t=[tt;tt+ts;tt+2*ts;tt+3*ts;tt+4*ts;tt+5*ts;tt+6*ts;tt+7*ts;tt+8*ts;tt+9*ts];?
m=length（tt）;
y=zeros（10m）;?
i=1;?
%FSK?modulate
while?i<=10?
y（i:）=x（i）*sin（2*pi*fh*t（i:））+~x（i）*sin（2*pi*fl*t（i:））;?
i=i+1;?
end?
t=reshape（t‘length（tt）*101）;?
y=reshape（y‘length（tt）*101）;?
subplot（221）;plot（ty）;?
title（‘信號時域‘）;?

n=length（