matlab二阶微分方程已知函Y=(1.082X)^1/2 用四次牛顿前插公式计算f(-207) 求matlab源程序 最好能列出差分表
matlab二阶微分方程已知函Y=(1.082X)^1/2 用四次牛顿前插公式计算f(-207) 求matlab源程序 最好能列出差分表
matlab二阶微分方程
已知函Y=(1.082X)^1/2 用四次牛顿前插公式计算f(-207) 求matlab源程序 最好能列出差分表
matlab二阶微分方程已知函Y=(1.082X)^1/2 用四次牛顿前插公式计算f(-207) 求matlab源程序 最好能列出差分表
>> clear
>> syms a b c d e;
>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e')
y =
(2*a*e + b*d + d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b - (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2)) - (2*a*e + b*d - d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b + (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2))
>> a=1;b=1;c=1;d=1;e=1;%若常数已知
>> t=1;%desolve中没指定自变量x,这里默认为t,带入任意一个x值
>> y=eval(y)
y =
1.1932 + 0.0000i
这是解方程的方法,画图也差不多,我在命令窗写的改麻烦,下面是我在editor写的.
syms a b c d e;
y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e','x');
x=200:0.1:400;%取值范围,步长
a=1;b=1;c=1;d=1;e=1;%带入参数
y=eval(y);
plot(x,y)
你补充问题的方程解得也是一个空的,并报错,或许就是无解的,我也很疑问.
打字不易,