早教吧作业答案频道 -->其他-->
求大神支招,我想用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...的网友还看了以下:
求大神支招,我想用ode45来求解一个微分方程的数值解,但是MATLAB老是说有错,怎么办?程序是 2020-05-16 …
已知函数f(x)=kx+k,x≤0lnx,x>0(其中k≥0),若函数y=f[f(x)]+1有4个 2020-06-12 …
在fluent中,残差曲线纵坐标上的1.12e-03,或是1e+05之类的都是什么数呀,是不是表示 2020-06-21 …
M是一种治疗艾滋病的新药(结构简式见图),下列关于M的说法正确的是()A.M易溶于水,且可以发生水 2020-07-05 …
下列对应:①M=RN=N*对应关系f:对集合M中的元素,取绝对值与N中元素对应②M={1,2,-1 2020-07-09 …
(2014•湖南二模)设x=a和x=b是函数f(x)=lnx+12x2-(m+2)x的两个极值点, 2020-08-02 …
(2014•上海模拟)定义一个对应法则f:P/(m,n)→P(m,n),(m≥0,n≥0).现有点 2020-08-02 …
定义一个对应法则f:P(m,n)→P′(m,n),(m≥0,n≥0).现有点A(2,6)与点B(6 2020-08-02 …
定义一个对应法则f:P(m,n)→P′(m,n),(m≥0,n≥0).现有点A(1,3)与点B(3 2020-08-02 …
∫区间0到X的e的-u次du我感觉应该是求e的-u次的原函数也就是e的-u次,然后属于复合函数在求- 2021-02-16 …