Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2009-03-18, 20:58   #1
xqwangnk
初级会员
 
注册日期: 2008-12-09
年龄: 45
帖子: 2
声望力: 0
xqwangnk 正向着好的方向发展
默认 [求助]关于参数推断(最优问题)

小弟初学Matlab, 想求下面问题最优,
function f = eushibor(x)
r=readshibor;
n=numel(r);
y=zeros(n,1);
T = 250;
h = 1/T;
i = 1:n-1;
y(i) = (-1/( 2*log(2*pi*r(i)*h*x(3)^2) ) - (r(i+1)-r(i)-x(1)*(x(2)-r(i))*h)^2/(2*x(3)^2*r(i)*h));
f = -sum(y);
function r=readshibor
r = xlsread('Shibor');

由于对Matlab的一些规则还不是很晓得, 请帮忙指出问题和修改方法,多谢多谢!!
xqwangnk 当前离线   回复时引用此帖
旧 2009-03-18, 20:59   #2
xqwangnk
初级会员
 
注册日期: 2008-12-09
年龄: 45
帖子: 2
声望力: 0
xqwangnk 正向着好的方向发展
默认 回复: [求助]关于参数推断(最优问题)

运行时总显示 Matrix must be square.
xqwangnk 当前离线   回复时引用此帖
旧 2009-03-24, 16:45   #3
云龙九现
初级会员
 
注册日期: 2009-03-20
年龄: 42
帖子: 26
声望力: 17
云龙九现 正向着好的方向发展
默认 回复: [求助]关于参数推断(最优问题)

function f = eushibor(x)
r=readshibor;
n=numel(r);
y=zeros(n,1);
T = 250;
h = 1/T;
i = 1:n-1;
y(i) = (-1/( 2*log(2*pi*r(i)*h*x(3)^2) ) - (r(i+1)-r(i)-x(1)*(x(2)-r(i))*h)^2/(2*x(3)^2*r(i)*h));
f = -sum(y);
function r=readshibor
r = xlsread('Shibor');
上面是你的程序
我帮你修改一下成这样的
function f = eushibor(x)
r = xlsread('Shibor');
r=readshibor;
n=numel(r);
y=zeros(n,1);
T = 250;
h = 1/T;
i = 1:n-1;
y(i) = (-1/( 2*log(2*pi*r(i)*h*x(3)^2) ) - (r(i+1)-r(i)-x(1)*(x(2)-r(i))*h)^2/(2*x(3)^2*r(i)*h));
f = -sum(y);
一般来说,函数调用是要传递和返回参数的。你是想通过readshibor函数来打开一个excel文件,然后将其传回给eushibor函数中的r。但其实没必要这么麻烦。当然,如果你真想这么做应该是这样的
------------------------------------eushibor.m------------------------------------------------------------------------
function f = eushibor(x)
r=readshibor();
n=numel(r);
y=zeros(n,1);
T = 250;
h = 1/T;
i = 1:n-1;
y(i) = (-1/( 2*log(2*pi*r(i)*h*x(3)^2) ) - (r(i+1)-r(i)-x(1)*(x(2)-r(i))*h)^2/(2*x(3)^2*r(i)*h));
f = -sum(y);
-----------------------------------------readshibor.m----------------------------------------------------------
function f1=readshibor()
f1=xlsread(‘Shibor’);
--------------------------------------------main.m-------------------------------------------------------------
x=5;
y=eushibor(x);
---------------------------------------------------------------------------------------------------------------------
:ft:
__________________
不到马里亚纳,不知道大海的深度
不到珠穆朗玛,不知道大山的高度
不到Labfans,不知道自己的无知
云龙九现 当前离线   回复时引用此帖
回复

主题工具
显示模式

发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码



所有时间均为北京时间。现在的时间是 16:57


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.