1.产生信号.m

来自「自适应线性预测编码问题」· M 代码 · 共 20 行

clear;
clc;
syms w;
1/int(1/(2*pi)/(abs(1+0.1*exp(-i*w)-0.8*exp(-i*w*2))^2),w,-pi,pi)
                                    % 计算噪声方差 结果为0.2700
V=sqrt(0.27)*randn(500,1);          % 产生白噪声序列V(n)
var(V)                              % 检验方差

X=[0+V(1);1+V(2)];
for n=3:500                         % 产生准信号序列
    x=-0.1*X(n-1)+0.8*X(n-2)+V(n);
    X=[X;x];
end
X=X(201:500);                       % 截取信号序列X(n)
var(X)                              % 检验方差

figure(1)
plot(V)                             % 绘制噪声波形
figure(2)
plot(X)                             % 绘制信号波形