用matlab做随机模型的参数估计出现错误,求帮助!-经管之家官网!

人大经济论坛-经管之家 收藏本站
您当前的位置> 软件培训>>

Matlab软件培训

>>

用matlab做随机模型的参数估计出现错误,求帮助!

用matlab做随机模型的参数估计出现错误,求帮助!

发布:白塔湖123 | 分类:Matlab软件培训

关于本站

人大经济论坛-经管之家:分享大学、考研、论文、会计、留学、数据、经济学、金融学、管理学、统计学、博弈论、统计年鉴、行业分析包括等相关资源。
经管之家是国内活跃的在线教育咨询平台!

经管之家新媒体交易平台

提供"微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯"等虚拟账号交易,真正实现买卖双方的共赢。【请点击这里访问】

提供微信号、微博、抖音、快手、头条、小红书、百家号、企鹅号、UC号、一点资讯等虚拟账号交易,真正实现买卖双方的共赢。【请点击这里访问】

最近在做毕业论文,想用matlab做极大似然估计,从网上下载到了Princeton大学yacineait-sahalia教授的一些matlab程序。可是弄到matlab里运行会出现错误,对matlab这个软件不大熟,不懂如何改,希望大家帮帮忙。用的模 ...
免费学术公开课,扫码加入


最近在做毕业论文,想用matlab做极大似然估计,从网上下载到了Princeton大学yacine ait-sahalia教授的一些matlab程序。可是弄到matlab里运行会出现错误,对matlab这个软件不大熟,不懂如何改,希望大家帮帮忙。
用的模型是3/2随机模型,dx=(a*x+b*x^2)dt+k*x^(3/2)dW其中我是有x的历史数据的,认为x的变化服从这个随机模型,现在需要用极大似然估计来估计参数a,b,k
下载到的m文件如下:
function output= ModelU1(x,x0,del,param)
% mu(x) = a*x + b*x^2;
%
%
% s(x) = sigma*x^(3/2);
%
%
% g(x) = 2/(sqrt(x)*sigma);
%
%
% lnpX(del, x, x0, 2, exact) = -log(2*Pi*del)/2 - log(s(x)) + cY(g(x), g(x0), -1, exact)/ del + cY(g(x), g(x0), 0, exact) + cY(g(x), g(x0), 1, exact)*del + cY(g(x), g(x0), 2, exact)*(del^2/2);
%
%
% pX(del, x, x0, 2, exact) = exp(lnpX(del, x, x0, 2, exact));
%
%
% cY(y, y0, -1, exact) = (-(1/2))*(y - y0)^2;
%
%
% cY(y, y0, 0, exact) = (a*(-y^2 + y0^2)* sigma^2 + (-8*b + 6* sigma^2)*log(y/y0))/(4* sigma^2);
%
%
% cY(y, y0, 1, exact) = -((1/(24*y*y0* sigma^4))* (48*b^2 + 24*b*(-2 + a*y*y0)* sigma^2 + (9 - 24*a*y*y0 + a^2*y*y0*(y^2 + y*y0 + y0^2))* sigma^4));
%
%
% cY(y, y0, 2, exact) = -((48*b^2 - 48*b* sigma^2 + (9 + a^2*y^2*y0^2)* sigma^4)/ (24*y^2*y0^2* sigma^4));
a = param(1);
b = param(2);
sigma = param(3);
y = 2/(sqrt(x)*sigma);
y0 = 2/(sqrt(x0)*sigma);
sx = sigma*x^(3/2);
output = -log(2*pi*del)/2 - log(sx) ...
+ (-(1/2))*(y - y0)^2 / del ...
+ (a*(-y^2 + y0^2)* sigma^2 + (-8*b + 6* sigma^2)*log(y/y0))/(4* sigma^2)...
+ ( -((1/(24*y*y0* sigma^4))* (48*b^2 + 24*b*(-2 + a*y*y0)* sigma^2 + (9 - 24*a*y*y0 + a^2*y*y0*(y^2 + y*y0 + y0^2))* sigma^4)) )*del ...
+ ( -((48*b^2 - 48*b* sigma^2 + (9 + a^2*y^2*y0^2)* sigma^4)/ (24*y^2*y0^2* sigma^4)) )*(del^2/2);
运行之后显示
??? Input argument "param" is undefined.
Error in ==> ModelU1 at 28
a = param(1);
程序的后面部分应该是做极大似然估计,可是为什么运行后一直在说param没有定义?
我应该怎么改呢?急求帮助!
感激不尽啊!
「经管之家」APP:经管人学习、答疑、交友,就上经管之家!
免流量费下载资料----在经管之家app可以下载论坛上的所有资源,并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库,各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里,你可以找到任何学科方向、有共同话题的朋友。
经管之家(原人大经济论坛),跨越高校的围墙,带你走进经管知识的新世界。
扫描下方二维码下载并注册APP
本文关键词:

本文论坛网址:https://bbs.pinggu.org/thread-2952323-1-1.html

人气文章

1.凡人大经济论坛-经管之家转载的文章,均出自其它媒体或其他官网介绍,目的在于传递更多的信息,并不代表本站赞同其观点和其真实性负责;
2.转载的文章仅代表原创作者观点,与本站无关。其原创性以及文中陈述文字和内容未经本站证实,本站对该文以及其中全部或者部分内容、文字的真实性、完整性、及时性,不作出任何保证或承若;
3.如本站转载稿涉及版权等问题,请作者及时联系本站,我们会及时处理。
经管之家 人大经济论坛 大学 专业 手机版