bianma.m

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

function y=bianma()
[y,fs,bit]=wavread('你的文件名')
a=y;
n=length(a);
b(1)=a(1);
for i=1:n-1,
    if b(i)<a(i+1)
     s(i)=1;
  else s(i)=0;
end
    if s(i)==1;
      b(i+1)=b(i)+0.02;
  else b(i+1)=b(i)-0.02;
end
end
k=length(s);
c=s+1;
for t=1:k,
  for j=1:16,
      y(16*(t-1)+j)=cos(2*pi*4000/c(t)*j/16000);
  end
end