matlab在电路中的全响应程序

合集下载

MATLAB实验：同步发电机三相短路响应过程仿真

MATLAB实验：同步发电机三相短路响
应过程仿真
简介
本实验通过使用MATLAB软件对同步发电机的三相短路响应过程进行仿真。

同步发电机是发电厂中常见的发电设备，了解其短路响应过程对于保证电网的稳定性很重要。

实验过程
1. 确定同步发电机的参数，包括额定功率、额定电压、额定频率、定子、转子等参数。

2. 通过在MATLAB中建立电力系统模型，将同步发电机和其他电网元件（例如传输线、发电机组）连接起来。

3. 设计相应的控制策略，实现短路故障的切除、保护与恢复。

4. 运行仿真程序，模拟同步发电机在发生三相短路故障时的响应过程。

5. 分析仿真结果，包括电流、电压、转速等参数的变化情况，以及系统的稳定性和安全性。

实验目标
通过进行该实验，我们可以了解同步发电机在发生三相短路故障时的响应过程，包括电流和电压的变化情况。

这有助于我们更好地了解电力系统的稳定性和安全性，并为实际电力系统中的故障分析和保护设计提供参考。

结论
通过对同步发电机三相短路响应过程的仿真实验，我们可以获得电流和电压的变化情况，并进一步分析电力系统的稳定性和安全性。

这对于电力系统的运行和维护至关重要，有助于提高电网的可靠性和安全性。

【注意】本文档中的内容仅供参考，详细的实验步骤和具体参数需根据实际情况进行确定和调整。

用MATLAB对RC、RL电路进行分析.

题目：用MATLAB 对RC 、RL 电路进行分析摘要： MATLAB 是美国Mathworks 公司开发的大型软件包，是MATrix LABoratory 的缩略语。

目前，MATLAB 广泛应用于线性代数、高等数学、物理、电路分析、信号与系统、数字信号处理、自动控制等众多领域，是当前国际上最流行的科学与工程计算的工具软件。

MATLAB 功能强大并且同其它高级语言相比具有语法规则简单、容易掌握、调试方便等特点。

Simulink 是MATLAB 软件的扩展，它是实现动态系统建模和仿真的一个软件包。

MATLAB 具有强大的图形处理功能、符号运算功能和数值计算功能。

其中系统的仿真（Simulink ）工具箱是从底层开发的一个完整的仿真环境和图形界面。

在这个环境中，用户可以完成面向框图系统仿真的全部过程，并且更加直观和准确地达到仿真的目标。

本次主要介绍基于MATLAB 的一阶动态电路特性分析。

关键字：MATLAB ；仿真；图形处理；一阶动态电路。

一. RC 串联电路1.1 RC 串联电路的零输入响应动态电路中无外施激励电源，仅由动态元件初始储能所产生的响应，称为动态电路的零输入响应。

在图1所示的RC 电路中，开关S 打向2前，电容C 充电，U u u C R =+。

当开关S 打向2后，电压C R u u =，电容储存的能量将通过电阻以热能的形式释放出来【2】。

图1 RC 电路的零输入响应电路分析：由图可知 t RC o e R U i 1-=， t RC o C R e U u u 1-== t RC o R e R U R I p 222-==，t RC o C C e R U iu p 22-== 在MATALAB 的M 文件编写以下程序:U0=40;R=10;C=0.5; %输入给定参数U1=10;R1=5;C1=0.5; %输入给定参数t=[0:0.1:10]; %确定时间范围Uc1=U0*exp(-t/(R*C));Uc2=U1*exp(-t/(R*C)); %电容电压值Ur1=U0*exp(-t/(R*C));Ur2=U1*exp(-t/(R*C)); %电阻电压值I1=U0/R*exp(-t/(R*C));I2=U1/R*exp(-t/(R*C)); %计算电流值Pc1=U0^2/R*exp(-2*t/(R*C));Pc2=U1^2/R*exp(-2*t/(R*C)); %电容功率值 Pr1=U0^2/R*exp(-2*t/(R*C));Pr2=U1^2/R*exp(-2*t/(R*C)); %电阻功率值 figuresubplot(5,1,1);plot(t,Uc1,t,Uc2); title('Uc(t)的波形图')subplot(5,1,2);plot(t,Ur1,t,Ur2); title('Ur(t)的波形图')subplot(5,1,3);plot(t,I1,t,I2); title('I(t)的波形图')subplot(5,1,4);plot(t,Pc1,t,Pc2); title('Pc(t)的波形图')subplot(5,1,5);plot(t,Pr1,t,Pr2); title('Pr(t)的波形图')波形仿真图：图2 RC 串联电路零输入响应特性曲线蓝线表示U0=40;R=10;C=0.5情况下的特性曲线绿线表示U1=10;R1=5;C1=0.5情况下的特性曲线1.2 RC 串联电路的直流激励的零状态响应零状态响应就是电路在零初始状态下（动态元件初始储能为零）由外施激励引起的响应。

