MATLAB电路方程组求解代码无法打印输出问题求助

问题原因 #

1. 方程组逻辑矛盾，无解
   第二个方程化简后存在矛盾：
   原第二个方程：

   Q == 1/(sqrt((C2*R2)/(C1*R1))*sqrt((C1*R1)/(C2*R2))*sqrt((C1*R2)/(C2*R1)))
   

   其中sqrt((C2*R2)/(C1*R1))*sqrt((C1*R1)/(C2*R2))计算结果为1，方程简化为Q == 1/sqrt((C1*R2)/(C2*R1))。结合第三个方程Tp == Q*(sqrt((C1*R2)/(C2*R1)))，代入后会得到Tp == 1，但代码中计算的Tp≈0.297，二者矛盾，导致vpasolve找不到有效解，sR1、sR2、sC2均为空值。

2. 变量赋值顺序错误
   代码提前定义sol = [sC2,sR1,sR2]时，sC2、sR1、sR2还未通过vpasolve赋值，后续即使求解成功，sol依然是空数组，fprintf调用空数组自然无输出。

3. 输出逻辑错误
   即使求解得到结果，fprintf使用的是提前定义的空sol，而非求解得到的sR1、sR2、sC2。

解决方法 #

1. 修正方程组逻辑
   对照电路理论中的Q值公式，替换矛盾的第二个方程（需根据你的实际电路拓扑确认正确公式）。例如常见二阶电路Q值公式可参考Q = 0.5*sqrt((C1*R2^2)/(C2*R1))。

2. 调整变量赋值与输出顺序
   在vpasolve求解完成后，再将结果存入sol，或直接使用求解得到的变量输出。

3. 增加解的有效性判断
   加入判断逻辑，当求解无解时给出明确提示，避免无输出的困惑。

修正后的示例代码 #

f0 = 1200;
Q = 3;
W0 = 2*pi*f0;
TpdB = -10.5;
Tp = 10^(TpdB/20);
C1 = 680e-9;
syms R1 R2 C2 positive;

% 替换为符合电路逻辑的方程组（此处Q的公式为示例，需根据实际电路调整）
[sR1,sR2,sC2] = vpasolve(...
    W0 == 1/sqrt(C1*C2*R1*R2),...
    Q == 0.5*sqrt((C1*R2^2)/(C2*R1)),...
    Tp == Q*(sqrt((C1*R2)/(C2*R1))),...
    [R1,R2,C2]);

% 判断解的有效性并输出
if isempty(sR1) || isempty(sR2) || isempty(sC2)
    disp('未找到有效解，请检查方程组是否正确');
else
    sol = [sR1, sR2, sC2];
    fprintf('R1: %e Ω, R2: %e Ω, C2: %e F\n', sol);
end

额外建议 #

• 新手可先使用solve尝试符号解，若符号解不存在，再用vpasolve数值求解，更易排查方程问题。
• 求解前单独化简每个方程，验证逻辑是否自洽，避免矛盾方程导致无解。

内容的提问来源于stack exchange，提问作者Old Jabba Jabba