乘法鉴相器的matlab仿真,Matlab仿真costas环代码问题。

clc;

clear all;

fs=32e6;    %采样频率

fc=6e6;     %载波频率

df=2000;    %初始频偏

fi=fc+df;   %输入信号载波频率

num=0.5e6;  %数据长度

SNR=30;     %输入数据信噪比(dB)

ts=1/fs;

%本地VCO的增益，也相当于最大捕获频偏

K0=10e3;

%环路滤波器系数

c1=2^(-5);

c2=2^(-11);

%环路输入信号为单载波信号，信噪比为SNR(dB)

din=sin(2*pi*fi*(0:num-1)*ts);

datain=awgn(din,SNR);

%设计IIR低通滤波器

N=3;    %滤波器阶数

R=60;   %阻带衰减(dB)

Wn=9e6; %截止频率

[lpf_b,lpf_a]=cheby2(N,R,Wn*2/fs);

%估计滤波器的增益(求取零点系数绝对值之和)

glpf=sum(abs(lpf_b))

%本地载波输出信号初值

o_sin=sin(2*pi*fc*(0:num-1)*ts);

o_cos=cos(2*pi*fc*(0:num-1)*ts);

%乘法运算输出信号初值

mult_i=zeros(1,num);

mult_q=zeros(1,num);

%鉴相器输出信号初值

pd=zeros(1,num