24小时热门版块排行榜    

查看: 1922  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

rhinoxixi

金虫 (正式写手)

[交流] 【求助】递推公式的参数拟合 已有2人参与

自己建立的一个模型,是一个递推公式:

Yn=a+b*(Xn)^c+d*Yn-1


要拟合参数,除了自己编程之外,有没有什么软件可以做到的?
matlab行吗?how?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
dongdong3881(金币+2): 多谢指导 2011-02-21 15:17:50
用1stOpt处理,结果更好更稳定,也不需要猜初值:

均方差(RMSE): 0.0154322799500792
残差平方和(SSE): 0.00238155264457617
相关系数(R): 0.971104117667667
相关系数之平方(R^2): 0.943043207351097

参数                  最佳估算
--------------------        -------------
a        0.100663466550055
b        -2.12112269614125E-24
c        -23.8930576866545
d        0.69639916154168
5楼2011-02-21 09:57:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

fspdlh

金虫 (正式写手)

给点数据什么的,说清楚点,好让大家帮你搞啊
2楼2009-04-11 16:41:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rhinoxixi

金虫 (正式写手)

谢谢热心的fspdlh!

具体说明如下 :

n      Xn                        Yn
1     0.82900        0.12500
2     0.81600        0.14600
3     0.64300        0.24300
4     0.57000        0.24700
5     0.47200        0.26800
6     0.36500        0.28400
7     0.24500        0.30700
8     0.11800        0.27100
9     0.37200        0.30000
10   0.55800        0.32400

Y0=0
Y1=a+b*X1^c+d*Y0
Y2=a+b*X2^c+d*Y1
.......

我要拟合出a、b、c、d

其实我知道编程用最小二乘法可以求出来,但是我不是很懂编程,所以问有没有软件算。如果有人fortran很厉害,也可以帮我编个程序来拟合参数,谢谢啦!
3楼2009-04-11 22:31:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fspdlh

金虫 (正式写手)

★ ★ ★ ★ ★
rhinoxixi(金币+5,VIP+0):你实在是太太太好人了,像天使一样~谢谢你! 4-12 00:46
  用MATLAB做的,beta为结果,不知是否符合要求,可能换别的初值会更好,楼主可以多试试。

function [beta r]=datafit(beta)
close all
hold on
load data.txt
X=data(:,2);
Y=data(:,3);
beta=nlinfit(X,Y,'fun',beta);
Y1=fun(beta,X);
r=corrcoef(Y,Y1);
beta=beta';
r=r(1,2);

function Y=fun(beta,X)
a=beta(1);
b=beta(2);
c=beta(3);
d=beta(4);
n=length(X);
Y(1)=a+b*X(1)^c;
for i=2:n
    Y(i)=a+b*X(i)^c+d*Y(i-1);
end
Y=Y';

>> [beta r]=datafit([0 1 0 1])
Warning: Iteration limit exceeded.  Returning results from final iteration.
> In nlinfit at 220
  In datafit at 7

beta =

   -4.4118
    4.5129
    0.0035
    0.7410


r =

    0.9689

>>

[ Last edited by fspdlh on 2009-4-11 at 23:29 ]
4楼2009-04-11 23:26:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 申请回稿延期一个月,编辑同意了。但系统上的时间没变,给编辑又写邮件了,没回复 10+3 wangf9518 2026-03-17 4/200 2026-03-19 23:55 by babero
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 307求调剂 +9 冷笙123 2026-03-17 9/450 2026-03-19 22:44 by 学员8dgXkO
[考研] 288求调剂 +15 于海海海海 2026-03-19 15/750 2026-03-19 22:41 by 学员8dgXkO
[考研] 271材料工程求调剂 +6 .6lL 2026-03-18 6/300 2026-03-19 15:41 by haoshis
[考研] 304求调剂 +3 曼殊2266 2026-03-18 3/150 2026-03-19 14:42 by peike
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 一志愿西南交大,求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 有没有道铁/土木的想调剂南林,给自己招师弟中~ +3 TqlXswl 2026-03-16 7/350 2026-03-17 15:23 by TqlXswl
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
信息提示
请填处理意见