| 查看: 3629 | 回复: 11 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
ju5200木虫 (正式写手)
|
[求助]
现在由一组离散的数据 具有周期性，怎么通过这些数据预测下面的情况？ 已有2人参与
|
|
|
比如有4000个数据 每40个数据为一组 每一组数据代表一天的测量值 所有数据截止到十月一日 数据具有周期性 就是说每天测量值都是在增加的 第二天又从一个较低的数增加 以此类推 我想要通过这写数据 预测十月一日以后的情况 应该怎么操作？ 需要先剔除一部分数据吗 ？ 具体预测要用神经网络还是拟合比较好？ 需要用到傅立叶变化吗？ |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
范进中举一文的中心思想
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
-
filecode，4个jtjc了
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 离散数据的求和 怎么求？
已经有6人回复
- matlab可否求离散三维函数（即一组（x，y，z）数据）的偏导数（如：Dz/Dy）？
已经有5人回复
- 【sobereva个人文集】电子定域性的图形分析
已经有56人回复
- matlab离散数据积分
已经有1人回复
4楼2013-12-17 12:26:51
【答案】应助回帖
★
感谢参与，应助指数 +1
xzhdty: 金币+1, 谢谢参与 2013-12-16 15:26:54
感谢参与，应助指数 +1
xzhdty: 金币+1, 谢谢参与 2013-12-16 15:26:54
|
超级菜鸟求助，能不能帮我把这道题转化为lsqcurvefit来解？急！ function KineticsEst1_Diff clear all clc % 动力学数据 t = [0 20 40 60 120 180 300]; CAm = [10 8 6 5 3 2 1]; % 用最小二乘样条拟合法计算微分dCA/dt--使用不经过实验点的B样条插值函数 knots = 3; K = 3; % 三次B样条 sp = spap2(knots,K,t,CAm); pp = fnder(sp); % 计算B样条函数的导函数 dCAdt = fnval(pp,t) % 计算t处的导函数值 rAm = dCAdt; % 绘制浓度拟合曲线 ti = linspace(t(1),t(end),200); CAi = fnval(sp,ti); plot(t,CAm,'ro',ti,CAi,'b-') xlabel('t') ylabel('C_A') legend('实验值','B样条拟合') % 非线性拟合 beta0 = [0.0053 1.39]; [beta,resnorm,residual,exitflag,output,lambda,jacobian] = ... lsqnonlin(@OptObjFunc,beta0,[],[],[],rAm,CAm); ci = nlparci(beta,residual,jacobian); % 参数辨识结果 fprintf('Estimated Parameters:\n') fprintf('\tk = %.4f ± %.4f\n',beta(1),ci(1,2)-beta(1)) fprintf('\tn = %.2f ± %.2f\n',beta(2),ci(2,2)-beta(2)) fprintf(' The sum of the squares is: %.1e\n\n',sum(residual.^2)) % 绘制反应速率拟合曲线 figure plot(t,rAm,'ro',t,Rate(CAm,beta),'b*') xlabel('t') ylabel('dC_Adt') legend('Experiment','Kinetic Model') % ------------------------------------------------------------------ function f = OptObjFunc(beta,rAm,CAm) rAc = Rate(CAm,beta); f = rAc - rAm; % ------------------------------------------------------------------ function rA = Rate(CA,beta) rA = -beta(1)*CA.^beta(2); % -rA = -dCA/dt = k*CA^n, 其中k=beta(1), n=beta(2) |
2楼2013-12-16 13:00:19
ju5200
木虫 (正式写手)
- 应助: 78 (初中生)
- 金币: 1778.6
- 散金: 163
- 红花: 12
- 帖子: 511
- 在线: 65.2小时
- 虫号: 1734187
- 注册: 2012-04-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2013-12-17 09:24:00
ju5200
木虫 (正式写手)
- 应助: 78 (初中生)
- 金币: 1778.6
- 散金: 163
- 红花: 12
- 帖子: 511
- 在线: 65.2小时
- 虫号: 1734187
- 注册: 2012-04-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2013-12-17 19:43:52
40