24小时热门版块排行榜

>论坛更新日志 (3925)
>学术会议 (137)
>论文投稿 (106)
>导师招生 (78)
>文献求助 (55)
>休闲灌水 (49)
>博后之家 (46)
>虫友互识 (34)
>硕博家园 (27)
>考博 (27)
>基金申请 (26)
>教师之家 (20)
>考研 (20)
>找工作 (18)
>公派出国 (18)
>论文道贺祈福 (14)

返回列表

落雪啸痕

铜虫 (初入文坛)

应助: 0 (幼儿园)
金币: 95.1
帖子: 5
在线: 1.5小时
虫号: 3098508
注册: 2014-03-29
专业: 机械摩擦学与表面技术

[求助] 为什么我用matlab改变初始参数值，得不到最优解，程序哪里有错误？

syms T t i
N=8;
T=zeros(1,N);
q(1)=3/100*(t/100)^2;
cmin=inf;
for T(1)=1000:-5:200;                %改变T（1）的初值
for i=2:8                            %改变i，获得到其他的q值和T值
q(i)=1.1.^(i-1).*3./100.*((t+symsum(T(i),i,1,i-1))/100)^2;
R=int(q(1),t,0,T(1));
digits(2)
z=vpa(solve(int(q(i),t,0,T(i))-R,T(i)));    %获得T(i)的值
T(i)=z(find(imag(z)==0))                   %筛选，取T(i)的正实数解
end
  C=(31000+500*8*R)/(symsum(T(i),i,1,8)+40);    %目标函数
  if cmin>C
   cmin=C;                                     %将C的较小值赋值给cmin
   Tmin=T(1);                                  %将C较小时的T(1)赋值给Tmin
  end
end
cmin
Tmin

回复此楼

» 猜你喜欢

1楼 2016-12-07 15:22:54

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

相关版块跳转我要订阅楼主落雪啸痕的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3	j5mmowWXNWxv 2026-09-05	3/150	2026-09-06 13:09 by eObJPa4gWmp6
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3	j5mmowWXNWxv 2026-09-05	3/150	2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3	j5mmowWXNWxv 2026-09-05	4/200	2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3	j5mmowWXNWxv 2026-09-05	3/150	2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3	j5mmowWXNWxv 2026-09-05	4/200	2026-09-06 12:29 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3	jPDp0sc2B7zQ 2026-09-05	7/350	2026-09-06 11:40 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3	jPDp0sc2B7zQ 2026-09-05	7/350	2026-09-06 11:20 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题，而不是去凝练问题 +10	瞬息宇宙 2026-09-01	19/950	2026-09-06 11:15 by jiangxuan2004
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3	jPDp0sc2B7zQ 2026-09-05	8/400	2026-09-06 11:09 by eObJPa4gWmp6
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3	jPDp0sc2B7zQ 2026-09-04	7/350	2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3	jPDp0sc2B7zQ 2026-09-04	8/400	2026-09-06 11:00 by eObJPa4gWmp6
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3	jPDp0sc2B7zQ 2026-09-05	4/200	2026-09-06 07:20 by E1iiBLMU2XnD
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6	7ZYIEW2YWn72 2026-09-04	9/450	2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +8	CRJ0Aq2FtVH0 2026-09-03	15/750	2026-09-06 05:48 by E1iiBLMU2XnD
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3	jPDp0sc2B7zQ 2026-09-05	5/250	2026-09-06 02:39 by E1iiBLMU2XnD
[基金申请] 面上没中，邀请各位路过的虫友分析一下分数 +12	jackleilei 2026-09-03	13/650	2026-09-05 18:10 by fireguard
[文学芳草园] 初秋的晨风 +4	阿美_Lml888 2026-09-04	6/300	2026-09-05 13:51 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗？ +6	瞬息宇宙 2026-08-31	9/450	2026-09-05 13:11 by haizai99
[博后之家] 广西大学-广州大学招聘博士后欢迎广大优秀人才！！！ +4	SCSIOyxguo 2026-09-03	7/350	2026-09-05 10:39 by Honglianlian
[基金申请] 要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够 +9	瞬息宇宙 2026-08-31	17/850	2026-09-02 14:04 by anjeeshine