matlab在电路中的全响应程序
- 格式:doc
- 大小:46.50 KB
- 文档页数:4
MATLAB实验:同步发电机三相短路响
应过程仿真
简介
本实验通过使用MATLAB软件对同步发电机的三相短路响应过程进行仿真。
同步发电机是发电厂中常见的发电设备,了解其短路响应过程对于保证电网的稳定性很重要。
实验过程
1. 确定同步发电机的参数,包括额定功率、额定电压、额定频率、定子、转子等参数。
2. 通过在MATLAB中建立电力系统模型,将同步发电机和其他电网元件(例如传输线、发电机组)连接起来。
3. 设计相应的控制策略,实现短路故障的切除、保护与恢复。
4. 运行仿真程序,模拟同步发电机在发生三相短路故障时的响应过程。
5. 分析仿真结果,包括电流、电压、转速等参数的变化情况,以及系统的稳定性和安全性。
实验目标
通过进行该实验,我们可以了解同步发电机在发生三相短路故障时的响应过程,包括电流和电压的变化情况。
这有助于我们更好地了解电力系统的稳定性和安全性,并为实际电力系统中的故障分析和保护设计提供参考。
结论
通过对同步发电机三相短路响应过程的仿真实验,我们可以获得电流和电压的变化情况,并进一步分析电力系统的稳定性和安全性。
这对于电力系统的运行和维护至关重要,有助于提高电网的可靠性和安全性。
【注意】本文档中的内容仅供参考,详细的实验步骤和具体参数需根据实际情况进行确定和调整。
题目:用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 串联电路的直流激励的零状态响应零状态响应就是电路在零初始状态下(动态元件初始储能为零)由外施激励引起的响应。
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对一阶动态电路的响应分析MATLAB是一种高级的数学计算工具,它凭借着强大的计算能力和友好的编程界面,受到了众多工程师和科学家的欢迎。
在电路分析领域,MATLAB也可以用来模拟电路,分析电路的响应以及优化电路的设计。
对于一阶动态电路的响应分析,可以使用MATLAB来模拟电路、计算电路的响应并进行优化。
首先,我们需要通过MATLAB中的电路分析工具箱(Circuit Analysis Toolbox)来建立一阶动态电路模型。
比如下图所示:通过上图,我们设置了一个简单的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)');```通过上述代码,我们得到了电路的响应图表,如下图所示:通过该图表,我们可以清晰地看到,输入的方波信号被电路滤波后,输出信号成为了一种类似于正弦波的响应,并随着时间的推移逐渐稳定。
课程设计任务书学生姓名: 专业班级: 指导教师: 工作单位: 题目 : 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 工具相等的开发设计。
序言 (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软件的应用,以便以后更加容易地解决各种问题。
用MATLAB 求解线性电路的正弦稳态响应一、 阻抗和导纳.一个含线性电阻、电感和电容等元件,但不含独立元的一端口,但它在角频率为w的正弦电压或电流激励下处于稳态时,端口的电流或电压将是同频率的正弦量,应用向量法,端口的电压向量.U 与电流向量.I 的比值定义为该一端口的阻抗Z 。
Z 的模值|Z|称为阻抗模,它的辐角Z ϕ称为阻抗角.则i u Z ϕϕϕ-=。
如果一端口含有元件R 、L 、C ,则其对应的阻抗为:Z=R,Z=jwL,Z=1/jwc.其导纳为Y=1/Z=1/R,Y=1/jwL,Y=jwc. 则 Z=R+Jx Y=G+Jb对于n 个阻抗串联而成的电路,其等效阻抗为 n eq Z Z Z Z +++= 21各个阻抗的电压分配为 n k U Z Z U eqk k ,,2,1,...........== 同理,对于n 个导纳并联而成的电路,其等效导纳n eq Y Y Y Y +++= 21各个导纳的电流分配为 n k I Y Y I eqk k ,2,1,...........== 例1, RLC 串联电路如图,如R=5Ω、L=3H 、C=0.25F 、tV u s 2cos 210=,求电流以及各元件的电压.用matlab 解为:% 已知: us=10*sqrt(2)*cos(2*t) VR=5; L=3; C=0.25;w=2; Us=10;Zl=j*w*L; Zc=1/(j*w*C);I=Us/(Zl+R+Zc)Uab=I*ZlUbc=I*RUcd=I*Zccompass([Uab,I,Ucd]) %画向量图gtext('Uab','color','b'); %用不同的颜色表示gtext('I','color','m');gtext('Ucd','color','r')abs(I) %取模angle(I)*180/pi %degreeangle(I) % rad%画坐标图t=-2:0.01:6; %横坐标从-2到6,每隔0。
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对曲线进行注解。