EF类和E/F类功率放大器（能量转换器）的波形推导和理想仿真--基于Matlab和ADS

EF类和E/F类功率放大器（能量转换器）的波形推导和理想仿真–基于Matlab和ADS

参考论文：Modeling and Analysis of Class EF and Class E/F Inverters With Series-Tuned Resonant Networks(2016)

这篇文章的思路和MTT的文章A Generalized High-Efficiency Broadband Class-E/F3 Power Amplifier Based on Design Space Expanding of Load Network还有CAS的文章Novel Design Space of Broadband High-Efficiency Parallel-Circuit Class-EF Power Amplifiers的底层理论是相似的。

全部工程下载（Matlab+ADS）：EF类和E/F类的波形推导和理想电路仿真工程-基于Matlab和ADS

这个论文我当时手推非常慢，花了一周搞推导吧，搞了个基于Matlab的syms推出的解析表达式和分析流程

众所周期，EF类和E/F类都是混合的模式类型，EF是在E类的基础上控制了偶次谐波使其短路，E/F类是在E类的基础上控制奇次谐波短路。但是，EF类和E/F类都是有其对应的原始电路实现的，下面对Matlab推导和ADS验证过程进行介绍。

目录

  • EF类和E/F类功率放大器（能量转换器）的波形推导和理想仿真--基于Matlab和ADS
    • 1 、EF类和E/F类的简单理论
    • 2 、使用MATLAB计算EF类和E/F类的波形与参数
    • 3 、ADS的理想仿真
    • 4 、EF类和E/F类的基于MATLAB推导过程

1 、EF类和E/F类的简单理论

尽管 E 类PA可以用相对简单的电路拓扑实现100%的理论效率，但它面临着高峰值开关电压的问题。将E类开关模式和谐波调谐技术相结合，可以得到 Class-EF 或 E/F 模式。经典的Class-EF与Class-E/F的电路图与其等效电路如下图所示：
在这里插入图片描述
其中 L 1 L1 L1为射频扼流圈， C 1 C1 C1为并联电容， C 3 {C_3} C3​和部分 L 3 {L_3} L3​谐振在基频（和L谐振在基频）， C 2 L 2 {C_2}{L_2} C2​L2​起到谐波调谐的作用。对于Class-EF类PA，其需要控制偶次谐波，因此n为偶数。类似的，对于Class-E/F类PA，其需要控制奇次谐波，因此n为奇数。

但是，虽然EF类和E/F类控制的谐波次数不同，但都需要满足ZVDS和ZVS条件，这也是推导的关键步骤之一。

但是，在开始推导之前，我们需要关注一些需要预先设置的变量。参考上面的论文，我们需要优先设置要控制的谐波次数q，占空比D，和电容C1与C2的比值 k = C 1 C 2 k=\frac{C_{1}}{C_{2}} k=C2​C1​​。当然，这些参数也会影响波形：

在这里插入图片描述

2 、使用MATLAB计算EF类和E/F类的波形与参数

main.m，运行前设置q,D,k参数！！！
这部分的推导参考 （4 、EF类和E/F类的推导过程），这边直接给出计算代码的效果。首先，我们需要确定的是设置要控制的谐波次数q，占空比D，和电容C1与C2的比值 k = C 1 C 2 k=\frac{C_{1}}{C_{2}} k=C2​C1​​。此处以q=2，D=0.375，k=0.869为例，得到的结果为：
在这里插入图片描述
同时，命令行也会给出对应的设计参数，这些参数对应的是归一化电压和电流，也就是理想的输出功率为1W：
在这里插入图片描述

3 、ADS的理想仿真

使用上面的计算参数，q=2，D=0.375，k=0.869：
在这里插入图片描述
运行仿真，得到的仿真波形如下所示：
在这里插入图片描述
对应的阻抗空间如下所示，其中二次谐波短路正是EF2类的典型特征：
在这里插入图片描述

4 、EF类和E/F类的基于MATLAB推导过程

