用matlab求解微分方程

matlab使用笔记
1.Excel读数据：
x = xlsread('路径') 读数据
xlswrite('路径'，变量) 将变量的值保存进Excel
路径的格式是cmd中cd命令跟的路径

Ctrl+R多行注释
Ctrl+T取消多行注释

2.求解一阶线性常微分方程
clear all
clc

syms y(x) %符号声明y是x的函数

ode=diff(y,x)-y==x;%diff（y，x）指的是y对x求一阶导，ode指的是常微分方程，常微分方程的等号用双等号

cond=y(0)==1;%cond指的是初始条件

dsolve(ode,cond)%dsolve（方程，初始条件）指的是求解常微分方程的函数

%matlab中赋值需要用分号隔开，运行函数不用

3.一阶非线性
同上
幂运算用^这个符号，例如y^2
dsolve的返回值就是解的函数，例如ys（x）=dsolve（ode，cond），用ode求出来的y（x）就是给ys（x）赋的值

4.二阶微分方程
diff（y,x,2）y对x求二阶导
初始条件有两个，写在方括号里向量的形式
dsolve(ode,cond,'implicit',true)解方程的隐式解
如果显式解和隐式解都无法得到解，则只能用数值解法

5.三阶微分方程
dsolve函数会运用简化规则，例如arcsin（sinx）=x，于是nπ周期性解会被省略掉，严格意义上会漏解，设置‘IgnoreAnalyticConstraints'为false则不应用简化规则（用得不多）

带未知参数运算要在syms后面声明是符号变量，例如syms a b c d e，然后这五个字母就可以在ode中参与运算了
simplify（ys(0)）可以简化在0处的运算值

ode求解同上

6.解微分方程组
ode=[ode1，ode2]
dsolve返回的值是多个函数，以向量的形式，例如[xs(t),ys(t)]=dsolve(ode,cond)
也可以sol=dsolve(ode,cond)，访问的时候sol.x sol.y

7.没有解析解，采用数值解的办法，运用ode45库
数值解不会返回y与x的解析关系，会返回x和y一一对应的一些值
只能求解一阶微分方程，且一阶项不是平方

一阶
[x,y]=ode45(@ode,tspan,y0,opts)，ode45是求解器，返回x，y的一一对应的数组，@ode是微分方程，tspan是积分区间，y0是初始值，opts是选项，用来规定相对容差和绝对容差，可以帮助提高计算速度或提高计算精度，例如opts=odeset('reltol',1e-2,'abstol',1e-4)

tspan=[0 5]，定义积分区间从0到5，中间是空格
y0=1;，定义初始值
opts=[];，定义选项

function dy=ode(x,y)
dy=x+y
end，用局部函数来定义ode，有格式要求：dy（y的一阶导）在左，其他在右

最后要把ode45这行命令放在最下面

不是一阶
首先要化成一阶的微分方程组，可以把y拆分成两个变量y(1)和y(2)，即拆成两个一阶方程
例如求解y''-(1-y^2)y'+2y=0，令y1=y，y2=dy1，方程变为y2’-（1-y1^2)y2+2y1=0；y2=dy1；两个方程组成的方程组

tspan=[0 5];
y0=[2 0]  %0处y1，y2的值

[x,y]=ode45(@ode,tspan,y0,opts)  %此时y=[y1 y2]

function dy=ode(x,y)
dy=zeros(2,1)  %声明dy是两行一列的0矩阵
dy(1)=y(2)
dy(2)=(1-y(1)^2)*y(2)-2*y(1)  %转化成的方程组
end

8.用数值解，求特定x的函数值
sol=ode45(@ode,tspan,yo,opts)，得到的是解得结构体
x=0.5
y3=deval(sol,x)，即可得到x=0.5时候的函数值

如果x在积分区间外面，用扩展函数，或者直接改积分区间
sol2=odextend(sol,@ode,8)，将积分区间扩展到8
x=6
y4=deval(sol2,x)

9.其他ode求解器及其特点
ode45，综合性能相对较好
ode15s，对于刚性问题求解能力较好，比如说图像有棱有角的这种趋于多段线性的这种
ode23，求解速度相对较快，精度相对较低
ode113，求解速度相对较慢，精度相对较高

10.隐形方程用数值法求解
tspan=[0 5];
x0=0;
y0=[2 0];
yp0=[0 0];
[y0,yp0]=decic(@ode,x0,y0,[2 0],yp0,[0 0])
opts=[];
sol=ode15i(@ode,tspan,y0,yp0,opts)
function f=ode(x,y,yp)
f=zeros(2,1);
f(1)=yp(1)-y(2);
f(2)=yp(2)-(1-y(1)^2)*y(2)+2*y(1)
end

11.乘除法
.*（点乘），矩阵每一项对应相乘，同./  .^