(c) 根据(a)中Ud，采用吸收功率为1kW的直流电流源Id。 解: Matlab程序修改如下: 第二行
clc,clg,clear 将clg改为clf 第七行
改为A=[-rs/ls(空格) -1/ls;1/cd (空格) -1/(cd*rload)]; 书中的空格不明显； 第十五行 Alfa0 删掉 第二十行
X=[il(1) vc(1)]; 改为: X=[il(1) ;vc(1)]
Pspice 输入文件:
SUBCKT DIODE_WITH_SNBU 101 102 改为:
SUBCKT DIODE_WITH_SNUB 101 102 运行结果如图:
4-13．图4.6(b)所示的单相整流电路中，i = I，计算THD，DPF，PF和CF。
d
d
图4.6 Ls = 0时二极管桥式整流电路的理想结构
解：由：id?Id，可知：
is为180?标准方波，波形图如下所示：
将is傅里叶展开：
is?4Id11(sin?1t?sin3?1t?sin5?1t?...) ?35Is1?4Id 2?因为Is是180?的交流方波，is的有效值Is=Id
IdisIs2?Is12%THDi?100%??100%?Is1Is14Id2总谐波畸变率： Id2?()2??100%??66.1%4Id2?假设负载为纯电阻负载，所以：
?co?s1?1 位移功率因数 DPF4IdIs1?2??1?0.9 则负载功率因数 PF?DPFIsId4Id波形因数（通常为峰值与有效值之比） CF?Is.peak???1.27 IsId4-14. 利用Matlab程序求解题4.12，计算出THD，DPF，PF和CF。
题图 桥式二极管整流电路
解：将整流电路图在simulink中搭建仿真模块，仿真模块如下图所示：
按照题图设置各元件的仿真参数，Us=120V，60Hz，Ls=1mH，Rs=0.4???PL=1kW，Start time为0，stop time为0.05，算法选取“ode23tb”。载入powergui模块，离散时间设置为“1e-5”，在FFT analysis窗口中，可以看到基波分量和各次谐波分量的含量，其中DC component为平均电压Ud 的值。
在“structure”菜单下可以改变要观察的各测量量；在“Display style”中分别选择“list”和“bar”，分别显示各次谐波清单列表和柱状图。显示结果如下：（见下页）
(1)输出电流Id傅里叶分析如下图：
从图中可知：
（THD）给出了电流Id的畸变率：THD=104.77%。 （2）输出电压Ud的傅里叶分析如下：
从图中可知：
（THD）给出了电流Ud的畸变率：THD=112.28%。
说明：因为输出电流电压直流成分极大，所以谐波含量极高。 （3）交流侧电流Is傅里叶分析如下：
从图中可知：
（THD）给出了电流Ud的畸变率：THD=42.95%。
说明：因为输入电流电压漏感影响不大，所以谐波含量较低。 各次谐波列表如下：
由FFT分析谐波列表可知，电流的基波分量相位θi=-27.1°、θv=0°。故其相位差为
Φ=-27.1?（滞后），所以DPF=cosΦ=0.89。 傅里叶分析可知电流的基波分量Is1=120.6A 由谐波畸变率公式
%TDH?100?可求得：Is=131.25A，故
PF?IS?IS1IS122
Is1120.6DPF??0.89?0.818 Is131.254-15.图4.20所示的单相整流电路中，Us = 120V，频率50Hz，Ls = 2mH，Rs =
0.4?，负载的瞬时功率pd(t) = 1kW。利用Pspice软件，做出Cd分别为：200、500、
1000和1500?F时，THD、DPF、PF以及换相压降△Ud(p-p)的函数曲线，并分析直流侧滤波电容的作用。
解：当Cd分别为：200、500、1000和1500?F时，电源侧电流THD变化情况分别如下图所示：