应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 请教一个非线性方程非线性项转化为多项式的方法
323/4返回列表上一页1234下一页
查看: 3357  |  回复: 31
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

guojiashun

铁杆木虫 (正式写手)
你自己写一个最小二乘的程序来做，应该也不难。
赞一下 回复此楼
21楼2014-08-14 14:31:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guojiashun

铁杆木虫 (正式写手)
又用待定系数法a=inv(X)*y’做了几组不同取点的结果，发现误差不是由matlab的逆矩阵算法误差造成的！这里我犯了一个严重且低级的数学错误！这里不能用待定系数法做！高次多项式的曲线无法用这个方法拟合，不同的取点会导致得到的多项式系数完全不同。看来最小二乘法等拟合方法才是正途！
赞一下 回复此楼
22楼2014-08-14 15:09:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

businiao412

金虫 (正式写手)
引用回帖:
20楼: Originally posted by guojiashun at 2014-08-14 14:28:41
由此可见数值方法里面还是最小二乘法最准确

我昨天晚上也用您说的方法测试了，发现选用不同的X值组成的四个点，得出的四个系数确实有差别的，不过没有仔细整理，我再看看，另外如果取多于五点时，显示五个方程四个未知数求解不出结果，您知道maple软件可以做这样带有位置参量的最小二乘拟合吗？
赞一下 回复此楼
23楼2014-08-14 15:11:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guojiashun

铁杆木虫 (正式写手)
引用回帖:
23楼: Originally posted by businiao412 at 2014-08-14 15:11:56
我昨天晚上也用您说的方法测试了，发现选用不同的X值组成的四个点，得出的四个系数确实有差别的，不过没有仔细整理，我再看看，另外如果取多于五点时，显示五个方程四个未知数求解不出结果，您知道maple软件可以做 ...

可以的，“Maple为我们提供了11个不同的命令实现最小二乘拟合计算”
赞一下 回复此楼
24楼2014-08-14 15:16:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

businiao412

金虫 (正式写手)
引用回帖:
22楼: Originally posted by guojiashun at 2014-08-14 15:09:31
又用待定系数法a=inv(X)*y’做了几组不同取点的结果，发现误差不是由matlab的逆矩阵算法误差造成的！这里我犯了一个严重且低级的数学错误！这里不能用待定系数法做！高次多项式的曲线无法用这个方法拟合，不同的取点 ...

其实我觉得，咱们用最小二乘法拟合的过程中，咱们选的点位置不同，数目不同得到的结果也会有有些差别吧？
赞一下 回复此楼
25楼2014-08-14 16:00:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guojiashun

铁杆木虫 (正式写手)
引用回帖:
25楼: Originally posted by businiao412 at 2014-08-14 16:00:47
其实我觉得，咱们用最小二乘法拟合的过程中，咱们选的点位置不同，数目不同得到的结果也会有有些差别吧？...

应该是，但是数目应该影响最大。
赞一下 回复此楼
26楼2014-08-14 16:03:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

businiao412

金虫 (正式写手)
引用回帖:
26楼: Originally posted by guojiashun at 2014-08-14 16:03:06
应该是，但是数目应该影响最大。...

所以我觉的即使使用最小二乘拟合，也很难得到文章中的结果。
赞一下 回复此楼
27楼2014-08-14 16:39:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guojiashun

铁杆木虫 (正式写手)
引用回帖:
27楼: Originally posted by businiao412 at 2014-08-14 16:39:14
所以我觉的即使使用最小二乘拟合，也很难得到文章中的结果。...

但论文中写的是用最小二乘得到的。
赞一下 回复此楼
28楼2014-08-14 18:21:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guojiashun

铁杆木虫 (正式写手)
我用最小二乘法做了一下：
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))
但是其它几项太复杂了，化简不动.........
赞一下(1人) 回复此楼
29楼2014-08-14 19:03:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guojiashun

铁杆木虫 (正式写手)
把取点密度减小：x=[1-Q:0.1*Q:1]; 可得：
eval(a(4))=-23980501898594795/90564573756653568/conj(Q)^3*conj(Q^(3/2))
             =-0.2648/Q^(3/2);
eval(a(3))=-1/1716*conj(Q^(3/2))*(-2997562737324347/2199023255552-12000*conj(Q)-500*conj(Q)*2^(1/2)+3008*conj(Q)*5^(1/2)+1116*conj(Q)*10^(1/2)+721*conj(Q)*7^(1/2)*10^(1/2)-704*conj(Q)*2^(1/2)*5^(1/2)-381*conj(Q)*3^(1/2)*10^(1/2)+624*conj(Q)*3^(1/2)*5^(1/2))/conj(Q)^3
             =-(1.6841e+003*Q-1.3631e+003)/(1716*Q^(3/2))
             =-(0.9814*Q-0.7943)/Q^(3/2) ≈ -(72*Q-56)/(77*Q^(3/2))

后面实在整不动了，说明最小二乘可定是可行的，但是这个论文作者采用何等手段获得如此整齐精准的表达式呢？
赞一下(1人) 回复此楼
30楼2014-08-14 19:28:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 businiao412 的主题更新
323/4返回列表上一页1234下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 u1xkze8ybe 2026-07-30 7/350 2026-07-31 12:10 by syhhffkrym
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 3/150 2026-07-31 12:08 by syhhffkrym
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 u1xkze8ybe 2026-07-30 6/300 2026-07-31 12:05 by syhhffkrym
[教师之家] 基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊 +6 锦衣卫寒战 2026-07-28 6/300 2026-07-31 11:27 by syl200707
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 小木虫看见有人已经知道结果了 +15 1234567wang 2026-07-25 19/950 2026-07-31 08:30 by ZJTJZ
[基金申请] 2026年国自然面上资助率 +8 布布和一二 2026-07-30 8/400 2026-07-31 01:51 by alongwaytogo
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 微信指数没变化，科研之友没阅读 +14 wangze12014 2026-07-28 17/850 2026-07-30 20:10 by kissu88
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了？ +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[有机交流] 这个自发加氧反应的机理是什么？ 50+4 YUAN2273 2026-07-25 7/350 2026-07-29 12:27 by xiaobenpu
[基金申请] 没消息就是被刷了呗 +13 lambert2014 2026-07-24 14/700 2026-07-29 10:38 by msh145800
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了，我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
关闭
请填处理意见
关闭