Plotting the output from differential equations

Question

onsagerian 2018-2-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/382510-plotting-the-output-from-differential-equations

评论： Walter Roberson 2018-2-14

Hello,

I am trying to design a program for making a plot using the outcome from differential equations. I have tried to construct a code several times for this, but I failed. Actually, documents on MathWorks for this issue are very confusing. Would you check my code designed to solve coupled differential equations and let me know what commands are needed to make plots? (Two plots should be generated. One for time vs S_1 the other for time vs S_2.)

I think the basic idea of the program is to save each outcome from given differential equations and use the data for plotting, but I don't know.

(Again, the following code is only for solving differential equations, so additional commands are required to plot data.)

syms S_1(t) S_2(t);
K=1;
k_1=20;
k_2=5;
k_3=5;
k_4=5;
k_5=2;
n=4;
dsolve(diff(S_1,t)==k_1*S_2^n/(K^n+S_2^n)-k_3*S_1-k_5*S_1,S_1(0)==0.0);
dsolve(diff(S_2,t)==k_2+k_5*S_1-k_4*S_2,S_2(0)==0.6);

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

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2018-2-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/382510-plotting-the-output-from-differential-equations#answer_305064

As I explained in https://www.mathworks.com/matlabcentral/answers/381943-solving-differential-equation-with-initial-conditions#comment_534054 those are related functions and must be dsolve() together. I also indicated there that the functions are effectively impossible to solve in closed form.

You will need to switch to numeric resolution. Star Strider posted links to odeFunction. The first example for odeFunction takes you through the entire workflow of functions needed to convert symbolic ode into inputs for the ode* routines.

2 个评论
显示无隐藏无

onsagerian 2018-2-14

编辑：Walter Roberson 2018-2-14

在 MATLAB Online 中打开

I am sorry, but I assumed once I construct differential equations and add "some commands" for plotting, the single program would automatically produce the plot I want to see. Of course, in general, once I get an output after solving differential equations, and use it to get plots in a separate file, then I would expect the result(plots). However, it is very inconvenient. So, is there any way to make a plot without solving differential equations? The following code is for solving a very simple differential equation with initial conditions. Again, without using a separate file, what commands are needed to get a plot from it?

syms y(x);
Dy=diff(y,x);
cond1=y(0)==2;
cond2=Dy(0)==5;
simplify(dsolve(diff(y,x,x)+4*diff(y,x)+4*y==(3+x)*exp(-2*x),cond1,cond2))

Walter Roberson 2018-2-14

In this simple example, dsolve is able to give you a closed form solution. You can store that closed form solution in a variable and fplot() that variable over the desired range.

In the case of the equations in your Question, above, there is no meaningful closed form solution and you will need to follow the work-flow show in odeFunction to translate your ode from symbolic to numeric and do a numeric simulation, which you can then plot.

请先登录，再进行评论。