二元序列10110010，画出2FSK波形

%对二元序列10110010，画出2FSK波形，设载波频率w1=2w2=2Rb（码元速率）。
function fskdigital(s,f0,f1)
%本程序实现FSK调制
%s-输入二进制序列，f0，f1-两个载波信号频率
%调用举例：（f0和f1必须是整数）fskdigital([1 0 1 1 0],1,2)
t=0:2*pi:2*pi;%初始定义
cp=[];mod=[];bit=[];
for n=1:length(s);%调制过程
if s(n)==0;
cp1=ones(1,100);
c=sin(f0*t);
bit1=zeros(1,100);
else s(n)==1;
cp1=ones(1,100);
c=sin(f1*t);
bits=ones(1,100);
end
cp=[cp cp1];
bit=[bit bit1];
mod=[mod c];
end
fsk=cp.*mod;
subplot(2,1,1);
plot(bit,'LineWidth',1.5);grad on;%分别画出原信号、已调信号示意
ylabel('Binary Signal');
axis([0 100*length(s) -2.5 2.5]);
subplot(2,1,2);
plot(bit,'LineWidth',1.5);grad on;
ylabel('FSK modulation');
axis([0 100*length(s) -2.5 2.5]);