| 查看: 705 | 回复: 3 | |||
[求助]
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); |
回复此楼» 猜你喜欢
346求调剂[0856]
已经有7人回复
-
一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！
已经有3人回复
-
310求调剂
已经有3人回复
-
机械专硕299求调剂至材料
已经有4人回复
-
070300化学319求调剂
已经有7人回复
-
08工科 320总分 求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有12人回复
-
307求调剂
已经有3人回复
-
265求调剂
已经有3人回复
-
301求调剂
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- opencv数据结构
已经有2人回复
- MATLAB程序设计与工程应用
已经有144人回复
- 求大神帮忙解题啊！！！
已经有4人回复
- 高等研究院（筹）信息科学与技术研究部招聘启事
已经有3人回复
- MATLAB2011b可以下载啦....
已经有21人回复
- MathWorks 宣布推出 2011b 版 MATLAB 和 Simulink
已经有5人回复
- 【分享】Matlab经典教程——从入门到精通【已搜索无重复】
已经有797人回复
- 上海高等研究院（筹）信息科学与技术研究部招聘启事
已经有8人回复
- 急用！！！新虫求助有关阶梯轴的优化设计用matlab编程
已经有1人回复
- 【资源】推荐一些经典的Matlab书籍、教程和工具箱免费下载【已搜索无重复】
已经有48人回复
- 【课件】最优化理论（加拿大维多利亚大学）
已经有18人回复
2楼2016-03-25 16:09:53
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.4小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
3楼2016-03-25 17:52:35
4楼2016-03-26 10:36:51
