matlab微分方程求解总结
[Matlab][隐式常微分方程组]ODE15i中odefun中的y应该如何输入？
[Matlab][隐式常微分方程组]ODE15i中odefun中的y应该如何输入？
根据ode15i文档中所说，odefun的格式必须为
作为参数odefun的应该是Function Handles ，即函数句柄。例子如下：
fhandle @(t,y,Dy)(t*y^2*(Dy)^3-y^3*(Dy)^2 t*(t^2 1)*Dy-t^2*y)
fhandle
@(t,y,Dy)(t*y^2*(Dy)^3-y^3*(Dy)^2 t*(t^2 1)*Dy-t^2*y)
[t,y] ode15i(fhandle,[1 10],sqrt(3/2),0);
ytrue sqrt(t.^2 0.5);
plot(t,y,t,ytrue,o);
主要问题在于fhandle，在这个例子中fhandle是以@符号作为标示符的匿名函数的句柄。匿名函数的写法普通函数类似，不过函数主体放在@后第二个小括号中，@后第一个小括号中列出用到的变量。
总之，按照示例中的格式来就可以了。
运行结果：
MATLAB常微分方程如何写一个二维数组？
如果你只是想统计某个数 比如10 出现多少次的话可以这样 sum( a(:) 10) 如果你想得到一个列表说明每个元素各出现了多少次的话可以这样 r hist( a(:), unique(a)) 结果是按unique(a)里面的顺序给出的每个元素的出现次数
您好!请问如何在Matlab中用四阶龙格库塔法解微分方程？
四阶龙格-库塔是求解常微分方程（常微分方程组）精度最高的一种数值方法。四阶龙格-库塔迭代公式为根据四阶龙格-库塔迭代公式，可用matlab软件编写得到如下计算代码：
调用方法：clcclose ally01h0.1a0b1[x,y] runge_kutta(@func,y0,h,a,b)disp(#39 x y#39)[x#39,y#39]x0:0.1:1 plot(x,y,#39:#39),grid onxlabel(#39x#39),ylabel(#39y（x）/y’(x)#39)title(#39函数图#39)
function f func(x,y)end