24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2350  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

jjliu2012

新虫 (小有名气)

[求助] 新手matlab求助,谢谢各位

最优控制中的最速下降法,若是求不出函数几乎函数梯度的解析表达式,最速下降法的程序该怎么样处理呢?
以下是我们平时用的最速下降法程序

function [x,val,k]=grad(fun,gfun,x0)
%功能: 用最速下降法求解无约束问题: min f(x)
%输入: x0是初始点, fun, gfun分别是目标函数和梯度
%输出: x, val分别是近似最优点和最优值, k是迭代次数.
maxk=5000; %最大迭代次数
rho=0.5;sigma=0.4;
k=0; epsilon=1e-5;
while(k¡maxk)
g=feval(gfun,x0); %计算梯度
d=-g; %计算搜索方向
if(norm(d)¡epsilon), break; end
m=0; mk=0;
while(m¡20) %Armijo搜索
if(feval(fun,x0+rho^m*d)¡feval(fun,x0)+sigma*rho^m*g’*d)
mk=m; break;
end
m=m+1;
end
x0=x0+rho^mk*d;
k=k+1;
end
x=x0;
val=feval(fun,x0);

%这里的fun是函数,f是函数,g是f函数的梯度,这里的f和g都有解析的表达式
function f=fun(x)
f=100*(x(1)^2-x(2))^2+(x(1)-1)^2;
function g=gfun(x)
g=[400*x(1)*(x(1)^2-x(2))+2*(x(1)-1), -200*(x(1)^2-x(2))]’;
回复此楼

» 猜你喜欢

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

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

virtualzx

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
你这里fun是定积分,不是t的函数,换个例子来吧

PS:如果换成不定积分fun(t)=int[x(s)^2.ds,a,t], 那么gfun(t)=x(t)^2=int[f,a0,t]^2
7楼2013-09-26 20:50:51
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

buhely

银虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
如果你的fun是用明确的数学语言描述的,但不能直接算出gfun,可以考虑用matlab的符号变量编程,然后算出gfun的表达式;
再不行貌似就只能再添加各种启发式搜素算法了
2楼2013-09-26 10:00:19
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jjliu2012

新虫 (小有名气)

引用回帖:
2楼: Originally posted by buhely at 2013-09-26 10:00:19
如果你的fun是用明确的数学语言描述的,但不能直接算出gfun,可以考虑用matlab的符号变量编程,然后算出gfun的表达式;
再不行貌似就只能再添加各种启发式搜素算法了

这里的x满足的是一个非线性的微分方程,x就没有解析解,fun是关于x的,也没有解析表达式啊
3楼2013-09-26 10:23:53
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)

引用回帖:
3楼: Originally posted by jjliu2012 at 2013-09-25 21:23:53
这里的x满足的是一个非线性的微分方程,x就没有解析解,fun是关于x的,也没有解析表达式啊...

你能不能把你得问题说得再明白点,贴的程序貌似和你的问题也不相关
4楼2013-09-26 11:24:25
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 283分材料与化工求调剂 +12 罗KAKA 2026-04-02 12/600 2026-04-02 16:50 by guanxin1001
[考研] 材料调剂 +8 一样YWY 2026-04-02 8/400 2026-04-02 12:30 by 173785543
[考研] 085601一志愿中山大学深圳材料工程330求调剂 +8 pipiver 2026-03-30 8/400 2026-04-02 12:01 by ms629
[考研] 085900土木水利336分求调剂 +4 Zhangjiangj 2026-03-31 6/300 2026-04-02 11:40 by 1753564080
[考研] 348求调剂 +6 小懒虫不懒了 2026-03-27 7/350 2026-04-02 11:18 by sanrepian
[考研] 324分 085600材料与化工 +20 呆鹅oor 2026-03-27 20/1000 2026-04-02 10:13 by oooqiao
[考研] 一志愿346上海大学生物学 +3 上海大学346调剂 2026-04-01 3/150 2026-04-02 08:36 by w虫虫123
[考研] 279求调剂 +7 莫xiao 2026-04-01 7/350 2026-04-01 22:05 by 客尔美德
[考研] 353求调剂 +4 拉钩不许变 2026-04-01 4/200 2026-04-01 18:10 by 记事本2026
[考研] 求调剂 +4 图鉴212 2026-03-30 5/250 2026-04-01 15:32 by 图鉴212
[考研] 安徽大学专硕生物与医药专业(086000)324分,英语已过四六级,六级521,求调剂 +10 美味可乐鸡翅 2026-03-26 12/600 2026-04-01 11:17 by syh9288
[考研] 289求调剂 +7 BrightLL 2026-03-29 7/350 2026-03-31 22:05 by 544594351
[考研] 277跪求调剂 +8 1915668 2026-03-27 13/650 2026-03-31 14:58 by 王亮_大连医科大
[考研] 英一数一总分334求调剂 +4 陈阳坤 2026-03-31 4/200 2026-03-31 14:22 by 记事本2026
[考研] 282求调剂 +3 呼吸都是减肥 2026-03-31 3/150 2026-03-31 10:00 by wangjy2002
[考研] 吉大生物学326分求调剂 +3 sunnyupup 2026-03-31 3/150 2026-03-31 09:28 by longlotian
[考研] 085701环境工程求调剂 +11 多久上课 2026-03-27 12/600 2026-03-30 21:21 by 研究僧导导
[考研] 调剂 +4 GK72 2026-03-30 4/200 2026-03-30 20:32 by dick_runner
[有机交流] 考研调剂 +8 watb 2026-03-26 8/400 2026-03-30 18:40 by 544594351
[考研] 352分 化工与材料 +5 海纳百川Ly 2026-03-27 5/250 2026-03-28 03:39 by fmesaito
信息提示
请填处理意见