基于MATLAB的RL并联电路频率响应特性分析

1理论分析1.1 MATLAB简介1.1.1 MATLAB的概况20世纪70年代，美国新墨西哥大学计算机科学系主任Cleve Moler为了减轻学生编程的负担，用FORTRAN编写了最早的MATLAB。

1984年由Little、Moler、Steve Bangert 合作成立了的MathWorks公司正式把MATLAB推向市场。

到20世纪90年代，MATLAB已成为国际控制界的标准计算软件。

MATLAB是矩阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。

除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。

MATLAB的基本数据单位是矩阵，它的指令表达式与数学,工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完相同的事情简捷得多。

MATLAB包括拥有数百个内部函数的主包和三十几种工具(Toolbox).工具包又可以分为功能性工具包和学科工具包.功能工具包用来扩充MATLAB的符号计算,可视化建模仿真,文字处理及实时控制等功能.学科工具包是专业性比较强的工具包,控制工具包,信号处理工具包,通信工具包等都属于此类.开放性使MATLAB广受用户欢迎.除内部函数外,所有MATLAB主包文件和各种工具包都是可读可修改的文件,用户通过对源程序的修改或加入自己编写程序构造新的专用工具包.1.1.2 MATLAB的特点一种语言之所以能如此迅速地普及，显示出如此旺盛的生命力，是由于它有着不同于其他语言的特点，正如同FORTRAN和C等高级语言使人们摆脱了需要直接对计算机硬件资源进行操作一样，被称作为第四代计算机语言的MATLAB，利用其丰富的函数资源，使编程人员从繁琐的程序代码中解放出来。

MATLAB最突出的特点就是简洁。

基于METLAB对一阶动态电路的响应分析

基于METLAB对一阶动态电路的响应分析MATLAB是一种高级的数学计算工具，它凭借着强大的计算能力和友好的编程界面，受到了众多工程师和科学家的欢迎。

在电路分析领域，MATLAB也可以用来模拟电路，分析电路的响应以及优化电路的设计。

对于一阶动态电路的响应分析，可以使用MATLAB来模拟电路、计算电路的响应并进行优化。

首先，我们需要通过MATLAB中的电路分析工具箱（Circuit Analysis Toolbox）来建立一阶动态电路模型。

比如下图所示：![image.png](attachment:image.png)通过上图，我们设置了一个简单的RC电路，其中R为100欧姆，C为100微法，输入信号为1V方波信号，周期为1ms。

使用MATLAB的Circuit Analysis Toolbox将电路参数输入，并建立电路模型。

具体过程如下：```MATLAB% 定义电路参数R = 100; % OhmC = 100e-6; % FVin = 1; % V% 定义电路模型f = @(t,Vout) (-Vout / (R*C)) + (Vin / (R*C));options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);[t, Vout] = ode45(f,[0,10],0,options);```通过上述代码，我们成功地建立了一个RC电路的数学模型，并计算得到了该电路的时域响应，其中输入信号的周期为1ms，最后得到输出电压Vout随时间的变化情况。

这里，我们使用了MATLAB中的ode45()函数，它可以用来求解微分方程。

接下来，为了更加清晰地观察电路响应，我们可以将电路的输入和输出信号绘制成图表。

代码如下：```MATLAB% 绘制输出信号图表plot(t,Vout,'LineWidth',2);grid on;title('RC Circuit Response');xlabel('Time (s)');ylabel('Voltage (V)');```通过上述代码，我们得到了电路的响应图表，如下图所示：![image-2.png](attachment:image-2.png)通过该图表，我们可以清晰地看到，输入的方波信号被电路滤波后，输出信号成为了一种类似于正弦波的响应，并随着时间的推移逐渐稳定。

Matlab课程设计-基于Matlab的RC串联电路频率响应特性分析

课程设计任务书学生姓名: 专业班级: 指导教师: 工作单位: 题目 : Matlab 应用课程设计-基于 Matlab 的 RC 串联电路频率响应特性分析初始条件:1. Matlab6.5以上版本软件;2. 先修课程:电路原理等;3. 2, 0.5R C F =Ω=。

