大师作文网MATLAB 解微分方程
来源：学生作业帮 编辑：大师作文网作业帮 分类：数学作业 时间：2024/11/10 11:10:44
MATLAB 解微分方程
V'=qi-q0;
m'=qi*pi-q0*p0;
m=55 65 69 72 78 85
qi=15 14 18 19 21 20
pi=2 1 2 4 3 4
v=20 22 25 27 31 33
数据采样时间间隔1S
求解p0 q0 的值
V'=qi-q0;
m'=qi*pi-q0*p0;
m=55 65 69 72 78 85
qi=15 14 18 19 21 20
pi=2 1 2 4 3 4
v=20 22 25 27 31 33
数据采样时间间隔1S
求解p0 q0 的值
这样可以吗?
t=1:6
m=[55 65 69 72 78 85]
qi=[15 14 18 19 21 20]
ppi=[2 1 2 4 3 4]
v=[20 22 25 27 31 33]
plot(t,[m;v],'o-')
vsp=spline(t,v)
dv=fnval(fnder(vsp),t)
msp=spline(t,m)
dm=fnval(fnder(msp),t)
q0=qi-dv
p0=solve('dm=qi*ppi-q0*p0','p0')
p0=subs(p0)
结果：
q0 =
14.7667 10.8667 15.7667 16.0667 16.9667 21.0667
p0 =
1.0414 0.7188 2.1170 4.4613 3.2954 3.4974
t=1:6
m=[55 65 69 72 78 85]
qi=[15 14 18 19 21 20]
ppi=[2 1 2 4 3 4]
v=[20 22 25 27 31 33]
plot(t,[m;v],'o-')
vsp=spline(t,v)
dv=fnval(fnder(vsp),t)
msp=spline(t,m)
dm=fnval(fnder(msp),t)
q0=qi-dv
p0=solve('dm=qi*ppi-q0*p0','p0')
p0=subs(p0)
结果：
q0 =
14.7667 10.8667 15.7667 16.0667 16.9667 21.0667
p0 =
1.0414 0.7188 2.1170 4.4613 3.2954 3.4974