早教吧作业答案频道 -->其他-->
求大神支招,我想用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...的网友还看了以下:
v=s/t是用单位时间内通过的路程来表示运动快慢的。能不能用单位路程所用的时间来表示运动的快慢? 2020-05-14 …
怎么用MATLAB解带参数的一元三次方程解方程:(-8)*(x^3)+9*(x^2)-16*((t 2020-05-14 …
椭圆方程式的题椭圆方程x^2/a^2+y^2/b^2=t怎么用y=f(u),x=f(u)表示如果椭 2020-05-16 …
一道比较简单的一阶线性方程应用题题目如下(英文的),只考虑第一小题即可.(1)200.32feet 2020-06-06 …
用方程来用华氏温度99是摄氏度的多少 2020-06-10 …
spss16.0中的近似t检验是怎么用的啊,有没有教程或者图示,可以直观点是这样的,我有两组数据, 2020-06-13 …
怎么计算二阶微分方程(mathematrica)我想用mathematica5来书写,但是我又不知 2020-07-31 …
一根铁棒在40℃时长12.506m,在80℃时长12.512m,已知长度l(cm)和温度t(℃)的 2020-08-01 …
填t填.(1)测量角t大小用;度量角t单位是,用符号“”表示,把半圆平均分成180份.每t份所对t 2020-08-02 …
利用变换x=lnt将微分方程d2y/dx^2-dy/dx+e^(2x)y=0化为关于t的微风方程.利 2020-12-12 …