友友们，我按照这个编出来不正确

[zt326542506]

 

clc; clear;

% 计算 H(z) 的部分分式展式
B = [1]; % 分子
A = [1 -0.75 0.125]; % 分母

[r, p, k] = residue(B, A); % 部分分式展式

% 显示结果
disp('部分分式展式结果：');
disp(['r: ', num2str(r')]);
disp(['p: ', num2str(p')]);
disp(['k: ', num2str(k')]);

% 通过 Z 变换求单位脉冲响应
n = 0.5; % 假设 n = 0.5

F1 = sym('2/(1-1/2*z^{-1})'); % 第一个部分
F2 = sym('1/(1-1/4*z^{-1})'); % 第二个部分

% 反转 Z 变换
f1 = iztrans(F1);
f2 = iztrans(F2);

% 计算单位脉冲响应 h
h = f1 + f2;

% 使用 impz 计算系统的单位脉冲响应
hn = impz(B, A, n); % 计算单位脉冲响应

% 绘制单位冲击响应图
figure;
stem(hn);
xlabel('n');
ylabel('h[n]');
title('单位脉冲响应');
grid on;