| 查看: 3174 | 回复: 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 ] |
回复此楼» 猜你喜欢
285求调剂
已经有5人回复
-
289求调剂
已经有7人回复
-
求调剂
已经有4人回复
-
【求调剂】085601材料工程专硕 | 总分272 |
已经有4人回复
-
286求调剂
已经有8人回复
-
083000学硕274求调剂
已经有3人回复
-
085701环境工程求调剂
已经有4人回复
-
070300求调剂306分
已经有3人回复
-
266求调剂
已经有12人回复
-
化学调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab解非线性方程组
已经有3人回复
- matlab解非线性方程组程序纠错与正确运行
已经有6人回复
- 非线性方程的线性化问题，大神们进来指教下！！！！感谢了！！！！！！！！！
已经有9人回复
- 1stopt/MATLAB 解非线性方程组问题
已经有4人回复
- 一个关于有约束非线性方程组的求解问题
已经有7人回复
- 求大神帮忙拟合一个非线性方程，求出模型参数
已经有15人回复
- 求助解复杂非线性方程组的好的方法
已经有24人回复
- 用matlab求解一个非线性方程组的解
已经有4人回复
- 非线性方程求解问题，求助各位兄弟姐妹们~~
已经有7人回复
- sos!!1stopt解非线性方程组，一直出现常数定义错误的问题！求大神帮助！
已经有12人回复
- MATLAB求解非线性方程组
已经有5人回复
- 求Matlab解非线性方程代码
已经有12人回复
- matlab无约束非线性方程求最小值的问题
已经有7人回复
- 【求助】非线性方程组的求解问题
已经有6人回复
- 求助matlab 非线性方程组全部解！
已经有6人回复
- 使用fsolve求解非线性方程问题
已经有9人回复
- 需要换迭代方法吗？
已经有7人回复
- 非线性数据拟合的数学原理！！！！！？？？？
已经有30人回复
- 急切求助！一个mathematica求解一个简单一元非线性方程的问题
已经有4人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- matlab求解非线性方程组，求助！
已经有6人回复
- 【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！
已经有5人回复
- 非线性方程组的迭代法（数值计算高手请进）
已经有7人回复
- 【求助】求解非线性方程
已经有8人回复
businiao412
金虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 1431.3
- 红花: 9
- 帖子: 469
- 在线: 58.9小时
- 虫号: 1496676
- 注册: 2011-11-17
- 专业: 石油天然气开采
8楼2014-08-11 16:33:54
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
