Effects of parameter on oscillatory behaviour.

Question

James Marlom 2018-3-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/388563-effects-of-parameter-on-oscillatory-behaviour

评论： Star Strider 2018-3-16

How can I keep changing only Kc value (1 to 5) and using it in the Differential Equation.

if true
  clc; clear all;
 A=2; %Cross Sectional Area of the Tank [m^2]
 Kc=1; %Proportional gain [m^2/min]
 t1=0.1; %Integral Time Constant [min]
 syms t y(t)
 ode= A*diff(y,t,2)+Kc*diff(y,t,1)+(Kc/t1)*y==0;
 Dy=diff(y);
 cond1 = y(0) == 0;
 cond2 = Dy(1e-10) == 2;
 conds=[cond1 cond2];
 sol=dsolve(ode,conds,t)
 ezplot(sol,[0,20])
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2018-3-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/388563-effects-of-parameter-on-oscillatory-behaviour#answer_310251

在 MATLAB Online 中打开

Leave ‘Kc’ as a symbolic variable. Then you can vary it in your fsurf call:

A=2; %Cross Sectional Area of the Tank [m^2]
%  Kc=1; %Proportional gain [m^2/min]
t1=0.1; %Integral Time Constant [min]
syms Kc t y(t)
ode= A*diff(y,t,2)+Kc*diff(y,t,1)+(Kc/t1)*y==0;
Dy=diff(y);
cond1 = y(0) == 0;
cond2 = Dy(1e-10) == 2;
conds=[cond1 cond2];
sol(t,Kc)= dsolve(ode,conds,t)
fsurf(sol,[0,20,1,5])
xlabel('t')
ylabel('Kc')
view(30,30)