用matlab编写程序，仿真高斯色噪声下正弦信号频率估计问题，高斯色噪声由AR1模型产生，模型参数为0.8，用蒙特卡洛方法，给出估计频率随信噪比变化的方差

% 高斯色噪声下正弦信号频率估计问题仿真
clear all;
clc;

% 信号参数
f = 0.1; % 正弦信号频率
fs = 100; % 采样频率
N = 1000; % 采样点数
t = (0:N-1)/fs; % 时间序列

% 生成正弦信号
x = sin(2*pi*f*t);

% 噪声参数
ar1 = 0.8; % AR1模型参数
sigma = 1; % 噪声标准差

% 蒙特卡洛模拟次数
M = 100;

% 信噪比范围
snr_range = -20:5:20;

% 初始化方差数组
variance = zeros(size(snr_range));

for k = 1:length(snr_range)
% 当前信噪比
snr = snr_range(k);

% 噪声方差
noise_variance = (10^(-snr/10))*(std(x))^2;

% 计算噪声
noise = zeros(N,1);
noise(1) = randn;
for n = 2:N
noise(n) = ar1*noise(n-1) + sqrt(1-ar1^2)*randn;
end
noise = noise*sqrt(noise_variance);

% 添加噪声
y = x + noise;

% 估计频率
f_hat = zeros(M,1);
for m = 1:M
% 加窗
w = hamming(N).*y;

% 快速傅里叶变换
Y = fft(w);

% 峰值索引
[~,idx] = max(abs(Y));

% 估计频率
f_hat(m) = (idx-1)*fs/N;
end

% 计算方差
variance(k) = var(f_hat);
end

% 画图
plot(snr_range,variance);
xlabel('信噪比/dB');
ylabel('频率估计方差');
title('高斯色噪声下正弦信号频率估计方差仿真');