美赛建模需要matlab吗,美赛(matlab自学)之微分方程建模

人口模型

malthus模型

美国人口的预报模型

建模与求解

参数估计

(1)非线性最小二乘估计

clc, clear

a = textread('data4.txt'); %把原始数据保存在纯文本文件data4 .txt中

x=a([2:2:6],:)';%提出人口数据

x = nonzeros(x);%去掉后面的0,并变成列向量

t =[1790:10 :2000]';

t0 =t(1); x0 =x(1);

fun=@ (cs,td)cs(1)./1 +(cs(1)/x0 -1) *exp( -cs(2) * (td-t0)));%cs(1) =xm,cs(2) =r

cs = lsqcurvefit( fun, rand(2,1) ,t(2:end) ,x(2 :end) , zeros(2,1))

xhat = fun(cs,[t;2010])%预测已知年代和2010年的人口

(2)线性最小二乘法

clc, clear

a = textread( 'data4.txt');%把原始数据保存在纯文本文件data4.txt中

x=a([2:2:6],:)'; x = nonzeros(x);

t =[1790:10 :2000]';

a=[ones(21,1)，-x(2:end)];

b=diff(x) ./x(2 :end) /10;

cs =a\b;

r =cs(1), xm=r/cs(2)

matlab求微分方程

微分方程通解

syms y(x)

dsolve(x^2+y+(x-2*y)*diff(y)==0)

运行结果:

微分方程的初边值问题

求解微分方程组

求解方程组：

clc, clear

syms f(x) g(x)%定义符号变量

df =diff(f);%定义f的一阶导数,用于初值或边值条件的赋值

[ f1,g1] =dsolve(diff(f,2) +3 *g= =sin(x),diff(g) +df= =cos(x))%求通解

f1 =simplify(f1), g1 =simplify(g1)%对符号解进行化简

[f2 ,g2] =dsolve(diff(f,2) +3 *g= =sin(x),diff(g) +df= =cos(x),df(2)= =0,f(3) = =3,g(5) = =1)

f2 =simplify(f2), g2 =simplify(g2)%对符号解进行化简

求解线性微分方程组

齐次线性微分方程组

clc,clear

syms x(t) y(t) z(t)

X=[x;y;z];

A=[2 1 3;0 2 -1;0 0 2];

B=[1 2 1]';

[x,y,z]=dsolve(diff(X)==A*X,X(0)==B)

非齐次线性微分方程组

clc,clear

syms x(t) y(t) z(t) %定义符号变量

X=[x;y;z];

A=[1,0,0;2,1,-2;3,2,1];B=[0;0;exp(t) *cos(2*t)];

X0=[0;1;1];%初值条件

X =dsolve(diff(X)==A*X +B,X(0)==X0)%求符号解

X=simplify([X.x;X.y;X.z])%显示解的各个分量

pretty(x)%分数线居中的显示方式