问题标题:
怎么用ode23求解一阶方程?看了一本很基础的书,就是不知道怎么输入表达式才行书是这么说的:“dy/dt=cos(t)我们很容易求得方程的积分为y(t)=sin(t)+C,它将可以用来检验我们得到的数值解.
问题描述:
怎么用ode23求解一阶方程?
看了一本很基础的书,就是不知道怎么输入表达式才行
书是这么说的:
“dy/dt=cos(t)
我们很容易求得方程的积分为y(t)=sin(t)+C,它将可以用来检验我们得到的数值解.
首先我们定义函数,我们创建一个.m文件,输入下面的内容.
functionydot=eq1(t,y)
ydot=cos(t);
我们的函数称为eq1.让我们求0≤t≤2π内的解并假设y(0)=2.调用的语句是:
>>[t,y]=ode23('eq1',[02*pi],2);”
我按说着的步骤File-New-M-File,然后输入
functionydot=eq1(t,y)
ydot=cos(t);
接着保存到安装目录下的work文件夹下,左侧的currentdirectory也切换到work文件夹了,然后就是在commandwindow下输入:
[t,y]=ode23('eq1',[02*pi],2);
提示没有定义eq1,直接输入eq1文件名,又提示Inputargument"t"isundefined.
请问如何输入命令?
康建斌回答:
大哥,这个是matlab程序,你怎么发到C++这里了?不过我还是可以帮你看看:functionydot=eq1(t,y)ydot=cos(t);这个文件要保存为:eq1.m接着保存到安装目录下的work文件夹下:不要这么做,只需要放在你当前运行的文件...
点击显示
其它推荐