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.乘除法
.*(点乘),矩阵每一项对应相乘,同./ .^