帮我用MATLAB写一个prony算法的算例

以下是一个简单的Prony算法的MATLAB代码示例：

% 生成信号
n = 100; % 信号长度
t = linspace(0,1,n); % 时间向量
f = [2, 3, 5]; % 频率向量
a = [1, 0.5, 0.2]; % 幅值向量
x = sum(a.'*sin(2*pi*f.'*t), 1); % 信号

% Prony算法
p = 3; % 信号中包含的正弦波的数量
M = n-p+1; % Toeplitz矩阵的大小
X = toeplitz(x(p:n), x(p:-1:1)); % Toeplitz矩阵
Y = -x(p+1:n).'; % Y向量
a_hat = Y/X; % Prony估计的系数向量

% 绘制结果
t_est = linspace(0,1,M); % 估计的信号长度
x_est = sum(a_hat.'*sin(2*pi*f.'*t_est), 1); % 估计的信号
figure;
subplot(2,1,1);
plot(t,x);
xlabel('时间');
ylabel('幅值');
title('原始信号');
subplot(2,1,2);
plot(t_est,x_est);
xlabel('时间');
ylabel('幅值');
title('Prony算法估计的信号');

% 显示结果
disp(['原始信号的频率：', num2str(f)]);
disp(['Prony算法估计的频率：', num2str(sort(flipud(abs(angle(roots([1 -a_hat.'].'))))*(n-1)/(2*pi)), '%.3f ')]);

注释：此算例生成一个长度为100的信号，包含3个正弦波成分，使用Prony算法估计正弦波的频率和幅值。然后，该算法绘制原始信号和估计信号的图形，并将原始频率和估计频率列出来。