基于Matlab的二阶电路的动态电路分析!

算法思路

1.一阶电路

1.只有电感:先得到电感所在端口的戴维宁等效电路;
接着代入公式计算,得到 i-t图像;
2.只有电容:先得到电容所在端口的戴维宁等效电路;
接着代入公式计算,得到u-t图像

2二阶电路

1.电感电容串联:将电感电容当成一个整体得出对应端口的戴维宁等效电路;接着代入公式,得到u-t图像;
2.电感电容并联:将电感电容当成一个整体得出对应端口的诺顿等效电路;接着代入公式,得到i-t图像。

一阶电路

//将初始电压值和电流值赋值给UStart和IStartUStart=app.UStartEditField.Value;IStart=app.IStartEditField.Value;
//如果是零输入 的话,电路里没电源自己加一个(即判断是否有源,如果无源执行if语句如下)if(app.VNumber==0&&app.INumber==0)NoSource=true;app.INumber=app.INumber+1;NewI=zeros(1,4);NewI(1,1)=app.INumber;NewI(1,2)=2;if(app.CNumber~=0)NewI(1,3)=app.CNP(1,3);NewI(1,4)=app.CNP(1,4);elseif(app.LNumber~=0)NewI(1,3)=app.LNP(1,3);NewI(1,4)=app.LNP(1,4);endapp.INP=[app.INP;NewI];  end
//只有电容的情况
 if(app.CNumber~=0&&app.LNumber==0)app.RNP2=[app.RNP;app.CNP];%将电容矩阵加到电阻矩阵中DynamicCircuit(app,app.RNP2);if(NoSource==true)syms  t R C ustart;
%求解微分方程app.U=dsolve('R*C*Du+u=0','u(0)=ustart','t');ustart=UStart;R=app.EqualResistance;C=app.CNP(1,2);app.U(t)=subs(app.U);tao=app.EqualResistance*app.CNP(1,2);%绘图fplot(app.UIAxes,app.U,[0,7*tao]);else%有源的情况syms  t R Vs C ustart;app.U=dsolve('R*C*Du+u=Vs','u(0)=ustart','t');ustart=UStart;R=app.EqualResistance;Vs=app.EqualVoltage;C=app.CNP(1,2);app.U(t)=subs(app.U);tao=app.EqualResistance*app.CNP(1,2);app.TAOEditField.Value=tao;fplot(app.UIAxes,app.U,[0,7*tao]);end
//只有电感的情况
elseif(app.CNumber==0&&app.LNumber~=0)app.RNP2=[app.RNP;app.LNP];%将电容矩阵加到电阻矩阵中%调用函数 计算动态电路中的戴维宁DynamicCircuit(app,app.RNP2);if(NoSource==true)syms  t R istart L;app.I=dsolve('(L/R)*Di+i=0','i(0)=istart','t');R=app.EqualResistance;istart=IStart;L=app.LNP(1,2);app.I(t)=subs(app.I);tao=L/R;app.TAOEditField.Value=tao;fplot(app.UIAxes,app.I,[0,7*tao]);elsesyms  t R Is L istart;app.I=dsolve('(L/R)*Di+i=Is','i(0)=istart','t');istart=IStart;R=app.EqualResistance;Is=app.EqualCurrent;L=app.LNP(1,2);app.I(t)=subs(app.I);tao=L/R;app.TAOEditField.Value=tao;%            clear(app.UIAxes);%            cla(app.UIAxes);fplot(app.UIAxes,app.I,[0,7*tao]);endapp.UIAxes.YLabel.String='电感电流/A';app.UIAxes.Title.String='I-t';

二阶电路