推导的过程主要是解方程，需要满足下面几个ZVDS和ZVS条件：
在这里插入图片描述
直接解也太麻烦了，所以我们可以使用Matlab的辅助推导来进行求解（classEF_solve.m文件），首先要解出文章中定义变量（A1,B1,A2,B2,p,phi）的表达式，结合上面的条件可以使用solve进行求解，我也是在注释给出了每个等式对应哪个方程：

%% 基本变量计算
clc
clear all
% 已知变量q1 D k，未知变量A1 B1 A2 B2 p phi
syms q1 D k
syms A1 B1 A2 B2 p phi
% Eqn14
q2=(q1*sqrt((k+1)/k));% Eqn5、Eqn12、Eqn16---连续性
eqn1=A1*cos(q1*2*pi*D)+B1*sin(q1*2*pi*D)==A2*cos(q2*2*pi*D)+B2*sin(q2*2*pi*D)-q2^2*p/(q2^2-1)*sin(2*pi*D+phi)+1/(k+1);
% Eqn5、Eqn12、Eqn17---连续性
eqn2=A1==A2*cos(q2*2*pi)+B2*sin(q2*2*pi)-q2^2*p/(q2^2-1)*sin(phi)+1/(k+1);
% Eqn5、Eqn12、Eqn18---导数连续性
eqn3=-A1*q1*sin(q1*2*pi*D)+B1*q1*cos(q1*2*pi*D)==-A2*q2*sin(q2*2*pi*D)+B2*q2*cos(q2*2*pi*D)-q2^2*p/(q2^2-1)*cos(2*pi*D+phi);
% Eqn5、Eqn12、Eqn19---导数连续性
eqn4=B1*q1==-A2*q2*sin(q2*2*pi)+B2*q2*cos(q2*2*pi)-q2^2*p/(q2^2-1)*cos(phi);
% Eqn23
eqn5=2*pi*k*(1-D)/(k+1)+p*(cos(2*pi*D+phi)-cos(phi))*(q2^2/(q2^2-1)-(k+1))+A2/q2*(sin(2*pi*D*q2)-sin(2*pi*q2))+B2/q2*(cos(2*pi*q2)-cos(2*pi*D*q2))==0;
% Eqn24
eqn6=k/(k+1)+p*sin(phi)*(q2^2/(q2^2-1)-(k+1))-(A2*cos(q2*2*pi)+B2*sin(q2*2*pi))==0;[A1,B1,A2,B2,p,phi]=solve(eqn1,eqn2,eqn3,eqn4,eqn5,eqn6,A1,B1,A2,B2,p,phi);simplify(A1)
simplify(B1)
simplify(A2)
simplify(B2)
simplify(p)
simplify(phi)

解出了这些变量就可以继续求目标的波形了，对应的公式也给出了：

%% 波形的计算
clearsyms q1 D k theta
syms A1 B1 A2 B2 p phi
syms tao
% Eqn14
q2=(q1*sqrt((k+1)/k));
assume((D>0)&(D<1))% Eqn12
il2_norm=A2*cos(q2*tao)+B2*sin(q2*tao)-q2^2*p/(q2^2-1)*sin(tao+phi)+1/(k+1);
% Eqn21
ic1_norm=1-p*(k+1)*sin(tao+phi)-il2_norm;
% Eqn26
beta_temp=simplify(int(ic1_norm,tao,2*pi*D,theta));
% Eqn28
vds_norm=simplify(2*pi*beta_temp/(int(beta_temp,theta,2*pi*D,2*pi)));
% Eqn29
is_norm=1-p*(k+1)*sin(theta+phi)-(A1*cos(q1*theta)+B1*sin(q1*theta));

最后是需要计算电路的相关参数，对应的公式也给出在注释，值得注意的是，vx算出来是分段函数，所以最终实现需要用if else语句实现：

