jiema.m

来自「PSK编码解码的基于MATLAB的完整程序」· M 代码 · 共 23 行

function y=jiema()
load('modulation_fsk')
n=length(modulation_fsk)/16;
y(1)=0.5;
f1=4000;
f0=8000;
for i=1:n
    for k=1:16;
        m1(i,k)=cos(2*pi*f1*k/16000);
        m2(i,k)=cos(2*pi*f0*k/16000);
        fy1(k)=m1(i,k)*modulation_fsk((i-1)*16+k);
        fy2(k)=m2(i,k)*modulation_fsk((i-1)*16+k);
    end
        y1=sum(fy1);
        y2=sum(fy2);
    if y1>y2
        y(i+1)=y(i)+0.02;     
    else
        y(i+1)=y(i)-0.02;
    end
end
plot(y)
wavwrite(y,'你所要另存为的文件名')