资源描述：
《MATLAB,SIMULINK在阻抗匹配网络中的分析与应用》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、武汉理工大学《基础技能强化训练》课程设计说明书MATLAB简介MATLAB是由美国MathWorks公司推出的用于数值计算和图形处理计算系统环境,除了具备卓越的数值计算能力外,它还提供了专业水平的符号计算,文字处理,可视化建模仿真和实时控制等功能.MATLAB的基本数据单位是矩阵,它的指令表达式与数学,工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言简捷得多.MATLAB是国际公认的优秀数学应用软件之一.概括地讲,整个MATLAB系统由两部分组成,即MATLAB内核及辅助工具箱,两者的调用构成了MATLAB的强大功能.MA
2、TLAB语言以数组为基本数据单位,包括控制流语句,函数,数据结构,输入输出及面向对象等特点的高级语言,它具有以下主要特点:1)运算符和库函数极其丰富,语言简洁,编程效率高,MATLAB除了提供和C语言一样的运算符号外,还提供广泛的矩阵和向量运算符.利用其运算符号和库函数可使其程序相当简短,两三行语句就可实现几十行甚至几百行C或FORTRAN的程序功能.2)既具有结构化的控制语句(如for循环,while循环,break语句,if语句和switch语句),又有面向对象的编程特性.3)图形功能强大它既包括对二维和三维数据可视化,图像处理,动画制作等高层次的绘图
3、命令,也包括可以修改图形及编制完整图形界面的,低层次的绘图命令.4)功能强大的工具箱工具箱可分为两类:功能性工具箱和学科性工具箱.功能性工具箱主要用来扩充其符号计算功能,图示建模仿真功能,文字处理功能以及与硬件实时交互的功能.而学科性工具箱是专业性比较强的,如优化工具箱,统计工具箱,控制工具箱,小波工具箱,图象处理工具箱,通信工具箱等.7武汉理工大学《基础技能强化训练》课程设计说明书1含耦合电感电路的计算电路图如图1所示，已知C1=-j3Ω，L1=j8Ω，L2=j6Ω，R1=4Ω，R2=5Ω。若Us=100v，L1、L2为反向耦合，互感阻抗为j2Ω；求I1
4、，I2的值。I2I1图1非线性电路1.1对含耦合电感电路进行原理分析：以上电路中负载已知，且已知L1、L2反向耦合，让我们求出电流I1、I2的值。这里我们可以用回路电流法直接将其求出。我们以L1、L2、R1、C1三个元件组成的电路为树，以V1那条线为一支路，以R2负载为另外一条支路，从而组成两回路，根据这两条回路我们可以求得两个回路的回路电流，根据回路电流可求得两个回路的回路电流，即I1、I2.1.2建模：I1、I2如图已设出，在两个回路中有一个公用的负载L2，由于在两个回路中如果流经两个回路中公用的阻抗的电流方向相同就去取正，如果流经的电流方向相反则为负
5、，在本题中流经L2的两个电流是相反的所以应取负，在含耦合电感的电路中，电感同向耦合互感阻抗取正，电感反向耦合互感阻抗取负，这里电感反向耦合所以取负。因而根据这两个7武汉理工大学《基础技能强化训练》课程设计说明书回路我们可以写出两个相关的方程；即我们可以把这两个方程转换成矩阵，设A矩阵为I1和I2的序数，B矩阵为右边的【100，0】，设I1和I2就为要求的未知量I，由此我们可以写出一个矩阵方程AI=B；两边同时左乘A的逆，就可以得到I=A的逆乘以B.在A矩阵中在I矩阵中在B矩阵中;通过计算我们就可以求出两支路的电流电流的值，即I1、I2。我们求得I1=20.
6、2621+1.2398j；I2=8.2182+2.8339j1.3应用MATLAB对上面的题目编程clear;clc;V1=100:R1=4;R2=5;C1=-1i*3;L1=1i*8;L2=1i*6;%为给定原件赋值A=[R1+C1+L2-1i*2-L2;-1i*2-L2L1+L2+1i*2*2+5];%对A矩阵各元素赋值B=[V1;V1*0];%对B矩阵各元素赋值I=inv(A)*B;%方程求解I1=I(1)%解出电流值I1I2=I(2)%接触电流值I2运行程序后得到的结果如图1-2所示：7武汉理工大学《基础技能强化训练》课程设计说明书图1-2运行结果
7、1.4绘图通过编程我们就可以求出I1、I2的值，但是当负载一定时流经负载的电流与电源电压有什么关系呢？为此我们可以运用MATLAB强大的绘图功能，通过设置电源电压的变化来得到各组不同的流经负载的电流值。为此我们可以这样编程：clear;clc;V1=50:100;R1=4;R2=5;C1=-1i*3;L1=1i*8;L2=1i*6;A=[R1+C1+L2-1i*2-L2;-1i*2-L2L1+L2+1i*2*2+5];B=[V1;V1*0];I=inv(A)*B;I1=I(1,V1-50+1);I2=I(2,V1-50+1);subplot(2,1,1);
8、plot(V1,I1);xlabel('V1');ylabel('