%% 电路参数计算
clearsyms q1 D k theta f
syms A1 B1 A2 B2 p phi
syms tao
% Eqn14
q2=(q1*sqrt((k+1)/k));
assume((D>0)&(D<1))
% Eqn12
il2_norm=A2*cos(q2*tao)+B2*sin(q2*tao)-q2^2*p/(q2^2-1)*sin(tao+phi)+1/(k+1);
% Eqn21
ic1_norm=1-p*(k+1)*sin(tao+phi)-il2_norm;
% Eqn26
beta_temp=simplify(int(ic1_norm,tao,2*pi*D,theta));
% Eqn51
beta_integral=int(beta_temp,theta,2*pi*D,2*pi);
% Eqn35
RDC_RL=0.5*((2*pi*(1-D)+A1/q1*sin(2*pi*D*q1)+2*B1/q1*sin(pi*D*q1)^2)/(cos(2*pi*D+phi)-cos(phi)))^2;
% Eqn52
reactance_C1=simplify(2*pi*RDC_RL/beta_integral);
% Eqn58
vx=simplify(2*(int(beta_temp*cos(theta+phi),theta,2*pi*D,2*pi))/(int(beta_temp,theta,2*pi*D,2*pi)));
% Eqn32
im=(2*pi*(1-D)+A1/q1*sin(2*pi*D*q1)+2*B1/q1*sin(pi*D*q1)^2)/(cos(2*pi*D+phi)-cos(phi));
% Eqn51-60
reactance_C2=reactance_C1*k;
reactance_L2=reactance_C2/q1^2;
reactance_Lx=vx/im*RDC_RL;Q_RLC=100;
RL=1/RDC_RL;
C1=1/(2*pi*f*RL*reactance_C1);
C2=1/(2*pi*f*RL*reactance_C2);
L2=reactance_L2*RL/(2*pi*f);
Lx=reactance_Lx*RL/(2*pi*f);
C3=1/(2*pi*f*Q_RLC*RL);
L3=Q_RLC^2*RL^2*C3;

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.hqwc.cn/news/525289.html

如若内容造成侵权/违法违规/事实不符，请联系编程知识网进行投诉反馈email:809451989@qq.com，一经查实，立即删除！

相关文章

【敬伟ps教程】文字处理工具

【敬伟ps教程】文字处理工具

文章目录 文字工具使用方式文字图层文字工具选项字符面板段落面板文字工具使用方式 文字工具（快捷键T），包含横排和直排两种类型 创建文本两种类型：点式文本、段落文本 创建文字方式 1、在画面上单击，出现文字光标，可输入文字，然后需要在工具栏中点击“√”或者 Ctrl+…
阅读更多...
简洁的链式思维（CCoT）提示

简洁的链式思维（CCoT）提示

原文地址&#xff1a;Concise Chain-of-Thought (CCoT) Prompting 传统的CoT导致了输出令牌使用的增加&#xff0c;而CCoT提示是一种旨在减少LLM响应的冗长性和推理时间的提示工程技术。 2024 年 1 月 24 日 Areas where Chain-Of-Thought-like methodology has been introd…
阅读更多...
Docker常用基础指令

Docker常用基础指令

目录 1 前言 2 常用指令 2.1 获取帮助 2.2 拉取镜像到本地 2.3 对本地镜像进行打包 2.4 对本地镜像的删除 2.5 通过tar包加载本地镜像 2.6 查看所有镜像 2.7 创建新的容器 2.8 查看容器 2.9 停止容器运行 2.10 运行容器 2.11 删除容器 2.12 查看容器日志 2.13 进…
阅读更多...
arcgis 栅格数据处理2——栅格转地级市（栅格转矢量图）

arcgis 栅格数据处理2——栅格转地级市（栅格转矢量图）

1. 获取空间分析权限&#xff08;解决无法执行所选工具问题&#xff09; 选中“自定义”中的“扩展模块” 在弹出的模块中选中能选的模块&#xff0c;此处需要选择“spatial analysis”以进行下一步分析 3. 将栅格数据转为整数型&#xff08;解决无法矢量化&#xff09; 选…
阅读更多...
STM32 HAL库RTC复位丢失年月日的解决办法

STM32 HAL库RTC复位丢失年月日的解决办法