要求完成的主要任务:(包括课程设计工作量及其技术要求,以及说明书撰写等具体要求1、在 RC 串联电路中 , 求该电路的频率响应10( ( ( C U j H j U j ωωω=和 20( ( ( R U j H j U j ωωω=,并绘出其特性曲线; 2、画出程序设计框图,编写程序代码,上机运行调试程序,记录实验结果 (含计算结果和图表等 ,并对实验结果进行分析和总结;3、课程设计说明书按学校统一规范来撰写,具体包括:⑴目录; ⑵理论分析;⑶程序设计; ⑷程序运行结果及图表分析和总结; ⑸课程设计的心得体会(至少 500字 ;⑹参考文献(不少于 5篇。

时间安排:周一、周二查阅资料,了解设计内容;周三、周四程序设计,上机调试程序;周五、整理实验结果,撰写课程设计说明书。

指导教师签名: 年月日系主任(或责任教师签名: 年月日目录1.Matlab 软件简介 (1)2.RC 串联电路频率响应特性分析 (2)3. 程序设计 (3)4. 程序运行结果及结果分析 (5)5. 课程设计的心得会 (7)6. 参考文献 (8)1.Matlab 软件简介1.1Matlab 语言的历史70年代后期 , 身为美国 New Mexico大学计算机系系主任的 Cleve Moler发现学生用 FORTRAN 编写接口程序很费时间 , 于是他开始自己动手 , 利用业余时间为学生编写 EISPACK 和 LINPACK 的接口程序。

Cleve Moler给这个接口程序取名为 Matlab 。

1984年, 为了推广 Matlab 在数值计算中的应用, Cleve Moler、 Johon Little 等正式成立了 Math works公司, 从而把 Matlab 推向市场, 并开始了对 Matlab 工具相等的开发设计。

MATLAB在《电路理论》课程教学中的应用

图５ＲＣ一阶电路图６一阶电路全响应波形解：根据一阶电路的三要素法可知， “ ，０）＝ “ （ｃｏ）＋［ “ （Ｏ＋）一
一
“ （（ｃｏ）］
ｃｌｃ；
，编制简单的ＭＡＴＬＡＢ程序如下：
利用板书讲解时，为了能体现出不同频率的正弦波叠加逼近非正弦周期信号的效果，往往先画出一次谐波的图形，与方波信号做对比，或是多画出几次谐波叠加的效果与方波信号做对比，此方法费时且不直观。此时在课堂教学中就可用ＭＡＴＬＡＢ来验证，首先通过简单的编程，画出方波信号的一次谐波波形，如图１所示，图２为方波信号前９次谐波叠加的效果。利用ＭＡＴＬＡＢ的绘图命令ｐｌｏｔ，画出方波信号各次谐波叠加后的效果，如图３所示。
１／２＋１／２］；＞＞Ｂ＝［一１５／３；１５／３ — ５＋１０；５１；
＞＞Ｕ＝Ａ＼Ｂ
图３方波信号各次谐波叠加二维图图４方波信号各次谐波叠加ｊ维图通过以上波形图可以很容易的看出，谐波次数取得越高，各次谐波叠加的波形就越接近方波。为了能更直观的观察、对比各次谐波叠加的效果，还可利用ＭＡＴＬＡＢ的绘图命令画出叠加效果的三维图形，如图
提供了一种方法。
［关键词］ＭＡＴＬＡＢ电路
一
辅助分析
辅助计算以一阶动态电路为例，运用ＭＡＴＬＡＢ对电路中的响应进行分析计算，并且可利用ＭＡＴＬＡＢ的绘图功能更加直观的观察各响应的变化情

基于MATLAB的电路频率响应分析-施晨程

3.2MATLAB
3.2.1
（1）程序设计框图：
（2）获取函数的频率响应的编程：
Clear
fz=[0.2 0.2 0];% 表达式分子的系数向量
fm=[0.35 3.23];% 表达式分母的系数向量
w=0:0.1:50;%设定频率的变化范围
g=freqs(fz,fm,w);%求其频率响应
x=angle(g);%求其相频响应
先求： =1.5+ =1.5+ ，又因为给出了电源的ω=1，所以 =1.5+ =1.5061+j0.1025=1.5061 ，所以根据欧姆定律得到 = = =37.5596
然后求出： = - * =3.85 ，于是：
可以求出： = =38.5
又根据基尔霍夫电流定律可以求出： = - =1.355
3
3.1
（1）对于第一小题，我采用编程的方法来求得频率响应的曲线，具体的程序设计思路如下：首先输入H的关于频率的表达式，然后分别列出分子分母的系数行列式，调用angle函数来计算电路的相频特性，然后再调用abs函数来计算电路的幅频特性，调用freqs函数计算频率特性，最后用plot函数来描绘曲线。
（2）对于第二小题，由第一小题求出了和函数，根据其关系式，将ω=1带入可以得到和的向量值，接着可以求出所要求的各个值，所得到的表达式用matlab化简，最后得出结果。
关键词：电路原理，传递函数，MATLAB，向量图
基于MATLAB的电路频率响应分析
1
1.1
本次基础强化训练，我将利用MATLAB软件进行对频率响应电路进行原理分析、建模以及必要的推导和可行性分析。通过这一过程，掌握MATLAB软件的基本操作，体会MATLAB在实际中的应用价值，并且在以后的学习工作中利用MATLAB为自己服务，解决自己遇到的问题。

MATLAB 在电路原理中的应用

序言 (2)1 基础强化训练的要求 (3)2 基础强化训练的内容 (3)2.1 所选择要解决的问题 (3)2.2 问题解决分析 (4)2.2.1 设计分析 (4)2.2.2 程序编写 (5)2.2.3 程序分析 (6)2.3 用simulink仿真 (6)2.3.1 仿真 (6)2.3.2 仿真过程中遇到的问题及分析解决 (8)3 心得体会 (9)参考文献 (10)MATLAB 在电路原理中的应用序言MATLAB是由美国MathWorks公司推出的用于数值计算和图形处理计算系统环境,除了具备卓越的数值计算能力外,它还提供了专业水平的符号计算,文字处理,可视化建模仿真和实时控制等功能.MATLAB的基本数据单位是矩阵,它的指令表达式与数学,工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言简捷得多.MATLAB是国际公认的优秀数学应用软件之一。

由于MATLAB的基本数据单位是矩阵，所以其可以运用于线性代数、电路基础等课程中解决复杂的矩阵方程，而且用MATLAB编程比用C编程更为灵活简单，所以许多电路原理上需要解复杂矩阵方程的问题用MATLAB可轻松地得到解决。

本次基础强化训练的目的即在此。

通过学习了解MATLAB软件的运用，将MATLAB的功能运用于解决电路题目中，在此过程中体会到MATLAB软件功能的强大及运用起来的方便。

学习运用MATLAB 解决电路基础问题的同时还可熟悉电路基础的知识，灵活运用线性代数的知识于电路基础之中。

MATLAB软件在今后的工作中会对我们有很大的帮助，所以现在了解它，有利于我们日后接触到这类问题时可很快得想出解决问题的方案。

通过本环节，巩固和深化已学课程的知识，培养学生综合运用这些知识，分析和解决实际问题，逐步树立正确的设计思想；培养严谨认真的科学态度和严谨务实的工作作风。

1 基础强化训练的要求此次基础强化训练主要以学习MATLAB软件为目的，了解相关问题的处理方法和步骤，基本掌握MATLAB的软件运用，让同学们更多的接触学习MATLAB软件的应用，以便以后更加容易地解决各种问题。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

A=[1 0 0 1 1 0;0 0 1 0 -1 1;0 1 0 -1 0 -1] %输入相关矩阵A
Is=[2;0;-5;0;0;0] %输入电流源列向量
Us=[0;3;0;-2;0;0] %输入电压源列向量
Y=[1 0 0 0 0 0 ;-2 1/3 0 0 0 0 ;0 0 1/4 0 0 0;0 0 0 1/5 0 1/3 ;0 0 0 0 1/2 0;0 0 0 0 0 1/8] %输入导纳矩阵Y
Z=inv(Y)%求解阻抗矩阵Z
Yn=A*Y*A' %求解节点导纳矩阵Yn
Un=inv(Yn)*(A*Is-A*Y*Us) %求解节点电压Un
U=A'*Un %求解之路电压U
I=Y*(U+Us)-Is %求解支路电流I
结果如下：
Z =
1.0000 0 0 0 0 0
6.0000 3.0000 0 0 0 0
0 0 4.0000 0 0 0
0 0 0 5.0000 0 -13.3333
0 0 0 0 2.0000 0
0 0 0 0 0 8.0000
Yn =
1.7000 -0.1667 -0.5333
-0.5000 0.8750 -0.1250
-2.2000 -0.4583 0.9917
Un =
-7.8339
-13.7855
-25.1626
U =
-7.8339
-25.1626
-13.7855
17.3287
5.9516
11.3772
I =
-9.8339
8.2803
1.5536
6.8581
2.9758
1.4221
A=input('相关矩阵A=:')
Is=input('电流源Is= :')
Us=input('电压源Us=:' )
display ('是否含有受控源？是（a=1)否（a=0)')
a=input('输入a=:' )
if(a==0) %判断语句
Z=input('阻抗矩阵Z=:')
Y=inv(Z) %求解导纳矩阵Y
else (a==1)
Y=input('导纳矩阵Y=:')
Z=inv(Y) %求解阻抗矩阵Z
end
Yn=A*Y*A' %求解节点导纳矩阵Yn
Un=inv(Yn)*(A*Is-A*Y*Us) %求解节点电压Un
U=A'*Un %求解之路电压U
I=Y*(U+Us)-Is %求解支路电流I
含受控源的运行结果如下
相关矩阵A=:[1 0 0 1 1 0;0 0 1 0 -1 1;0 1 0 -1 0 -1]
电流源Is= :[2;0;-5;0;0;0]
电压源Us=:[0;3;0;-2;0;0]
是否含有受控源？是（a=1)否（a=0)
输入a=:1
导纳矩阵Y=:[1 0 0 0 0 0 ;-2 1/3 0 0 0 0 ;0 0 1/4 0 0 0;0 0 0 1/5 0 1/3;0 0 0 0 1/2 0;0 0 0 0 0 1/8]
Z =
1.0000 0 0 0 0 0
6.0000 3.0000 0 0 0 0
0 0 4.0000 0 0 0
0 0 0 5.0000 0 -13.3333
0 0 0 0 2.0000 0
0 0 0 0 0 8.0000
Yn =
1.7000 -0.1667 -0.5333
-0.5000 0.8750 -0.1250
-2.2000 -0.4583 0.9917
Un =
-7.8339
-13.7855
-25.1626
U =
-7.8339
-25.1626
-13.7855
17.3287
5.9516
11.3772
I =
-9.8339
8.2803
1.5536
6.8581
2.9758
1.4221
A=[1 -1 0;-1 0 1] %输入相关矩阵A
Is=[-2;-1;0] %输入电流源列向量
Us=[2;0;-3] %输入电压源列向量
Z=[2 0 0;0 0.8 0;0 0 3] %输入阻抗矩阵Z
Y=inv(Z) %求解导纳矩阵Y
Yn=A*Y*A'%求解节点导纳矩阵Yn
Un=inv(Yn)*(A*Is-A*Y*Us) %求解节点电压Un U=A'*Un%求解支路电压U
I=Y*(U+Us)-Is%求解支路电流I
运行结果如下：
Yn =
1.7500 -0.5000
-0.5000 0.8333
Un =
0.2759
4.9655
U =
-4.6897
-0.2759
4.9655
I =
0.6552
0.6552
0.6552
(要有表达式RU194页)
A=input('相关矩阵A=:')
Is=input('电流源Is= :')
Us=input('电压源Us=:' )
display ('是否含有受控源？是（a=1)否（a=0)')
a=input('输入a=:' )
if(a==0) %判断语句
Z=input('阻抗矩阵Z=:')
Y=inv(Z) %求解导纳矩阵Y
else (a==1)
Y=input('导纳矩阵Y=:')
Z=inv(Y) %求解阻抗矩阵Z
end
Yn=A*Y*A' %求解节点电压Un
Un=inv(Yn)*(A*Is-A*Y*Us) %求解节点电压Un %求解支路电流I
U=A'*Un %求解支路电压U
I=Y*(U+Us)-Is %求解支路电流I
不含受控源的运行结果如下：
是否含有受控源？是（a=1)否（a=0)
输入a=:0
Yn =
1.7500 -0.5000
-0.5000 0.8333
Un =
0.2759
4.9655
U =
-4.6897
-0.2759
4.9655
I =
0.6552
0.6552
0.6552
(要有表达式RU194页)
封面设计题目：电路信号处理综合课程设计课程设计任务书
一、解决问题
1．大规模电路的计算机辅助分析
2．动态电路的计算机辅助分析
3．谐振电路频率响应的分析
设计说明书写作：
一、课程设计任务书（第1页）四设计总结五参考文献
二、目录（第3页）
三、设计第一题
A程序B例子C最后运行结果D对曲线进行注解。