| 查看: 1929 | 回复: 18 | |||
| 【悬赏金币】回答本帖问题，作者Rosalind01将赠送您 10 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
用matlab对微分方程参数拟合已有2人参与
|
|||
|
论文遇到难题，走投无路，希望大神帮忙看一下 需要使用matlab进行参数拟合的，模型为dr=k*(θ-r)dt+σ*r^0.5dW 拟合的参数是k，θ，σ 数据x=[0,0.08,0.17,0.25,0.5,0.75,1] y=[1.805,2.1493,2.1641,2.2051,2.3147,2.3466,2.3873] 用最小二乘法怎么做？或者其他能拟合出来的程序 |
回复此楼» 猜你喜欢
声子谱计算问题
已经有2人回复
-
已结束
已经有23人回复
-
物理学I论文润色/翻译怎么收费?
已经有278人回复
-
Error in routine lambda (100): wrong # or too many modes
已经有0人回复
-
祈福国自然面上
已经有60人回复
-
请问castep计算出的能带图横轴数据缩放因子怎么计算
已经有1人回复
-
Result in Physics 从见刊开始 一般多久可以被SCI收录
已经有21人回复
-
面青地会评时间？？？
已经有8人回复
-
科研求助
已经有0人回复
|
公式没有问题，就是下面这程序怎么用呢，我输入MATLAB就显示未定义变量 % CIR initial parameters estimation x = Model.Data(1:end-1); % Time series of interest rates observations dx = diff(Model.Data); dx = dx./x.^0.5; regressors = [Model.TimeStep./x.^0.5, Model.TimeStep*x.^0.5]; drift = regressors\dx; % OLS regressors coefficients estimates res = regressors*drift - dx; alpha = -drift(2); mu = -drift(1)/drift(2); sigma = sqrt(var(res, 1)/Model.TimeStep); InitialParams = [alpha mu sigma]; % Vector of initial parameters |
15楼2019-02-22 13:43:26
2楼2019-02-22 00:40:17
3楼2019-02-22 00:42:32
hollya
至尊木虫 (文坛精英)
葛城?{美
- 应助: 744 (博后)
- 金币: 43445.1
- 散金: 1143
- 红花: 127
- 沙发: 779
- 帖子: 24904
- 在线: 834.9小时
- 虫号: 458790
- 注册: 2007-11-13
- 性别: MM
- 专业: 计算机软件
4楼2019-02-22 04:13:16
5
