| 查看: 1733 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
新手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))]’; |
回复此楼» 猜你喜欢
科研之友最近访客那出现来自北京的未知用户，什么情况？
已经有20人回复
-
75批博士后面上项目公布了
已经有5人回复
-
名单中没有，评审结果反馈也没有信息，有没有相同的小伙伴指教
已经有6人回复
-
放榜是先短信还是先系统？
已经有9人回复
-
提交了只要等着就好，这样的评审会有吗？
已经有4人回复
-
博后基金意难平……
已经有20人回复
-
2024国自然会评日期
已经有4人回复
-
no news is good news
已经有9人回复
-
两个月二审还没完成
已经有7人回复
-
收到科研之友的邮件
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MATLAB高手帮我看看这个数据存的格式，谢谢
已经有3人回复
- matlab图像处理求光点中心纵坐标
已经有7人回复
- 求基于DCT域的数字图像水印Matlab代码，谢谢各位大神！
已经有3人回复
- matlab解决重心法选址问题，运行不了啊，求助各位大神
已经有10人回复
- 关于Matlab求孔隙率(找最佳阈值)
已经有6人回复
- matlab如何得到图中点对应元素的的位置，请具体点。谢谢。
已经有9人回复
- 纯新手，全方面求助关于MATLAB........
已经有9人回复
- 下载了MATLAB2012可是安装不了。大侠帮帮忙谢谢
已经有13人回复
- 三方协议的问题，请教！万分感谢！
已经有10人回复
- 向各位matlab高手求助了，声音处理方面。
已经有3人回复
- 新手求助Matlab求解系数
已经有10人回复
- 求助各位大侠高手平板污染问题！！！小女子感激不尽
已经有14人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
- 求助一个关于萃取问题，各位大侠帮帮我这新手吧！
已经有10人回复
- 求助各位大侠一个去羟基化的问题，谢谢！
已经有8人回复
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【求助】新手自学GIS,迫切需要帮助,谢谢各位亲啊!有重金······
已经有27人回复
- Matlab画的图经过转换成PDF后不清楚，怎么办？
已经有13人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 求助Origin8.0拟合方程
已经有6人回复
- 【求助】帮我编个程啊~C语言，C++，Matlab，VB...都行~
已经有14人回复
- 【求助】关于概率累计曲线
已经有6人回复
- 【求助】请问Comsol能模拟表面等离子体激元吗
已经有41人回复
virtualzx
木虫 (著名写手)
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
8楼2013-09-26 20:58:19
buhely
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1108.5
- 散金: 151
- 红花: 2
- 帖子: 52
- 在线: 70.3小时
- 虫号: 932783
- 注册: 2009-12-24
- 专业: 控制理论与方法
2楼2013-09-26 10:00:19
3楼2013-09-26 10:23:53
virtualzx
木虫 (著名写手)
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
4楼2013-09-26 11:24:25
50