//电容电感并联的情况如下
if((app.CNP(1,3)==app.LNP(1,3)&&app.CNP(1,4)==app.LNP(1,4))||(app.CNP(1,4)==app.LNP(1,3)&&app.CNP(1,3)==app.LNP(1,4)))app.RNP2=[app.RNP;app.CNP];%将电容矩阵加到电阻矩阵中DynamicCircuit(app,app.RNP2);if(NoSource==true)syms  t G C L ustart istart;app.I=dsolve('L*C*D2i+G*L*Di+i=0','i(0)=istart','Di(0)=ustart/L','t');istart=IStart;ustart=UStart;G=1/app.EqualResistance;C=app.CNP(1,2);L=app.LNP(1,2);app.I=subs(app.I);%求解特征根%                tao=EqualResistance*app.CNP(1,2);%                app.TAOEditField.Value=tao;tao=C/G;fplot(app.UIAxes,app.I,[0,10*tao]);else%无源情况syms  t G C L Is ustart istart;app.I=dsolve('L*C*D2i+G*L*Di+i=Is','i(0)=istart','Di(0)=ustart/L','t');istart=IStart;ustart=UStart;G=1/app.EqualResistance;Is=app.EqualCurrent;C=app.CNP(1,2);L=app.LNP(1,2);app.I=subs(app.I);%求解特征根%                tao=EqualResistance*app.CNP(1,2);%                app.TAOEditField.Value=tao;tao=C/G;fplot(app.UIAxes,app.I,[0,10*tao]);endapp.UIAxes.YLabel.String='电感电流/A';app.UIAxes.Title.String='I-t';
//电容电感串联的情况如下
  elseif(app.CNP(1,3)==app.LNP(1,3)||app.CNP(1,3)==app.LNP(1,4)||app.CNP(1,4)==app.LNP(1,4)||app.CNP(1,4)==app.LNP(1,3))LC=zeros(1,4);LC(1,2)=app.CNP(1,2)*app.LNP(1,2);if(app.CNP(1,3)==app.LNP(1,3))LC(1,3)=app.CNP(1,4);LC(1,4)=app.LNP(1,4);elseif(app.CNP(1,4)==app.LNP(1,4))LC(1,3)=app.CNP(1,3);LC(1,4)=app.LNP(1,3);elseif(app.CNP(1,3)==app.LNP(1,4))LC(1,3)=app.CNP(1,4);LC(1,4)=app.LNP(1,3);elseif(app.CNP(1,4)==app.LNP(1,3))LC(1,3)=app.CNP(1,3);LC(1,4)=app.LNP(1,4);endapp.RNP2=[app.RNP;LC];%将电容电感的假电阻新矩阵矩阵加到电阻矩阵中DynamicCircuit(app,app.RNP2);if(NoSource==true)syms  t R C L  ustart istart;app.U=dsolve('L*C*D2u+R*C*Du+u=0','u(0)=ustart','Du(0)=istart/C','t');ustart=UStart;istart=IStart;R=app.EqualResistance;C=app.CNP(1,2);L=app.LNP(1,2);app.U=subs(app.U);%求解特征跟%                tao=EqualResistance*app.CNP(1,2);%                app.TAOEditField.Value=tao;tao=2*L/R;fplot(app.UIAxes,app.U,[0,7*tao]);
//统一有源情况如下(以上二阶均无源)
elsesyms  t R C L Us ustart istart;app.U=dsolve('L*C*D2u+R*C*Du+u=Us','u(0)=ustart','Du(0)=istart/C','t');ustart=UStart;istart=IStart;Us=app.EqualVoltage;R=app.EqualResistance;C=app.CNP(1,2);L=app.LNP(1,2);app.U=subs(app.U);%求解特征跟
%                tao=EqualResistance*app.CNP(1,2);
%                app.TAOEditField.Value=tao;tao=2*L/R;     fplot(app.UIAxes,app.U,[0,7*tao]);endapp.UIAxes.YLabel.String='电容电压/V';app.UIAxes.Title.String='U-t';end%撤销电流源  让电流源矩阵少一行if(NoSource==true)app.INP(app.INumber,:)=[];app.INumber=app.INumber-1;endapp.UIAxes.XLabel.String='时间/s';app.TAOEditField.Value=tao;hold(app.UIAxes,"on");grid(app.UIAxes,"on");

人机界面

编者独白
1).本次完善在上次一阶电路基础上增加了诺顿定理的验证,最大功率的验证(下次将代码附上);
2).本次实验将一阶电路的动态电路分析全部解决,将二阶电容电感捆绑型电路全部解决(零状态、零响应、全响应)

