R = 1000;% 电阻的阻值为1000Ω C = 1e-6;% 电容的容值为1μF R1 = 2000;% 并联电阻1的阻值为2000Ω R2 = 4000;% 并联电阻2的阻值为4000Ω f = logspace（0， 7， 10000）;% 生成10000个对数分布的频率点 % 计算幅频特性和相频特性 H1 = 1./（1 + 1j f/fH）; % 高通滤波器的幅频特性 H2 = 1./（1 + 1j f/fL）;% 低通滤波器的幅频特性 H3 = 1./（1 + 1j f/f）; % 并联后滤波器的幅频特性 phi1 = atan（-f/fH）; % 高通滤波器的相频特性 phi2 = atan（-f/fL）; % 低通滤波器的相频特性 phi3 = atan（-f/f） - atan（f./（f^2 + fH fL））;% 并联后滤波器的相频特性 % 绘制幅频特性曲线 loglog（f， abs（H1）， 'r'， f， abs（H2）， 'b'， f， abs（H3）， 'g'）;legend（'高通滤波器'， '低通滤波器'， '并联后滤波器'）;xlabel（'频率 （Hz）'）;ylabel（'幅值'）;% 绘制相频特性曲线 semilogx（f， phi1， 'r'， f， phi2， 'b'， f， phi3， 'g'）;legend（'高通滤波器'， '低通滤波器'， '并联后滤波器'）;xlabel（'频率 （Hz）'）;ylabel（'相位 （rad）'）;没有定义fH和fL

非常抱歉，之前的代码中确实存在一些错误。以下是修改后的代码：

R = 1000; % 电阻的阻值为1000Ω
C = 1e-6; % 电容的容值为1μF
R1 = 2000; % 并联电阻1的阻值为2000Ω
R2 = 4000; % 并联电阻2的阻值为4000Ω
f = logspace(0, 7, 10000); % 生成10000个对数分布的频率点

% 计算截止频率
fH = 1/(2*pi*R*C); % 高通滤波器的截止频率
fL = 1/(2*pi*R*C); % 低通滤波器的截止频率
f = 1/(2*pi*C*(1/R1+1/R2)); % 并联后滤波器的截止频率

% 计算幅频特性和相频特性
H1 = 1./(1 + 1j*f./fH); % 高通滤波器的幅频特性
H2 = 1./(1 + 1j*f./fL); % 低通滤波器的幅频特性
H3 = 1./(1 + 1j*f./f); % 并联后滤波器的幅频特性
phi1 = atan(-f./fH); % 高通滤波器的相频特性
phi2 = atan(-f./fL); % 低通滤波器的相频特性
phi3 = atan(-f./f) - atan(f./(f.^2 + fH.*fL)); % 并联后滤波器的相频特性

% 绘制幅频特性曲线
loglog(f, abs(H1), 'r', f, abs(H2), 'b', f, abs(H3), 'g');
legend('高通滤波器', '低通滤波器', '并联后滤波器');
xlabel('频率 (Hz)');
ylabel('幅值');

% 绘制相频特性曲线
semilogx(f, phi1, 'r', f, phi2, 'b', f, phi3, 'g');
legend('高通滤波器', '低通滤波器', '并联后滤波器');
xlabel('频率 (Hz)');
ylabel('相位 (rad)');

这样运行代码后，就可以得到高通、低通和并联后滤波器的幅频特性曲线和相频特性曲线，并且计算出截止频率了。