24小时热门版块排行榜

>论坛更新日志 (3793)
>虫友互识 (886)
>论文投稿 (268)
>休闲灌水 (229)
>导师招生 (169)
>基金申请 (149)
>文献求助 (147)
>硕博家园 (99)
>考博 (78)
>博后之家 (72)
>公派出国 (72)
>教师之家 (61)
>考研 (58)
>论文道贺祈福 (41)
>招聘信息布告栏 (39)
>找工作 (35)

返回列表

悠然清风

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 3
在线: 1.2小时
虫号: 4540434
注册: 2016-03-25

[求助] matlab高维函数如何调用到主函数中？已有1人参与

%最速下降法
%最速下降法
function [x,k,tol,fmin]=minzs(fun,gfun,x0,eps)
format long;
tic
if nargin==3
eps=1.0e-6;
end
tol=1;
syms r;
k=0;
while tol>eps
g0=feval(gfun,x0);
tol=norm(g0);
m=x0+r*g0;
yf=feval(fun,m);
r1=minnd(yf,0);
x1=x0+r1*g0;
x0=x1;
k=k+1;
end
x=vpa(x0,5);
fmin=feval(fun,x0);
fmin=vpa(fmin,5);
format short;

fun1.m函数
function f=fun1(x)
y=ones(1,20);
var=vpa(ones(1,20));
for i=1:20;
x(i)=vpa([\'x\',num2str(i)]);
var(i)=x(i);
y(i)=1.0;
end

m=0;
for i=1:10
m=m+(x(2*i-1))^2+100*(x(2*i))^2;
end

f=0.5*m;

gfun1函数
function gf=gfun1(x)
var=vpa(ones(1,20));
for i=1:20;
x(i)=vpa([\'x\',num2str(i)]);
var(i)=x(i);
if mod(i,2)==0
x(i)=100*x(i);
end
end
gf=x;

clear
x0=ones(1,20);
[x,k,tol,fmin]=minzs(\'fun\',\'gfun\',x0)
为什么不能运行
The following error occurred converting from sym to double:
Error using mupadmex
Error in MuPAD command: DOUBLE cannot convert the input expression
into a double array.

If the input expression contains a symbolic variable, use the VPA
function instead.

Error in gfun1 (line 4)
x(i)=vpa([\'x\',num2str(i)]);

Error in minzs (line 12)
g0=feval(gfun,x0);

回复此楼

» 猜你喜欢

求个博导看看已经有16人回复
想换工作。大多数高校都是评职称时认可5年内在原单位取得的成果吗？已经有8人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
上海工程技术大学【激光智能制造】课题组招收硕士已经有5人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
需要合成515-64-0，50g，能接单的留言已经有4人回复
自荐读博已经有4人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
带资进组求博导收留已经有10人回复

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

opencv数据结构已经有2人回复
MATLAB程序设计与工程应用已经有144人回复
求大神帮忙解题啊！！！已经有4人回复
高等研究院（筹）信息科学与技术研究部招聘启事已经有3人回复
MATLAB2011b可以下载啦.... 已经有21人回复
MathWorks 宣布推出 2011b 版 MATLAB 和 Simulink 已经有5人回复
【分享】Matlab经典教程——从入门到精通【已搜索无重复】已经有797人回复
上海高等研究院（筹）信息科学与技术研究部招聘启事已经有8人回复
急用！！！新虫求助有关阶梯轴的优化设计用matlab编程已经有1人回复
【资源】推荐一些经典的Matlab书籍、教程和工具箱免费下载【已搜索无重复】已经有48人回复
【课件】最优化理论（加拿大维多利亚大学）已经有18人回复

1楼 2016-03-25 15:54:09

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

悠然清风

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 3
在线: 1.2小时
虫号: 4540434
注册: 2016-03-25

function [x,fmin]=minnd(f,x0,eps)
format long;
if nargin==2
eps=1.0e-6;
end
g1=diff(f);
g2=diff(g1);
k=0;
t=1;
while t>eps
h1=subs(g1,findsym(g1),x0);
h2=subs(g2,findsym(g2),x0);
t=abs(h1);
x0=x0-h1/h2;
k=k+1;
end
x=x0;
fmin=subs(f,findsym(f),x);
format short;