STM32 HAL库RTC复位丢失年月日的解决办法 0.前言一、实现方式1.CubeMX配置&#xff1a;2.MX_RTC_Init()函数修改2.编写手动解析函数 二、总结 参考文章&#xff1a;stm32f1 cubeMX RTC 掉电后日期丢失的问题 0.前言 最近在使用STM32F103做RTC实验时&#xff0c;发现RTC复位后时…
阅读更多...
浏览器工作原理与Javascript高级（前后端异步）

浏览器工作原理与Javascript高级（前后端异步）

总体介绍 浏览器运行是多进程的&#xff0c;包括主进程、渲染进行、网络进程和GPU进程等等 解析HTML时(渲染进程)用到两大引擎&#xff0c;一个是渲染引擎&#xff08;用于渲染页面&#xff09;、一个是JS引擎用于解析JS代码。又JS引擎运行是单线程的&#xff0c;所以渲染和 …
阅读更多...
基于Kronig-Penney能带模型的MATLAB求解与仿真

基于Kronig-Penney能带模型的MATLAB求解与仿真

目录 1.程序功能描述 2.测试软件版本以及运行结果展示 3.核心程序 4.本算法原理 5.完整程序 1.程序功能描述 基于Kronig-Penney能带模型的MATLAB求解与仿真.综合利用 MATLAB提供的求解常微分方程、矩阵行列式、代数表达式化简及绘图等函数 ,可使 Kronig-Penney能带模型分析…
阅读更多...
Roguelike游戏经久不衰，到底有什么魔力?

Roguelike游戏经久不衰，到底有什么魔力?

Roguelike作为诸多游戏类型之一&#xff0c;多年来一直蓬勃发展&#xff0c;深受广大玩家们的喜爱。 近几年游戏市场Roguelike游戏爆款频出&#xff0c;从经典的《以撒的结合》开始&#xff0c;到《杀戮尖塔》、《哈迪斯》、《元气骑士》&#xff0c;再到现在的爆款《吸血鬼幸存…
阅读更多...
生成对抗网络 (GAN)

生成对抗网络 (GAN)

生成对抗网络&#xff08;Generative Adversarial Networks&#xff0c;GAN&#xff09;是由Ian Goodfellow等人在2014年提出的一种深度学习模型。GAN由两部分组成&#xff1a;一个生成器&#xff08;Generator&#xff09;和一个判别器&#xff08;Discriminator&#xff09;&…
阅读更多...
【ARM Coresight | AMBA BUS | Cache | CoreLink | GCC 专栏导读】

【ARM Coresight | AMBA BUS | Cache | CoreLink | GCC 专栏导读】

请阅读【嵌入式开发学习必备专栏 】 文章目录 1. ARM Coresight SoC-400/SoC-600 专栏导读目录1.1 ARM Coresight 专题1.1.1 Performance Profiling1.1.2 ARM Coresight Debug 工具系列1.1.2.1 ARM DS5 系列1.1.2.2 劳特巴赫 Trace32 系列1.1.2.3 JTAG OpenOCD 系列 1.2 ARM Ca…
阅读更多...
这是谁的女儿？其母亲早已红过头了，现在小小年纪的她也爆红网络，没想到吧？

这是谁的女儿？其母亲早已红过头了，现在小小年纪的她也爆红网络，没想到吧？

这是谁的女儿&#xff1f;其母亲早已红过头了&#xff0c;现在小小年纪的她也爆红网络&#xff0c;没想到吧&#xff1f; 原来&#xff0c;作母亲的她在红极一时后似乎沉寂了下来&#xff0c;没想到她11岁的女儿近年来也在社交媒体上走红&#xff0c;她为何也成了小网红呢&…
阅读更多...
git分布式管理-头歌实验搭建Git服务器

git分布式管理-头歌实验搭建Git服务器

一、Git服务器搭建 任务描述 虽然有提供托管代码服务的公共平台&#xff0c;但是对一部分开发团队来说&#xff0c;为了不泄露项目源代码、节省费用及为项目提供更好的安全保护&#xff0c;往往需要搭建私有Git服务器用做远程仓库。Git服务器为团队的开发者们&#xff0c;提供了…
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023