主要内容

本页采用了机器翻译。点击此处可查看英文原文。

simulateFrequencyResponse

在频域中对电路进行仿真

自 R2025a 起

说明

[realImpedance,imaginaryImpedance] = simulateFrequencyResponse(circuitModel,frequency) 在频域中对等效电路模型 circuitModel 进行仿真,并返回实部和虚部阻抗。

示例

示例

全部折叠

本例演示了如何使用 simulateFrequencyResponse 函数在频域中仿真两种不同的电化学阻抗谱 (EIS) 模型。

建立一个分数阶等效电路模型,其中一个电阻与另一个电阻和一个电容的并联电路串联。

eisfom = eisModel("R0+(R1,C1)")
eisfom = 
  EISModel with properties:

     CircuitTopology: "R0+(R1,C1)"
       NumParameters: 3
       ParameterList: ["R0"    "R1"    "C1"]
     ParameterValues: [0.0250 0.0035 150]
    CircuitImpedance: "(R0 + ((R1 * (1/(i*w*C1))) / (R1 + (1/(i*w*C1)))))"
       ParameterList: ["R0"    "R1"    "C1"]

请指定要对模型进行仿真的频率点。

freq = logspace(-2,5,100)

在这些频率下对 EIS 电路进行仿真,并绘制频率响应曲线。

[simRealZ, simImagZ] = simulateFrequencyResponse(eisfom,freq);
plot(simRealZ,-simImagZ,linewidth=1.5)
axis equal

Nyquist plot of the frequency response for the R0+(R1,C1) circuit

现在在 EIS 电路拓扑中加入一个瓦尔堡阻抗项,并观察频率响应如何变化。

eisfom.CircuitTopology = "R0 + (R1,C1) + SIW1";
eisfom = 

  EISModel with properties:

     CircuitTopology: "R0 + (R1,C1) + SIW1"
       NumParameters: 4
       ParameterList: ["R0"    "R1"    "C1"    "SIW1"]
     ParameterValues: [0.0250 0.0035 150 1.0000e-03]
    CircuitImpedance: "((R0 + ((R1 * (1/(i*w*C1))) / (R1 + (1/(i*w*C1))))) + (SIW1*(w)^(-1/2) - SIW1*i*(w)^(-1/2)))"
       ParameterList: ["R0"    "R1"    "C1"    "SIW1"]

对该电路进行仿真,并绘制频率响应曲线。

[simRealZ, simImagZ] = simulateFrequencyResponse(eisfom,freq);
plot(simRealZ,-simImagZ,linewidth=1.5)
axis equal

Nyquist plot of the frequency response for the R0+(R1,C1) + SIW1 circuit

输入参数

全部折叠

用于在频域中进行仿真的等效电路模型,指定为 ECMEISModel 对象。

频率,指定为标量或向量。

数据类型: double

输出参量

全部折叠

电路阻抗的实部,以标量或向量形式返回。

数据类型: double

电路阻抗的虚部,以标量或向量形式返回。

数据类型: double

版本历史记录

在 R2025a 中推出