| 查看: 953 | 回复: 1 | ||
renwanbo123新虫 (初入文坛)
|
[求助]
Matlab求解钛卷径向导热系数-黄金分割法
|
|
各位大侠，小弟正在做钢卷退火的数值仿真，需要通过实验和编程来求解钢卷的径向导热系数，我参考的文献是《铝卷径向等效导热系数的仿真_李建平2009》、《钛带卷真空退火传热过程的数值模拟_陈临夏2012》，文献中采用的黄金分割法来优化目标函数，我搞不明白，跪求大侠指点！ 小弟现在只有5个金币，全部金币奉上，以示诚意！ 以下我找的住程序和黄金分割法Matlab代码，但是我不知道文献中的目标函数J(k)怎么输入，然后怎么进行迭代的。 主程序： syms x a b e; %定义变量 a=input('搜索区间的第一点\a='); %确定搜索区间 b=input('搜索区间的第二点\b='); e=input('搜索精度\ne='); %收敛精度 disp('需求的优化函数f=f(x)，调用xmin=golden(f,a,b,e)'); 黄金分割法代码： function xmin=golden(f,a,b,e) k=0; a1=b-0.618*(b-a); %插入点的值 a2=a+0.618*(b-a); while b-a>e %循环条件 y1=subs(f,a1); y2=subs(f,a2); if y1>y2 %比较插入点的函数值的大小 a=a1; %进行换名 a1=a2; y1=y2; a2=a+0.618*(b-a); else b=a2; a2=a1; y2=y1; a1=b-0.618*(b-a); end k=k+1; end %迭代到满足条件为止就停止迭代 xmin=(a+b)/2; fmin=subs(f,xmin) %输出函数的最优值 fprintf('k=\n'); %输出迭代次数 disp(k);  1.png  2.png |
回复此楼» 猜你喜欢
今天放榜吗？
已经有13人回复
-
filecode，4个jtjc了
已经有9人回复
-
今天基金会出结果吗？20260819
已经有8人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有115人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
-
时间戳变了，能看出什么问题？
已经有14人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有30人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有7人回复
-
什么时候开奖？
已经有6人回复
renwanbo123
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1095.4
- 帖子: 29
- 在线: 10.5小时
- 虫号: 1644270
- 注册: 2012-02-26
- 性别: GG
- 专业: 金属材料
2楼2015-07-31 09:30:36
