function?[y1?fs?fb?fc1?fc2]?=?fsk（xfs?fb?fc1?fc2flag）;
?
%?----------------------------------------------------%
%?x:?bit?stream;????????%
%?（fs:sampling?fre;?fb:bit?rate;?fc:carrier?fre）??????%
%?y:?已調(diào)信號?????????????????????????????????????????%
%-----------------------------------------------------%
if?nargin<5;?
????fc2=2.5;
end;
if?nargin<4;?
????fc1=2;?
end;
if?nargin<3;
????fb=1;???
end;
if?nargin<2;
????fs=32;
end;
if?nargin<1;?
????x=[0?0?1?1?1?0?0?1?0?1?1?0?0?0?1?0];?
end;
T=length（x）/fb;
m=fs/fb;
dt=1/fs;?t=0:dt:T-dt;
xx=2*reshape（repmat（x‘1m）‘1length（x）*m）-1;%RecPulse（x?m?2）;????????????%?基帶信號
y1=?（xx>=0）.*sin（2*pi*fc1*t）+（xx<0）.*sin（2*pi*fc2*t）;???????????????????????%?FS
?
%?------------------------圖像-----------------------------------?%
if?flag==1;
????figure（3）
????subplot（211）;?plot（t?y1?t?xx*.5?[0?T]?[0?0]?‘b:‘）;
????title（‘Signal‘）;
????
????n=length（y1）;?y=fft（y1）/n;?y=abs（y（1:fix（n/2）））*2;
????I=find（y<1e-04）;?y（I）=1e-04;?%x=20*log10（x）;????
????f1=m/n;?f=0:f1:（length（y）-1）*f1;
????subplot（212）;?plot（f?y?‘r‘）;?grid?on;
????title（‘Spectrum‘）;?xlabel（‘f/fb‘）;?zoom?xon;
????set（gcf‘num‘‘off‘‘name‘...
????????[‘FSK‘?blanks（10）]）;
end;
????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????448205??2013-05-09?09:05??大四上通信系統(tǒng)仿真\通信系統(tǒng)仿真sq.docx
?????目錄???????????0??2013-01-20?09:13??大四上通信系統(tǒng)仿真\
?????文件????????1257??2013-01-09?15:37??大四上通信系統(tǒng)仿真\fsk.m
?????文件????????1783??2013-01-09?15:37??大四上通信系統(tǒng)仿真\fskdet.m
?????文件????????2740??2013-01-10?17:16??大四上通信系統(tǒng)仿真\Utiaoz.m