一、实验目的
实现一个可分析任意还有动态元件电路的暂态响应过程的程序。通过程序设计和开发,深刻理解暂态响应这一过程的物理意义,通过数学表达式描述电路结构和响应过程,最终通过程序设计语言实现一个具有良好可视化界面的分析软件
二、设计思路
分析此次实验可以看出,电路的暂态响应过程是实验所要展现的重点,电路的构建过程可由之前线性电路的编程得到。因此,我们简化了电路模型的输入,构建戴维宁和诺顿等效电路代替实际电路模型。在此基础上,通过KCL、KVL的约束条件,以及动态电路元件自身特性和初始(特定)状态,求解电路方程的解。最后根据需要将电路的响应动态展示出来。
三、关键步骤:
1、在一阶电路中,使用pause函数加for循环动态描述电路参数的响应过程
2、二阶电路中过程较为复杂,学习起来也较为困难。所以我将四个参数放在了同一张图里便于分析比较。同时,为了美观性,将动态改为了静态展示。
3、人机交互界面使用了MATLAB中自带的guide,通过menu实现了用户对电路模型的建立。通过while循环和多个标志参数实现子菜单的管理,从而实现了一套程序的多次重载。
四、亮点:
1、功能强,在题目要求的基础上,增加了一阶电路的全响应,以及二阶电路的零状态和零输入分析功能;
2、注重细节效果,一阶电路中多次实验,为了更好的展示效果,最终确定步数叠加100tao次;
3、通过函数减少了代码量,使结构更加简练便捷;
软件介绍:MATLAB——一个很常规但是功能强大的软件
在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等,主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。
MATLAB的基本数据单位是矩阵,它的指令表达式与数学、工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的事情简捷得多,并且MATLAB也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C,FORTRAN,C++,JAVA的支持。
五、效果展示
主界面
子页
六、代码
下面展示一些 源码
。
syms Us U0 Req C Is I0 L high;
RoadFlag=1;
StateFlag=0;
init()
while RoadFlag==1
RoadType=menu('电路类型','RC','RL','RCL','恢复出厂设置','结束');
switch RoadType
case 1
StateFlag=1;
while StateFlag==1
StateFlag=menu ('选择输入类型','等效电压源','初始电容电压','等效电阻','电容','零输入响应','零状态响应','全响应','退出');
switch StateFlag
case 1;
Us=input('请输入Us的值:');
case 2
U0=input('请输入U0的值:');
case 3
Req=input('请输入Req的值:');
case 4
C=input('请输入C的值:');
case 5 %零输入响应
tao=Req*C;
x1=0;dx=10*tao;%确定起始点和增量
x2=x1+dx;%确定画图时的横坐标终止值x1
y1=0;dy=U0;%确定起始点和增量
y2=y1+dy+dy/5;%确定画图时的横坐标终止值x1
figure('NumberTitle','off','Name','一阶RC电路零输入动态分析');
hold on
axis([x1 x2 y1 y2]);
grid;
title('Uc-t');
for t=x1:0.01:x2;
Uc=U0*(exp(-t/tao));
plot(t,Uc,'k.','markersize',15);
fprintf('时间=%d 电容电压=%d\n',t,Uc);
pause(0.01);
end
fprintf('U0=%dV,Req=%dΩ,C=%d\n',U0,Req,C);
fprintf('时间常数=%ds\n' ,tao);
case 6 %零状态响应
tao=Req*C;
x1=0;dx=10*tao;
x2=x1+dx;
y1=0;dy=Us;
y2=y1+dy+dy/5;
figure('NumberTitle','off','Name','一阶RC电路零状态动态分析');
hold on
axis([x1 x2 y1 y2]);
grid;
title('Uc-t');
for t=x1:0.01:x2;
Uc=Us*(1-exp(-t/tao));
plot(t,Uc,'k.','markersize',15);
fprintf('时间=%d 电容电压=%d\n',t,Uc);
pause(0.01);
end
fprintf('Us=%dV,Req=%dΩ,C=%d\n',Us,Req,C);
fprintf('时间常数=%ds\n' ,tao);
case 7 %全响应
tao=Req*C;
x1=0;dx=10*tao;
x2=x1+dx;
y1=0;dy=Us+U0;
y2=y1+dy+dy/5;
figure('NumberTitle','off','Name','一阶RC电路全响应分析');
hold on
axis([x1 x2 y1 y2]);
grid;
title('Uc-t');
for t=x1:0.01:x2;
Uc=U0*(exp(-t/tao))+Us*(1-exp(-t/tao));
plot(t,Uc,'k.','markersize',15);
fprintf('时间=%d 电容电压=%d\n',t,Uc);
pause(0.01);
end
fprintf('Us=%dV,U0=%dV,Req=%dΩ,C=%d\n',Us,U0,Req,C);
fprintf('时间常数=%ds\n' ,tao);
case 8
StateFlag=0;
RoadFlag=0;
end
end
case 2
StateFlag=2;
while StateFlag==2
StateFlag=menu ('选择输入类型','等效电流源','初始电感电流','等效电阻','电感','零输入响应','零状态响应','全响应','退出');
switch StateFlag
case 1;
Is=input('请输入Is的值:');
case 2
I0=input('请输入I0的值:');
case 3
Req=input('请输入Req的值:');
case 4
L=input('请输入L的值:');
case 5 %零输入响应
tao=Req/L;
x1=0;dx=10*tao;%确定起始点和增量
x2=x1+dx;%确定画图时的横坐标终止值x1
y1=0;dy=I0;%确定起始点和增量
y2=y1+dy+dy/5;%确定画图时的横坐标终止值x1
figure('NumberTitle','off','Name','一阶RL电路零输入动态分析');
hold on
axis([x1 x2 y1 y2]);
grid;
title('Il-t');
for t=x1:0.01:x2;
Il=I0*(exp(-t/tao));
plot(t,Il,'k.','markersize',15);
fprintf('时间=%d 电感电流=%d\n',t,Il);
pause(0.01);
end
fprintf('U0=%dV,Req=%dΩ,C=%d\n',I0,Req,L);
fprintf('时间常数=%ds\n' ,tao);
case 6
tao=Req/L;
x1=0;dx=10*tao;
x2=x1+dx;
y1=0;dy=Is;
y2=y1+dy+dy/5;
figure('NumberTitle','off','Name','一阶RL电路零状态动态分析');
hold on
axis([x1 x2 y1 y2]);
grid;
title('Il-t');
for t=x1:0.01:x2;
Uc=Il*(1-exp(-t/tao));
plot(t,Il,'k.','markersize',15);
fprintf('时间=%d 电容电压=%d\n',t,Il);
pause(0.01);
end
fprintf('I0=%dV,Req=%dΩ,L=%d\n',I0,Req,L);
fprintf('时间常数=%ds\n' ,tao);
case 7 %全响应
tao=Req/L;
x1=0;dx=10*tao;
x2=x1+dx;
y1=0;dy=Is+I0;
y2=y1+dy+dy/5;
figure('NumberTitle','off','Name','一阶RL电路全响应分析');
hold on
axis([x1 x2 y1 y2]);
grid;
title('Il-t');
for t=x1:0.01:x2;
Il=I0*(exp(-t/tao))+Is*(1-exp(-t/tao));
plot(t,Il,'k.','markersize',15);
fprintf('时间=%d 电感电流=%d\n',t,Il);
pause(0.01);
end
fprintf('Is=%dV,I0=%dV,Req=%dΩ,L=%d\n',Is,I0,Req,L);
fprintf('时间常数=%ds\n' ,tao);
case 8
StateFlag=0;
RoadFlag=0;
end
end
case 3
StateFlag=3;
while StateFlag==3
StateFlag=menu ('选择输入类型','等效电源','初始Uc/Il','等效电阻','电容/电感','零输入响应','零状态响应','全响应','退出');
switch StateFlag
case 1;
Us=input('请输入Us的值:');
Is=input('请输入Is的值:');
case 2
U0=input('请输入Uc的初始值:');
I0=input('请输入Il的初始值:');
case 3
Req=input('请输入Req的值:');
case 4
C=input('请输入C的值:');
L=input('请输入L的值:');
case 5 %零输入响应
figure('NumberTitle','off','Name','二阶电路零输入动态分析');
y=dsolve('L*C*D2y+Req*C*Dy+y=0','y(0)=U0','Dy(0)=-I0/C','t');
x=-C*diff(y,1);
z=L*diff(x,1);
Ul=eval(z)+I0/C;
Ic=eval(x);
Il=eval(x);
Uc=eval(y);
t=0:0.01:10;
plot(t,Uc,t,Ic,t,Ul,t,Il);
legend('Uc-t','Ic-t','Ul-t','Il-t');
%测试用
% Req=1;
% % L=1;
% % C=1;
% % I0=1;
% % U0=10;
case 6 %零状态响应
figure('NumberTitle','off','Name','二阶电路零状态动态分析');
high=10;
hold on
axis([0 8 -high,high]);
grid;
y=dsolve('L*C*D2y+1/Req*L*Dy+y=Is','y(0)=0','Dy(0)=-U0/L','t');%求Il;
x=L*diff(y,1);
z=-C*diff(x,1);
Ic=eval(z);
Uc=eval(x);
Ul=eval(x);
Il=eval(y);
t=0:0.01:10;
plot(t,Uc,t,Ic,t,Ul,t,Il);
legend('Uc-t','Ic-t','Ul-t','Il-t');
case 7 %全响应
hold on
high=Us+U0;
axis([0 8 -high,high]);
grid;
y1=dsolve('L*C*D2y+Req*C*Dy+y=0','y(0)=U0','Dy(0)=-I0/C','t');
x1=-C*diff(y,1);
z1=L*diff(x,1);
y2=dsolve('L*C*D2y+1/Req*L*Dy+y=Is','y(0)=0','Dy(0)=-U0/L','t');%求Il;
x2=L*diff(y,1);
z2=-C*diff(x,1);
Ic=eval(x1+z2);
Uc=eval(y1+x2);
Ul=eval(z1+x2);
Il=eval(x1+y2);
t=0:0.01:10;
plot(t,Uc,t,Ic,t,Ul,t,Il);
legend('Uc-t','Ic-t','Ul-t','Il-t');
case 8
StateFlag=0;
RoadFlag=0;
end
end
case 4
init();
case 5
RoadFlag=0;
end
end