SOS:这样的微分方程组该怎么处理?
需要解一个方程组,模拟出A(x)随x的变化曲线:E(x)^2-2bA(x)A'(x)E(x)-ag^2=0 (1)
A''(x)+ cE(x)^2 A(x)-dA(x)= 0 (2)
其中A(x)是待求的光波振幅的空间分布函数;A'(x)表示A(x)对x坐标的一阶导;A''(x)表示A(x)对x坐标的二阶导;E(x)是电场的分布函数a,b,c,d,g都是个给定的常数。
从数学上分析,方程(1)就是一个关于E(x)的一元二次方程[ A(x)看作常系数],因为一元二次方程有两个根,所以从方程(1)中应该能解出两个E(x):E1(x)与E2(x),即得出两个E(x)与A(x)的对应关系。
并且可以看出E1(x)恒大于0,E2(x)恒小于0。而从物理上考虑,E(x)应该有时大于0,有时小于0(类似于正弦曲线)。也就是说在某个区间取E1(x),在另一个区间则会取E2(x)。
那么这个程序能编写吗?怎么处理呢?求高人指点,不胜感激! 可以编写,这是一种微分代数方程!
详细的,可以看一下薛定宇《高等应用数学问题的MATLAB求解》P231
页:
[1]