《带通滤波器电路分析.doc》由会员分享，可在线阅读，更多相关《带通滤波器电路分析.doc（10页珍藏版）》请在金锄头文库上搜索。
1、带通滤波器电路分析与优化一、实验目的1熟悉利用修正节点法求解电路传递函数的方法；2理解各种电路优化的方法，并能将其灵活应用于具体电路的优化；3掌握利用MATLAB工具来优化电路的方法，并能够用程序来实现优化方法，以利于对电路的研究。二、实验要求分析下面的带通滤波器电路图1. 高通滤波器电路图考虑电路的频率响应： 要求优化上面各个元件参数，使得有以下的理想响应：图2. 理想带通滤波器的幅频响应设计过程可以看成对下面的代价函数进行优化，求最小值：（其中是上图给出的理想的频率响应）可选的优化内容包括：1 在元件值没有约束的情况下寻找最优的元件参数，使得归一化后的频率响应满足上图。2 假设元件的值限制
2、在下面的范围内：,如何选择元件，满足设计要求3 假设元件的取值是离散的，即：1、1.1、1.2、1.3、1.5、1.6、1.8、2.0、2.2、2.4、2.7、3.0、3.3、3.6、3.9、4.3、4.7、5.1、5.6、6.2、6.8、7.5、8.2、9.1此时如何选择元件，满足设计要求4 假设所有元件值有的误差，（简单起见，假设元件值在这个范围内均匀分布）分析此时代价函数的变化（代价函数成为一个概率分布）5 设定一个代价函数变化极限，然后尝试反过来确定允许的元件误差范围三、实验内容1问题分析首先，利用修正节点法求解该电路的传递函数，其传递函数为，当输入为1时，所得到的输出即为传递函数的值
3、。为了求得与理想带通滤波器的幅频响应较一致的滤波器，我们设定了如下代价函数：，由于带通滤波器在滤波器边缘变化较快，而其他地方变化较慢，我们应当对边缘处进行着重考虑。因此，我们最终选定如下代价函数：，其中为频点处的权重，经过基本的筛选和分析，最后采取在之间的权重为15，之间权重为50，之间权重为100，其余频率段的权重为1。2代价函数的求取首先，要求得该电路传递函数的值，可以有两种方法供选择，第一是通过修改第一次大作业的读网表文件来求取该值，此方法通用性较好，可以在分析不同的电路时较为方便的得到结果。第二是直接列出修正节点法方程进行求取，此方法通用性不好，只能针对此电路，但是对于本次所研究的对象
4、，表达较为简单，实现起来较为容易。综合考虑，我们选择用第二种方法来求取。在实现时，首先进行数据转换，即电阻转为电导，电容和电感也进行相应的转换。将转换后的数据保存在xteg中，而修正节点法电路表达式中T矩阵的元素保存在数组T中，通过p=inv(T)*b;求得电路参数的值。为了得到多个频点的值，可用如下方法得到不同频率的传递函数的值：7f_start=10000;f_step=10000;f_end=8000000;for f=f_start:f_step:f_end s=2*pi*f*j; i=f/10000;H(i)=abs(p(4)/abs(p(1);其中，H(i)为传递函数。代价函数的求
5、取要用到所求得的传递函数和理想幅频响应，根据以上结论，可通过如下方法实现：for i=1:800 w(i)=1;endfor i=1:800 %此处为权重的设置 endendfor i=1:800 a(i)=w(i)*(H(i)-Hteg(i)2;endfor i=1:800 G=G+a(i);end该程序主要实现了权重的设置，和通过a(i)=w(i)*(H(i)-Hteg(i)2;和G=G+a(i);语句进行代价函数的求取。完整的代码见附录fun.m。3在元件值没有约束的情况下寻找最优的元件参数，使得归一化后的频率响应满足要求。1） 方法：此部分的关键在于最优元件参数的求取，使用MATLAB
6、自带的优化工具箱的fminsearch函数进行搜索，使用途中所给的参数做初值，通过如下语句进行搜索：options=optimset(MaxFunEvals,100000,MaxIter,100000);%函数参数设定x=fminsearch(x)fun(x),x0,options);具体的程序见附录中main.m。2） 结果：得到的x可通过调用G=fun(x);求取代价函数的值，得到的最优元件值和幅频响应如下：表1. 无约束优化下最优元件参数元件R1R2C1C2数值20691846945121929068399730304.4590029536225.4344914925元件C3L1L2L3
7、数值21351.34546833400.05927190308.01450174950.0844666000其中，电阻的单位为，电容的单位为pF，电感为H.图3. 无约束优化下的幅频响应3） 结论：通过电路的优化，得到了较为理想的幅频响应，其边缘的陡峭性也较为理想，但是中间出现了一定的纹波，为了使纹波得到改善，我们尝试了在代价函数中将抖动大的频率段加更大的的权重或者取更多的频点，这些方法都能使纹波得到改善。但抖动和上升沿缓慢是不能同时避免的，最后综合考虑了各种因素，选定了该权重。4假设元件的值限制在下面的范围内：，。如何选择元件，满足设计要求1） 方法：思路类似无约束条件，但是需要对元件的取值
8、做一定的约束。设定元件取值的上下限：lb=0.1,0.1,0.01,0.01,0.01,0.01,0.01,0.01;%下限设定ub=106,106,107,107,107,100,100,100;%上限设定选用优化工具箱的有约束条件优化函数fmincon来完成该优化：options=optimset(MaxFunEvals,100000,MaxIter,100000,TolCon,1e-9,TolFun,1e-9);x=fmincon(x)fun(x),x0,A,b,Aeq,Beq,lb,ub,nonlcon,options);2） 结果：得到的x可通过调用G=fun(x);求取代价函数的值
9、，得到的最优元件值和幅频响应如下：表2. 元件参数有约束时的最优元件参数元件R1R2C1C2数值184.28110122931193.880848521873006.2478299541025.76876139846元件C3L1L2L3数值3006.772764190800.5984924568670.14563259380.59823789471图4. 有取值范围约束下的最优幅频响应3） 结论：通过对比图3和图4的结果可以看出来，有约束时的幅频响应不如无约束时的结果，但是两者相差不大，可见此时对传递函数的幅频响应没有太大的影响。但是通过表1和表2的对比可以看出来，此时元件的取值有较大的差别。
10、5假设元件的取值是离散的，即：1、1.1、1.2、1.3、1.5、1.6、1.8、2.0、2.2、2.4、2.7、3.0、3.3、3.6、3.9、4.3、4.7、5.1、5.6、6.2、6.8、7.5、8.2、9.1此时如何选择元件，满足设计要求1） 方法：离散优化的主要想法是模仿模拟退火法，但由于模拟退火法中关于温度的设定更多时候是一个经验值，在此，我们就略作简化。整个离散优值的选取是在连续优值确定的情况之下，依据可选优值附近的离散值作为初值，对每一元素X（电阻电容电感）实施线性搜索迭代，不断在更优值的基础上获取，最终获取在这一搜索范围内的最优离散值。核心实现方式主体如下：for n=1:3
11、 %综合考虑时间与效率设定循环搜索三周 x=y; for k=1:8 for i=1:8 for j=1:24 x(k)=model(j)*gain(i); cur_val=fun1(x); if(cur_valpre_val) pre_val=cur_val; y(k)=x(k); end end end x(k)=y(k);%很重要 end在此值得注意的是，离散最优值的搜索往往会陷入局部最优的僵局而造成无法获取整体结果最优的优化值，对此主要有两种解决办法：1）可选取更加良好的搜索方式2）对于搜索初值首先进行一定的优化，保证足够好的搜索方向和起点。另外，对于优化函数频带特性和纹波要求显然是较
12、为矛盾的要求，因而，我们选择了对纹波要求较高的初值再配合对频带要求较高的目标函数获取兼顾两者的优化值。2） 结果：得到的x可通过调用G=fun(x);求取代价函数的值，得到的最优元件值和幅频响应如下：表3. 元件参数离散时的最优元件参数元件R1R2C1C2数值13000000000016.000000000056000.0000000000270.0000000000元件C3L1L2L3数值11000.00000000000.03300000006.80000000000.2000000000图5. 元件取值离散时的最优幅频响应3） 结论：通过与图3和图4的结果对比可以看出来，元件取值离散时的
13、幅频响应不如前两种情况时的结果，主要体现在幅频响应曲线的陡峭性不够好，但是纹波稍有改善。通过与表1和表2的对比可以看出来，此时元件的取值与前两种情况也有较大的差别。6假设所有元件值有的误差，（简单起见，假设元件值在这个范围内均匀分布）分析此时代价函数的变化（代价函数成为一个概率分布）1） 方法：容差设计的任务是设计电路标称值与分配电路中元器件的参数容差，使电路性能的偏差最小，或者在保证电路性能满足指标要求的条件下允许元器件参数的容差范围最大。统计过程中，我们采用蒙特卡洛分析法，先是使用随机数发生器按元件值的概率分布来选择元件，然后对整个输出量进行统计，得到最后的模拟分析结果。为了方便评估，我们
14、选定无约束优化所得到的结果作为偏差的基准：x1=1.0e+004 *0.00206918469451 0.00219290683997 3.03044590029536 0.02254344914925 2.13513454683340 0.00000592719030 0.00080145017495 0.00000844666000;%无约束优化最佳值主要通过如下语句来实现偏差的随机性：x(j)=x1(j)*0.95+x1(j)*(1.05-0.95)*rand(1);完整的代码见附件main4.m2） 结果：求得代价函数并通过i=0:10:100;hist(G,i);语句做得如下分布图：图6. 代价函数概率分布图3） 结论：由代价函数的概率分布可以看到，若元件值在最优解的基础上有一定的随机误差时，代价函数值将会增大。而且可以看到，代价函数的概率分布在大于10之后呈指数下降，也就是说当元件已取得最优解后，其本身的误差对电路造成极大影响的概率很小，并且系统误差越大，发生的概率越小，这也是我们进行电路优化的意义所在。6设定一个代价函数变化极限，然后尝试反过来确定允许的元件误差范围1） 方法：由容差分析得出
