24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

园中园

新虫 (正式写手)

应助: 2 (幼儿园)
金币: 33.8
红花: 6
帖子: 339
在线: 221.4小时
虫号: 4199962
注册: 2015-11-05
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 请高手指点：如何将子程序(.m文件)作为参数传递给函数已有2人参与

请高手指点：
   如何将子程序(.m文件)作为一个参数传递给函数

比如： function a=AS(子程序，其他参数)
                        子程序
                           。。。。。。。
            end

回复此楼

» 猜你喜欢

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

我的股市三讲：领会实质把握要领灵活运用稳中求胜已经有4人回复

追求，奋斗，再追求，再奋斗

1楼 2016-06-02 16:18:59

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

aprillf

木虫 (正式写手)

应助: 12 (小学生)
金币: 3720.2
散金: 45
红花: 5
帖子: 562
在线: 248.1小时
虫号: 1688600
注册: 2012-03-13
性别: MM
专业: 地理信息系统

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
月只蓝: 金币+3, 感谢应助 2016-06-30 13:17:18

引用回帖:

3楼: Originally posted by 园中园 at 2016-06-02 19:08:18
大神，能说详细一些吗？...

可以命令行输入 doc fsolve

然后查看fsolve函数说明文档中的例子

用法比如：先定义一个函数m文件待调用

CODE:

function [F,J] = bananaobj(x)

% Evaluate the vector function and the Jacobian matrix for 

% the system of nonlinear equations derived from the general 

% n-dimensional Rosenbrock function.

% Get the problem size

n = length(x);  

if n == 0, error('Input vector, x, is empty.'); end

if mod(n,2) ~= 0, 

   error('Input vector, x ,must have an even number of components.');

end

% Evaluate the vector function

odds  = 1:2:n;

evens = 2:2:n;

F = zeros(n,1);

F(odds,1)  = 1-x(odds);

F(evens,1) = 10.*(x(evens)-x(odds).^2); 

% Evaluate the Jacobian matrix if nargout > 1

if nargout > 1

   c = -ones(n/2,1);    C = sparse(odds,odds,c,n,n);

   d = 10*ones(n/2,1);  D = sparse(evens,evens,d,n,n);

   e = -20.*x(odds);    E = sparse(evens,odds,e,n,n);

   J = C + D + E;

end

然后调用的时候这样使用函数句柄：

CODE:

n = 64;  

x0(1:n,1) = -1.9; 

x0(2:2:n,1) = 2;

options = optimoptions(@fsolve,'Display','iter','Jacobian','on');

[x,F,exitflag,output,JAC] = fsolve(@bananaobj,x0,options);