版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

a364448242

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 17
在线: 4小时
虫号: 3004526
注册: 2014-02-28
专业: 结构工程

[求助] 请问一下matlab这个是什么问题已有2人参与

是要整成双精度么~~怎么破。菜鸟求教。谢谢

1.jpg

2.jpg

回复此楼

» 猜你喜欢

反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有70人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

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

matlab二重积分的错误提示，指的是什么已经有53人回复
求大神看一下我这个，感觉是软件设置的问题已经有7人回复
以下matlab代码出现问题已经有10人回复
求投过Materials Letters文章或者懂得MATLAB的人帮忙解释一下是什么意思已经有11人回复
以下提示是什么意思？请指教！已经有3人回复
Matlab这个函数是什么意思？已经有8人回复
请教MATLAB中下面的提示是什么意思？已经有4人回复
MATLAB, 代码错误，请大侠帮忙看看。已经有5人回复
如何用matlab或origin画曲线用圆圈隔一段距离标记一下？已经有11人回复
请matlab高手看一下这个问题怎样编程才好，附程序代码，谢谢！已经有9人回复
求助一个matlab关于分类讨论的问题已经有9人回复
帮忙看一下这条语句是什么意思已经有7人回复
请问大牛 'TRACKING' FEATURE 是什么意思，急求答案，谢谢已经有19人回复
问下这个问题怎么解决.....真心啊已经有4人回复
这是什么错误？已经有3人回复
关于Matlab作图的问题已经有14人回复
关于H谱的，这两句话什么意思？重谢！已经有15人回复
linux日志大家帮忙看一下是什么问题？已经有3人回复
【求助】弱弱的问一下Matlab数据提取的问题已经有5人回复
【求助】麻烦大家帮忙看一下德国学校的申请表上这一句话是什么意思呢？已经有18人回复
【求助】这是什么错误？已经有6人回复
牛人帮我看一下，编辑这是什么意思已经有3人回复

1楼 2014-06-18 18:04:35

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

回帖支持 ( 显示支持度最高的前 50 名 )

生如夏花ZJU

铜虫 (初入文坛)

应助: 4 (幼儿园)
金币: 6.1
帖子: 11
在线: 15.4小时
虫号: 2624968
注册: 2013-08-30
性别: GG
专业: 流体力学

【答案】应助回帖

★
感谢参与，应助指数 +1
fegg7502: 金币+1, 应助指数+1, 3ks 2014-06-19 08:04:35

这个调试一下应该就能看出来的。在出错的那一行前面设个断点，运行到断点处就可以在workspace里看x的属性。
因为你这个前面调用过fun(x0)没问题，那么问题最可能是在x= fminunc('fun',x0);这一句了。

赞一下(1人)

回复此楼

生如夏花之绚烂，死如秋叶之静美

7楼2014-06-18 23:03:05

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

普通回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68123.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★
感谢参与，应助指数 +1
fegg7502: 金币+1, 3ks 2014-06-19 08:04:07

fun 没定义
给出完整的程序看看。

赞一下

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

2楼2014-06-18 18:49:08

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

a364448242

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 17
在线: 4小时
虫号: 3004526
注册: 2014-02-28
专业: 结构工程

★
fegg7502: 金币+1, 鼓励交流 2014-06-19 08:04:15

引用回帖:

2楼: Originally posted by 月只蓝 at 2014-06-18 18:49:08
fun 没定义
给出完整的程序看看。

clc
tic
x0=input('请输入初始点(形式如[0 0])：');
x0=x0';
disp('变尺度法求解无约束最优化问题，本程序以Rosenbrock函数为例');
disp('1-----Armijo搜索');
disp('2-----Armijo_goldstein搜索');
tt=input('请您选择搜索方法\n');
disp('1-----DFP校正公式');
disp('2-----Symmetric Rank 1校正公式');
disp('3-----BFGS校正公式');
zz=input('请您选择所用的系数公式\n');
disp('正在计算，请稍等……');
f0=fun(x0);
g0=gfun(x0);
p0=-gfun(x0);
bcd(x0,tt,zz);
disp('---------调用系统函数相比较-----------');
x= fminunc('fun',x0);
f=fun(x);
disp('最小值时的迭代点 :');
      x
      disp('函数的最小值为:');
      f