基于Matlab的二阶电路的动态电路分析!相关推荐

  1. 逆谐波滤波matlab,基于MATLAB仿真的SPWM逆变电路谐波分析及滤波器设计

    第 7 卷第 3 期 2010 年 9 月 长 沙 理 工 大 学 学 报 ( 自 然 科 学 版 ) Journal of Changsha University of Science and Te ...

  2. 基于matlab的单相pwm逆变电路的仿真研究,基于MATLAB的单相PWM逆变电路的仿真研究...

    基于MATLAB的单相PWM逆变电路的仿真研究 [摘要]逆变电路是PWM控制技术最为重要的应用场合.这里在研究单相桥式PWM逆变电路的理论基础上,采用Matlab的可视化仿真工具Simulink建立单 ...

  3. 毕业设计 matlab led,基于Matlab的单相半波可控整流电路的设计与仿真论文和Matlab源码...

    摘 要 电力电子技术是一门诞生和发展于20世纪的崭新技术,在21世纪仍将以迅猛的速度发展.以计算机为核心的信息科学将是21世纪起主导作用的科学技术之一.本次单相半波可控整流电路设计是基于MATLAB的 ...

  4. matlab三相电路基波图形,毕业设计基于matlab的三相交流调压电路的设计与仿真.doc...

    本科毕业论文(设计) 题 目:基于Matlab的三相交流调压电路的设计与仿真 学 院:_____自动化工程学院_____ 专 业:___电气工程及其自动化___ 姓 名:_______________ ...

  5. PLC电路的稳态特性matlab,基于matlab的线性电路正弦稳态特性分析.doc

    题 目: 基于matlab的线性电路正弦稳态特性分析课题要求 利用matlab强大的图形处理功能,符号运算功能和数值计算功能,实现线性电路正弦稳态特性的仿真波形. 课题内容: 1 对R,L,C三种基本 ...

  6. 基于MATLAB/Simulink的三相感应电机动态数学建模仿真模型

    Three_Phase_Induction_Motor:基于MATLAB/Simulink的三相感应电机动态数学建模仿真模型 仿真条件:MATLAB/Simulink R2015b ID:698065 ...

  7. 基于MATLAB Simulink的三相感应电机动态数学建模仿真模型

    Three_Phase_Induction_Motor:基于MATLAB Simulink的三相感应电机动态数学建模仿真模型. 仿真条件:MATLAB Simulink R2015b YID:5380 ...

  8. 电路仿真matlab实验总结,基于matlab的boost电路仿真的实验报告分析.doc

    基于matlab的boost电路仿真的实验报告分析.doc Boost电路 1.实验名称:基于matlab的boost电路仿真的实验报告分析. 2.实验目的:学习matlab的基础知识和操作: 改变占 ...

  9. matlab降压启动,基于 Matlab 的笼形异步电动机降压启动分析

    原标题:基于 Matlab 的笼形异步电动机降压启动分析 基于 Matlab 的笼形异步电动机降压启动分析 陈滨掖 摘要:三相异步电动机全压启动时的瞬间大电流及启动转矩会对负载造成很大冲击.降压启动的 ...

最新文章

  1. 解决: AOSP 编译AndroidQ preview 失败
  2. PHP常用正则:手机 邮箱 网址
  3. discuz数据从godaddy主机中导出的mysql数据乱码变问号???的解决方法
  4. python向dict里添加_Python有条件地向Dict添加键
  5. 别让6种不良心理偷走你的好人缘
  6. linux纯内核直接用吗,Linux:为啥内核有的变量没有初始化就敢直接使用?
  7. java 排序算法面试题_面试题: java中常见的排序算法的实现及比较
  8. 备份图解 ---mysql 博客
  9. 阻止表单提交,超链接默认行为
  10. 黑马程序员ssm总结[大全版本,有对应pdf+源码](spring->springmvc-->springboot-->maven高级->cloud微服务)
  11. 腾跃英语计算机学院微信公众号,英语四级报名_微信还能这么玩:Geek大学生搭建英语课堂互动系统_沪江英语...
  12. ncm在线转换mp3格式
  13. wc与xargs命令
  14. uniapp canvas生成海报不显示问题
  15. WeUI 简明入门指南
  16. 检测PE文件加壳信息用的特征码
  17. @html.dropdownlist默认多选,Html.DropDownList()赋默认值
  18. 752_LaTeX基础软件使用测试
  19. 将秒转换为小时:分钟:秒
  20. win10 文件服务器 群晖,转帖!手把手教您用win10做NAS:搭配emby,VM虚拟群晖,直通......

热门文章

  1. R语言基础之用R语言绘制各类函数图像
  2. 写一个函数,计算一个整数数组的平均值
  3. 实验8:Problem D: 选举班干部了!
  4. 智能家居:以ZigBee技术实现控制器设计
  5. 整数划分 --- 一个老生长谈的问题 动态规划
  6. 为什么绘声绘影10下载好后用不了
  7. 基于微信小程序的在线考试系统源码
  8. 施乐背包服务器显示代码a0,富士施乐故障维修代码.doc
  9. 餐饮管理系统oracle实训,餐饮管理系统开题报告
  10. 使用“PC定时执行专家”,提高你的工作效率