早教吧作业答案频道 -->其他-->
求大神支招,我想用ode45来求解一个微分方程的数值解,但是MATLAB老是说有错,怎么办?程序是clear all;options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);[T,Y]=ode45(@guanxi,[0 12],[0 1 1],options);对应 guanxi.m 是funct
题目详情
求大神支招,我想用ode45来求解一个微分方程的数值解,但是MATLAB老是说有错,怎么办?
程序是
clear all;
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);
[T,Y]=ode45(@guanxi,[0 12],[0 1 1],options);
对应 guanxi.m 是
function dy=guanxi(t,y)
dy = zeros(2,1);
Dy(1)=Dy(2)*4*y(1)/3*((-2)*16*exp(-0.5*t)+85.5-y(2)) - Dy(2)*4*(-0.5)*16*exp(-0.5*t)/3*((-32)*exp(-0.5*t)+85.5-y(2));
自己认为可能的问题所在市我的微分式中有指数函数这个东西,所以不行~只是感觉~
程序是
clear all;
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4]);
[T,Y]=ode45(@guanxi,[0 12],[0 1 1],options);
对应 guanxi.m 是
function dy=guanxi(t,y)
dy = zeros(2,1);
Dy(1)=Dy(2)*4*y(1)/3*((-2)*16*exp(-0.5*t)+85.5-y(2)) - Dy(2)*4*(-0.5)*16*exp(-0.5*t)/3*((-32)*exp(-0.5*t)+85.5-y(2));
自己认为可能的问题所在市我的微分式中有指数函数这个东西,所以不行~只是感觉~
▼优质解答
答案和解析
初值不够呀,应该有四个,题目能不能贴出来?
看了 求大神支招,我想用ode45...的网友还看了以下:
已知y=f(x)是偶函数,且在[0,+∞)是减函数,求函数f(1-x²)的单调区间为什么我算的不对 2020-04-25 …
一道函数填空题求解!已知函数y=mx²+(m-3)x+1的值域是[0,正无穷),则实数m的取值范围 2020-05-14 …
[0,1]内的所有数字的平均数是多少?我是这样想的,假设[0,1]内有n个数,那么根据连续求和的那 2020-05-22 …
0到底是不是分数呢?我个人认为0是一个分数因为0/任何号码都可以.可以给我一个更清楚地接释吗?如果 2020-05-24 …
我们把0.00000000025简单记作0.00…025小数点后有9个0.下面有两个小数,a=0. 2020-06-17 …
0.5^x=x我想了好久!指数函数y=0.5^x和对数函数y=log0.5x的交点怎么求?就是方程 2020-07-20 …
非正整数是0吗他们说0是非正整数可是0.1是正数但他不是正整数而且0是整数啊我好矛盾啊0是非正整数可 2020-11-03 …
请问“投出考生数”是什么意思?如果一个学校在我们省计划招2个人投出考生数是“0”请问这是什么意思 2020-11-06 …
求找一名言那个名言是你有?我有?意思是你有对付我的招数我也有对付你的招数请补充问号里的是什么! 2020-11-17 …
我们老师说:不论什么时候,a和-a都互为相反数,a=0-a=-0-0有意义吗?0不是正数和负数,为什 2021-02-21 …