一.求解方程组函数:solve

示例:求解含有参数的方程组  ax+by=10.ax-by=20

编程实现:

syms a b x y;
s=solve(a*x+b*y==10,a*x-b*y==20,x,y);
sol_x=s.x
sol_y=s.y

sol_x =
 
15/a
 
 
sol_y =
 
-5/b

二.求解微分方程:dsolve

1. 用字符串描述微分方程及其初始条件;
2. 导数的表示规则:以未知函数y为例,Dy”表示 y1阶导数,“D2y”表示对y2阶导数,其他各阶导数类似.

示例:

 

示例1:

syms y t;
ans=dsolve('Dy=(10-0.02*t)*t','y(0)=4',t)

ans =
 
4 - (t^2*(t - 750))/150



示例二:



三.求解微分方程的ode23函数

 

 

 

 

 示例如下:

 代码如下:

function qiujiefangcheng=myfun
y0=[0;2];
tspan=[0,3];%求解的范围
[t,y]=ode23(@fun,tspan,y0)




function  dfun=fun(t,y)
dfun=[-exp(1-t)*y(1)+0.8*y(2);
    y(1)-y(2).^3]

代码如下:

 

 

Logo

魔乐社区(Modelers.cn) 是一个中立、公益的人工智能社区,提供人工智能工具、模型、数据的托管、展示与应用协同服务,为人工智能开发及爱好者搭建开放的学习交流平台。社区通过理事会方式运作,由全产业链共同建设、共同运营、共同享有,推动国产AI生态繁荣发展。

更多推荐