matlab 解二阶微分方程 >>
我解的方程是光在非均匀介质中的传播,介质折射率为位置的函数n=2-(x^2+y^2)
用费马定理得到二阶微分方程y''+2(y+xy')(1+y'^2)/(2-x^2-y^2)=0
我把它转换为两个一阶微分方程,t=x ,y(1)=y ,y(2)=y' ,
所以y(2)'=2(y(1)+ty(2))(1+y(2)^2)/((t^2+y(1)^2)-2)]
其中边界条件为x=-1,y=0
从上面的思路,
需要解的就是这个二阶微分方程,下面是matlab语句,报错
>> clear;close;
>> sinit=bvpinit(-1:1,[1;0]);
>> odefun=inline('[y(2);2*(y(1)+t*y(2))*(1+y(2)^2)/((t^2+y(1)^2)-2)]','t','y');
>> bcfun=inline('[ya(1)+1;yb(1)]','ya','yb');
>> sol=bvp4c(odefun,bcfun,sinit)
运行结果报错:
Error using ==> bvp4c at 203
Unable to solve the collocation equations -- a singular Jacobian encountered
这个二阶微分方程该怎么解呢,
请给正确的语句,
人气:127 ℃ 时间:2020-04-15 05:32:47
解答
先把方程化为一阶的,然后利用ode45等函数进行求解,自己试一下,如果还有疑问再说.
推荐
猜你喜欢
- 同义句转换
- 初一语文上册第一单元作文
- 填开字词语:()伟业,()荒地,()产品,()思想,()矿藏,()道路
- 9X+270=90乘45
- 神州七号 翟志刚太空漫步问题
- 甲乙两车从A丶B两地相向而行,己知甲车的速度是乙车的5\6.甲车先从A城开出55千米后,乙车才从B城出发.两车
- 在CAD中我画图时按1:1的比例画图,当我打印出图怎么才知道图纸与实际的比例是1:多少?
- 有一电动机,额定功率为300kw,接在1140v的交流电源上,通过电动机的额定电流为200A,试求电动机的功率因数.