| 查看: 3009 | 回复: 31 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
businiao412金虫 (正式写手)
|
[求助]
请教一个非线性方程非线性项转化为多项式的方法 已有2人参与
|
||
|
RT,楼主近日做振动问题时遇到一个二阶非线性微分方程如图片1，看文献中其有如下解法:在运用摄动法之前，先将其中的非线性项（1-q）3/2运用三阶最小二乘拟合的方法在区间1-Q<q<1上进行展开如图片2，及图片3，请问这是怎么做到的呢？我查了很多资料就是找不到这个方法的出处，本人数学知识较为薄弱，望不吝赐教。 图片1-3如下图所示：  方程1.png  2.png  3.png [ Last edited by businiao412 on 2014-8-10 at 15:16 ] |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有41人回复
-
回收溶剂求助
已经有7人回复
-
硝基苯如何除去
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
求助文献
已经有3人回复
-
三无产品还有机会吗
已经有6人回复
guojiashun
铁杆木虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 7191.2
- 散金: 82
- 红花: 12
- 帖子: 940
- 在线: 369.1小时
- 虫号: 933003
- 注册: 2009-12-25
- 性别: GG
- 专业: 机械动力学
|
解方程组的方法应该和点取多少没关系，四个未知量，凑够四个方程就行。我试了一下：取x=[1-Q,1-6*Q/10,1-3*Q/10,1]，y=[-Q^(3/2),-(6*Q/10)^(3/2),-(3*Q/10)^(3/2),0]， X=[x(1)^3,x(1)^2,x(1),1; x(2)^3,x(2)^2,x(2),1; x(3)^3,x(3)^2,x(3),1; x(4)^3,x(4)^2,x(4),1;] a=inv(X)*y' =25/7/Q^3*conj(Q^(3/2))-5/3/Q^3*3^(1/2)*5^(1/2)*conj(Q^(3/2))+10/21/Q^3*3^(1/2)*10^(1/2)*conj(Q^(3/2)) 15/14/Q^3*(-10+3*Q)*conj(Q^(3/2))-1/6/Q^3*(-30+13*Q)*3^(1/2)*5^(1/2)*conj(Q^(3/2))+2/21/Q^3*(-15+8*Q)*3^(1/2)*10^(1/2)*conj(Q^(3/2)) 3/14/Q^3*(50-30*Q+3*Q^2)*conj(Q^(3/2))-1/6/Q^3*(30-26*Q+3*Q^2)*3^(1/2)*5^(1/2)*conj(Q^(3/2))+2/21/Q^3*(15-16*Q+3*Q^2)*3^(1/2)*10^(1/2)*conj(Q^(3/2)) -1/14/Q^3*(50-45*Q+9*Q^2)*conj(Q^(3/2))+1/6/Q^3*(10-13*Q+3*Q^2)*3^(1/2)*5^(1/2)*conj(Q^(3/2))-2/21/Q^3*(5-8*Q+3*Q^2)*3^(1/2)*10^(1/2)*conj(Q^(3/2)) 根据论文知道Q必然大于0，所以conj(Q^(3/2))=Q^(3/2)，这样整理符号运算所得第一项：a3=(25/7-5*sqrt(15)/3+10*sqrt(30)/21)/Q^(3/2)=-0.2753/Q^(3/2)≈ -8/(33*Q^(3/2))，其它不再逐一整理。符号运算得到的a3表达式和论文中还是有误差，这应该是matlab矩阵求逆算法的误差造成的。四元一次方程组应该也不难解，如果手工解是不是应该完全精准？ |
19楼2014-08-14 14:27:18
herowolf
木虫 (正式写手)
- 应助: 25 (小学生)
- 金币: 3766.3
- 红花: 8
- 帖子: 731
- 在线: 381.1小时
- 虫号: 3283264
- 注册: 2014-06-19
- 性别: GG
- 专业: 流体力学
2楼2014-08-10 16:53:53
guojiashun
铁杆木虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 7191.2
- 散金: 82
- 红花: 12
- 帖子: 940
- 在线: 369.1小时
- 虫号: 933003
- 注册: 2009-12-25
- 性别: GG
- 专业: 机械动力学
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
businiao412: 金币+15, ★★★★★最佳答案, 衷心感谢您细致透彻的讲解，有种豁然开朗的感觉，我按照您的思路试做一下，有问题再向您请教。 2014-08-11 14:57:17
感谢参与，应助指数 +1
businiao412: 金币+15, ★★★★★最佳答案, 衷心感谢您细致透彻的讲解，有种豁然开朗的感觉，我按照您的思路试做一下，有问题再向您请教。 2014-08-11 14:57:17
|
对你的问题理解不是十分准确。仅仅就我的理解说一下： 摄动法的基本思想是把响应q看做多个不同量级的摄动的叠加，即：q(t,ε )=u0(t)+ε*u1(t)+ε ^2*u2(t)+......【1】（多尺度法还要把考虑不同的时间尺度，q(t,ε)=u0(t1,t2,...)+ε*u1(t1,t2,...)+ε ^2*u2(t1,t2,...)+.....），然后把q的近似表达式【1】代入原方程，通过整理和推导，根据方程两边 小参数 ε 的同次幂系数相等的条件，把原方程转化为分别关于u0、u1、u2.....的微分方程，然后分别逐次求得各摄动u0、u1、u2........最后，叠加可得方程q的近似解析解。 这个非线性微分方程的非线性项是（1-q）^3/2，而在运用摄动法把式【1】代入原方程中进行整理时，由于（1-q）^3/2这一项带有分数次幂，故【1】代入后其展开是非常困难的，即使展开其中也带有 ε 的分数次幂项，也就无法利用 ε 的同次幂系数相等的条件来求解（后面的套路步骤走不下去了）。所以，一般在求解此类非线性微分方程时，用多项式来近似表达诸如（1-q）^3/2这种难于处理的项，这样当式【1】代入后展开、整理就变得很简便了。类似的情况如齿轮动力学理论中，一个主要的非线性因素是齿侧间隙，间隙作用方式用分段函数来表达，而分段函数的存在也给摄动法、多尺度法等的运用带来困难，数学上难于展开处理，故也有学者用多项式来近似这个分段函数（当然还有其它处理方法，这个还要视造成困难的项的具体函数性质）。 这里提前预估了一个拟合区间q∈[1-Q, 1]，Why？第一，方程首先要有稳态解；第二，也是为了数学处理上的便利性。（1）根据（1-q）^3/2可知，q必须小于等于1，否则出现复数项，方程无稳态解，故已知q的上限；（2）采用多项式来拟合（1-q）^3/2，取多少阶合适？我们当然是希望阶数越低越好。阶数越高，式【1】代入后的数学展开越复杂繁琐。根据数学知识易知拟合区间越小，就可以用越低的阶数来获得足够的拟合精度，所以根据已知条件提前划定一个尽量小的区间十分有益和必要。接下来要确定q的下限，其中Q是q所能预期到的最大的峰值（上限）和谷值（下限）之间的差值（peak-to-peak value，我想这里的Q可以通过数值求解、实验或者实际观测得到），那么已知q的上限是1，下限自然就是1-Q了。 最后，文中确定用三次多项式来拟合，接着就要确定多项式的系数a0、a1、a2、a3。这里采用“最小二乘法”来确定，最小二乘法是一种数值方法，必须基于已知的数据进行拟合。上面已经得到了拟合区间q∈[1-Q, 1]，那么在这个区间内取一些点，如1-Q, 1-9*Q/10, 1-8*Q/10, 1-7*Q/10, ......, 1，依次代入（1-q）^3/2中，可得相应值Q^3/2, (9*Q/10)^3/2, (8*Q/10)^3/2, (7*Q/10)^3/2, ......, 0，已知了函数（1-q）^3/2的一组自变量和因变量，套用“最小二乘法”的套路解得多项式的系数a0、a1、a2、a3也就不难了吧？ 总结，不论是摄动法，还是多项式拟合，都是无限逼近思想的体现。摄动法本身得到的就是方程的无限近似解，因此用多项式无限近似方程中复杂项也就变得无可厚非了。 PS：周日下午无事，多敲了些废话，纯属个人见解，限于水平，有误之处请各位指正。 |
3楼2014-08-10 19:43:10
businiao412
金虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 1431.3
- 红花: 9
- 帖子: 469
- 在线: 58.9小时
- 虫号: 1496676
- 注册: 2011-11-17
- 专业: 石油天然气开采
4楼2014-08-11 14:59:13