disp('CPU运行时间：')
disp('-------------结束-------------');
toc

赞一下

回复此楼

3楼2014-06-18 19:46:28

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68123.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

引用回帖:

3楼: Originally posted by a364448242 at 2014-06-18 19:46:28
clc
tic
x0=input('请输入初始点(形式如)：');
x0=x0';
disp('变尺度法求解无约束最优化问题，本程序以Rosenbrock函数为例');
disp('1-----Armijo搜索');
disp('2-----Armijo_goldstein搜索');
tt=input('请 ...

gfun的定义呢

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

4楼2014-06-18 20:31:23

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68123.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

引用回帖:

以及bcd的定义

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

5楼2014-06-18 20:31:55

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68123.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★
fegg7502: 金币+1, 应助指数+1, 鼓励交流 2014-06-19 08:04:27

你这个函数涉及的函数比较多，你把所有的函数都放在matlab的工作目录下试试看看。

赞一下

回复此楼

» 本帖已获得的红花（最新10朵）

a364448242

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

6楼2014-06-18 20:33:35

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

a364448242

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 17
在线: 4小时
虫号: 3004526
注册: 2014-02-28
专业: 结构工程

引用回帖:

4楼: Originally posted by 月只蓝 at 2014-06-18 20:31:23
gfun的定义呢...

%(求梯度)
function g=gfun(x0)
syms x1 x2
f=100*(x2-x1^2)^2+(1-x1)^2;
gg=jacobian(f,[x1;x2]);
g=(subs(gg,{x1,x2},{x0(1),x0(2)}))';

赞一下

回复此楼

8楼2014-06-19 09:24:05

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

a364448242

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 17
在线: 4小时
虫号: 3004526
注册: 2014-02-28
专业: 结构工程

引用回帖:

5楼: Originally posted by 月只蓝 at 2014-06-18 20:31:55
以及bcd的定义...

function [x,f,times]=bcd(x0,tt,zz)
xx1=-1:0.1:2;
xx2=-1:0.1:2;
%----------绘图函数
set(0,'defaultfigurecolor',[1 .8 0]);
[x1,x2]=meshgrid(xx1,xx2);
f=100*(x2-x1.^2).^2+(1-x1).^2;
contour(x1,x2,f,20);
hold on
title('变尺度法','FontSize',16,'fontweight','bold');
strl(1)={'信科012班'};
strl(2)={'9号;23号'};
uicontrol('Style','text','FontSize',12,'Position',[75 320 120 38],...
'String',strl);
%---------------变尺度法搜索最小值函数
%---------------tt表选择搜索方法系数
%---------------zz表选择系数
x=x0;
h0=[1 0;0 1];
epis=0.000001;
f=fun(x0);
g=gfun(x0);
times=0;
plot(x0(1),x0(2),'k*')
while norm(g)>=epis
p0=-h0'*g;
times=times+1;
x0=x;
f0=f;
g0=g;
switch tt       %--------------------选择搜索方法
      case 1       %--------------------Armijo法
x=s1(x0,p0);
      case 2       %--------------------Armijo_goldstein法
x=s2(x0,p0);
  otherwise
return
end

plot(x(1),x(2),'g*')
plot([x0(1), x(1)],[x0(2), x(2)],'r')
f=fun(x);
g=gfun(x);
s0=x-x0;
y0=g-g0;
if norm(g)<epis
      disp('最后的迭代步数为:')
      times
      disp('最小值时的迭代点 :');
      x
      disp('函数的最小值为:');
      f
      break;
   end
switch zz
      case 1
      h=h0+(s0*s0')/(s0'*y0)-(h0*y0*y0'*h0)/(y0'*h0*y0);             %-------DFP校正公式
      case 2
      h=h0+((s0-h0*y0)*(s0-h0*y0)')/(y0'*(s0-h0*y0));                   %-------Symmetric Rank 1校正公式
      case 3
      h=h0+((1+(y0'*h0*y0)/(s0'*y0))*s0*s0'-h0*y0*s0'-s0*y0'*h0)/(s0'*y0); %----------BFGS校正公式
      otherwise
         return
   end
      h0=h;

   end

赞一下

回复此楼

9楼2014-06-19 09:26:33

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