二维非线性方程的处理方法

zjbwow

大家好，最近在搞一个二维非线性方程组的处理方法：
a1*p1+a2*p1*y1+a3*(y1)'+a4*y1+a5*y2+a6*y1^2+a7*y1*y2+a8*y1*(y1)''=0
b1*p1+b2*p1*y2+b3*(y2)'+b4*y2+b5*y1+b6*y2^2+b7*y1*y2+b8*y2*(y2)''=0
其中(y1)'代表一阶导数，(y1)''代表二阶导数
这个系统有一个输入p1，两个输出y1和y2，并且a、b系数相差很大，最大值a6是最小值a8的1e7倍。
我想先转化为模态坐标对方程解耦，再使用非线性的方法进行分析。
不知道思路对不对，还望大牛学霸们多多指教！