24小时热门版块排行榜    

Znn3bq.jpeg
查看: 717  |  回复: 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);
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

悠然清风

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by wurongjun at 2016-03-25 17:52:35
你的函数fun1写法不正确: x(i)=vpa();这块有问题!

谢谢,我已经解决了。
4楼2016-03-26 10:36:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

悠然清风

新虫 (初入文坛)

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;
2楼2016-03-25 16:09:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
悠然清风: 金币+5, 还是谢谢,我已经解决了 2016-03-26 18:12:32
你的函数fun1写法不正确: x(i)=vpa([\'x\',num2str(i)]);这块有问题!
善恶到头终有报,人间正道是沧桑.
3楼2016-03-25 17:52:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 308求调剂 +5 VvvvL 2026-04-10 5/250 2026-04-12 10:17 by babysonlkd
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 求调剂 +18 MAX怅惘 2026-04-09 20/1000 2026-04-11 23:31 by zhen~
[考研] 277 数一104,学硕,求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 291分调剂 +5 上岸小莹加油 2026-04-09 6/300 2026-04-11 21:06 by 逆水乘风
[考研] 考研求调剂 +3 ban班小七 2026-04-11 3/150 2026-04-11 20:48 by may_新宇
[考研] 269电子信息求调剂,可转专业 +11 独酌wl 2026-04-06 11/550 2026-04-11 11:12 by 逆水乘风
[考研] 087100初试311求调剂 +4 任雅琴 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[考研] 085402通信工程调剂,有4项学科竞赛国奖(电赛国二),硕士研究生调剂自荐信。 +5 m永o不v言o弃m 2026-04-09 5/250 2026-04-11 09:33 by zhq0425
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 一志愿京区985,085401,与本科专业一致,电子信息工程, +4 阳光开朗的男孩 2026-04-10 4/200 2026-04-10 18:27 by shenrf
[考研] 材料专业344求调剂 +16 hualkop 2026-04-10 21/1050 2026-04-10 17:28 by laoshidan
[考研] 机械还有还有名额吗?太难了 +6 笑笑袁 2026-04-10 6/300 2026-04-10 11:54 by 高维春
[考研] 环境专硕调剂 +16 会说话的肘子 2026-04-06 16/800 2026-04-10 10:30 by asy1wn
[考博] 博士自荐 +7 可可小胖 2026-04-08 7/350 2026-04-10 08:28 by kimhero
[考研] 284求调剂 +7 让我上岸吧阿西 2026-04-09 7/350 2026-04-09 18:59 by haironglove
[考研] 一志愿华东理工085601材料工程303分求调剂 +15 a1708 2026-04-06 15/750 2026-04-08 16:23 by luoyongfeng
[考研] 281求调剂 +10 椰子蘑菇 2026-04-06 10/500 2026-04-08 11:43 by zzucheup
[考研] 336求调剂,一志愿中科大 +9 墨彧 yuyu 2026-04-06 9/450 2026-04-08 11:24 by 想读书的菌菌
[考研] 277、学硕,求调剂 数一104, +11 瓶子PZ 2026-04-07 12/600 2026-04-07 23:30 by 一只好果子?
信息提示
请填处理意见