早教吧作业答案频道 -->其他-->
求大神支招,我想用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...的网友还看了以下:
求极限lim[cosx-e^(-x^2/2)]/x^4 其中x趋向于0.我的做法为什么错了:lim 2020-05-16 …
EXCLE表格中有四列需要做平均值运算,但是特殊情况为有些列的值为0,求救公式一条0不参与运算,如 2020-05-16 …
数据结构关于图的问题已知一个带权图的顶点集V和边集G分别为:V={0,1,2,3,4,5};E={ 2020-06-17 …
A,B均为三阶可逆矩阵,且A^3=0,则A:E-A,E+A均不可逆?B:E-A不可逆但E+A可逆? 2020-07-20 …
已知点A(0,0)B(2,3)C(2,4)D(5,5)E(1,4)F(0,6)连接ABCDEF将线 2020-08-02 …
积分的问题∫∫e^(x^2+y^2)∫∫e^(x^2+y^2)=∫(0,π/2)dθ∫(0,R)re 2020-10-31 …
求证e^i(4π/n)+e^i(8π/n)+...+e^i4(n-1)π/n+e^i(4nπ/n)= 2020-11-01 …
已知向量a≠e,|e|=1,对任意t∈R,恒有|a-te|≥|a-e|,则(a,e皆为向量)A.a⊥ 2020-11-02 …
n阶方阵A满足A^2=O,E是n阶单位阵,则A.|E-A|≠0,但|E+A|=0B|E-An阶方阵A 2020-11-02 …
母亲为红绿色盲,父亲正常,其四个儿子有()可能是色盲.A.1个B.2个C.3个D.4个E.0个 2020-11-25 …