matlab二阶微分方程
已知函Y=(1.082X)^1/2用四次牛顿前插公式计算f(-207)求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)
你补充问题的方程解得也是一个空的,并报错,或许就是无解的,我也很疑问。
打字不易,如满意,望采纳。
>> 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)
你补充问题的方程解得也是一个空的,并报错,或许就是无解的,我也很疑问。
打字不易,如满意,望采纳。
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询