最佳答案:lz需要给出y_0的初值~这个方程满足李普希兹条件,因此,解存在唯一并且可以唯一延拓到边界,应用lax等价定理可以知道,向前欧拉法具有二阶的收敛速度……所以可以
最佳答案:clearclc[t,x]=ode45(@(t,x)[-x(1)^2-x(2);x(1)-x(2)^3],[0 5],[1 0.5]);plot(t,x(:,1
最佳答案:配方法,使得ax〃+by〃=(a-b)x+(b-4a)y=k(ax+by),解出2a=b,k=-1所以x〃+2y〃=-(x+2y)可以解出x+2y的函数G(t)
最佳答案:最好贴全,谁调用了i_u_calculate,怎么更新的Ai值 查看原帖
最佳答案:挺难的,建议你发到学问社区上的mathematica小组,那高手挺多的
最佳答案:先用三角关系消去β,得到新方程组最后根据初值决定两个待定常数
最佳答案:按我的理解,由第二个表达式得到x=-[2cosy-sin(2t+y)]/[dy/dt-2]该式关于时间求导得到的dx/dt表达式代入第一个表达式,则得到关于y和
最佳答案:1.>>[x,y]=dsolve('D2x+2*Dx=x+2*y-exp(-t)','Dy=4*x+3*y+4*exp(-t)')或者>>a='D2x+2*Dx
最佳答案:知道dsolve函数就好求常微分方程或方程组了:)>> s=dsolve('Dv=-k*v-g*sin(a),Da=g*cos(a)/v','v(0)=v0,a
最佳答案:用ode15i试一试,我做过类似的,可以解决你的问题!
最佳答案:移项一下,将得到:dP1/dt=[c1*(p0-P1)-c2*(P1-P2)]/V1;dP2/dt=[c2*(p1-P2)-c3*P2]/V1;这个是常微分方程
最佳答案:matlab 返回的警告说的很清楚,你输入的微分方程组没有解析解.建议你尝试用ode45进行数值求解,最后查找q数组中接近0.8时所对应的h即为所求.