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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 请教一个非线性方程非线性项转化为多项式的方法
323/4返回列表上一页1234下一页
查看: 3163  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +3 冬十三 2026-03-24 3/150 2026-03-24 21:31 by peike
[考研] 食品专硕 一志愿双一流 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[考研] 0854电子信息求调剂 +6 α____ 2026-03-22 7/350 2026-03-24 19:46 by sunjie8888
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 招08考数学 +7 laoshidan 2026-03-20 16/800 2026-03-24 17:52 by 乌拉儿山脉
[考研] 材料学硕，求调剂 6+3 糖葫芦888ll 2026-03-22 7/350 2026-03-24 17:11 by hello七七
[考研] 材料292调剂 +8 橘颂思美人 2026-03-23 8/400 2026-03-24 16:33 by laoshidan
[考研] 300求调剂，材料科学英一数二 +5 leaflight 2026-03-24 5/250 2026-03-24 16:25 by laoshidan
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +5 niko- 2026-03-22 5/250 2026-03-23 22:01 by fuyu_
[考研] 化学308分求调剂 +3 你好明天你好 2026-03-23 3/150 2026-03-23 20:11 by macy2011
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 一志愿070300浙大化学358分，求调剂！ +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 280分求调剂 一志愿085802 +4 PUMPT 2026-03-22 7/350 2026-03-22 22:13 by 星空星月
[考研] 石河子大学（211、双一流）硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 289材料与化工（085600）B区求调剂 +3 这么名字咋样 2026-03-22 4/200 2026-03-22 17:56 by 云民大李老师
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 一志愿南昌大学，327分，材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
信息提示
关闭
请填处理意见
关闭