| 查看: 3186 | 回复: 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 ] |
回复此楼» 猜你喜欢
一志愿大连理工大学材料求调剂
已经有5人回复
-
材料与化工调剂一志愿大连海事085600，349
已经有6人回复
-
一志愿中海洋320化学工程与技术学硕求调剂
已经有4人回复
-
材料科学与工程调剂
已经有11人回复
-
311求调剂一志愿合肥工业大学
已经有4人回复
-
一志愿食品科学与工程083200求调剂
已经有4人回复
-
材料科学与工程求调剂
已经有10人回复
-
一志愿北化085600材料专硕275|有文章专利｜求调剂
已经有11人回复
-
一志愿 南京航空航天大学 ，080500材料科学与工程学硕
已经有4人回复
-
085600材料与化工调剂
已经有15人回复
guojiashun
铁杆木虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 7191.2
- 散金: 82
- 红花: 12
- 帖子: 940
- 在线: 369.1小时
- 虫号: 933003
- 注册: 2009-12-25
- 性别: GG
- 专业: 机械动力学
|
我用最小二乘法做了一下： syms Q x=[1-Q:0.001*Q:1]; y=-(1.-x).^(3/2); N=length(x); syms A for n=1:N for m=1:4 A(m,n)=x(n)^(m-1); end end a=(A*A')\(A*y'); 其中a(1)=a0; a(2)=a1; a(3)=a2; a(4)=a3;然后得到很长很长很长一串a关于Q的表达式，人工化简基本不可能了.............然后用eval函数化简最简单的 a3 项，得 eval(a(4))=-1230776300378990335736947704895/5070602400912917605986812821504/conj(Q)^3*conj(Q^(3/2))= -0.2427/Q^(3/2)= -8/(33*Q^(3/2)) 但是其它几项太复杂了，化简不动......... |
29楼2014-08-14 19:03:52
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
