function [k,en]=TDXJ(f,x,e,h0,m0)
syms x1 x2 m m1;
d=-[diff(f,x1);diff(f,x2)];
flag=1;
k=0;
while(flag==1)
d_temp=subs(d,x1,x(1));
d_temp=subs(d_temp,x2,x(2));
nor=norm(d_temp);
if(nor>= e)
x_temp=x+md_temp;
f_temp=subs(f,x1,x_temp(1));
f_temp=subs(f_temp,x2,x_temp(2));
K=0;
h=h0;
A=m0;
t=2;
T=0;
while T==1
m1=m0+h;
K=K+1;
f_temp1=subs(f_temp,m,m0);
f_temp2=subs(f_temp,m,m1);
if f_temp2<f_temp1
h=t*h;
A=m0;
m0=m1;
else
if K==1
h=-h;
else
T=0;
end
end
end
if m1>A
M=A;
N=m1;
else
M=m1;
N=A;
end
piancha=0.000000000000000000000000000000000000001;
while((N-M)>=piancha)
m3=N+0.382(N-M);
m4=N+0.618*(N-M);
f_temp3=subs(f_temp,m,m3);
f_temp4=subs(f_temp,m,m4);
if f_temp3<f_temp4
M=m4;
else
N=m3;
end
end
m_temp=(M+N)/2;
x=x+m_temp*d_temp;
k=k+1;
else
flag=0;
end
end
en=double(x);
end
无法从 sym 转换为 logical。 出错 TDXJ (line 36) if m1>A 这个问题怎么解决呀
- 写回答
- 好问题 0 提建议
- 追加酬金
- 关注问题
- 邀请回答
-
1条回答 默认 最新
- 技术专家团-Joel 2021-11-18 21:13关注
这是因为你的m1和A都是符号型变量,无法比较大小
可以尝试double(m1)>double(A)
类似的问题都可以这样解决
本回答被题主选为最佳回答 , 对您是否有帮助呢?解决 2无用
悬赏问题
- ¥15 R语言绘制生存曲线时出现警告:All aesthetics have length 1, but the data has 2 rows.
- ¥15 苍穹外卖拦截器token为null
- ¥15 pyqt怎么把滑块和输入框相互绑定,求解决!
- ¥20 wpf datagrid单元闪烁效果失灵
- ¥15 券商软件上市公司信息获取问题
- ¥100 ensp启动设备蓝屏,代码clock_watchdog_timeout
- ¥15 Android studio AVD启动不了
- ¥15 陆空双模式无人机怎么做
- ¥15 想咨询点问题,与算法转换,负荷预测,数字孪生有关
- ¥15 C#中的编译平台的区别影响