fskjietiao.m

来自「这是一个FSK ASK 升余弦的的调制程序」· M 代码 · 共 33 行

function fskjietiao
M=2;                                   %进行M元调制解调,M定为2
Fc=2;                                  %载波频率设定为10Hz
Fd=1;                                   %解调后数字信号速率为1Hz
Fs=4;                                  %模拟信号采样速率设定为50Hz
x=randint(10,1,M);                     %产生要调制的数字信号
[y,t]=dmod(x,Fc,Fd,Fs,'fsk',M);         %用M元FSK方式进行调制
subplot(3,1,1);plot(t,y)                %FSK调制图象
ynoisy=y+0.1*randn(Fs/Fd*10,1) ;        %加入高斯噪声
subplot(3,1,2);plot(t,ynoisy);           
z=ddemod(ynoisy,Fc,Fd,Fs,'fsk',M); %解调接受信号
f=2;
t=0:2*pi/99:2*pi;
cp=[];mod=[];bit=[];
for n=1:length(z)
    if z(n)==0;
        cp1=zeros(1,100);
        bit1=zeros(1,100);
    else z(n)==1;
        cp1=ones(1,100);
        bit1=ones(1,100);
    end
    c=sin(f*t);
    cp=[cp cp1];
    mod=[mod c];
    bit=[bit bit1];
end
ask=cp.*mod;
subplot(3,1,3);
plot(bit,'lineWidth',1.5);grid on;
ylabel('Binary Signal');
axis([0 100*length(z) -2.5 2.5]);
s=symerr(x,z)