| 5 | 1/1 | 返回列表 |
| 查看: 1801 | 回復: 7 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
matlab程序出錯調試 已有2人參與
|
||
|
我的程序出現(xiàn)問題,但是不知道原因請求各位高手幫忙,不勝感激! function f = statefunction(xx) global mu lam10 Tmp TmE global n x=xx(1:n); x=reshape(x,n,1); u=xx(end-1:end); %x1=x(1),x2=x(2),x3=x(3),x4=x(4),x5=x(5),x6=x(6),x7=x(7),x8=x(8),x9=x(9),x10=x(10),x11=x(11),x12=x(12) % lamb7'=x(13),lamb8'=x(14),lamb9'=x(15),lamb11'=x(16),lamb12'=x(17) %控制變量u1,u2 u1=u(1); u2=u(2); %控制變量u3,u4 uu1=controlE1(x); u3=uu1(1); u4=uu1(2); function y=controlE1(x) %x1=x(1),x2=x(2),x3=x(3),x4=x(4),x5=x(5),x6=x(6),x7=x(7),x8=x(8),x9=x(9),x10=x(10),x11=x(11),x12=x(12) % lamb7'=x(13),lamb8'=x(14),lamb9'=x(15),lamb11'=x(16),lamb12'=x(17) %x1=x(1),x2=x(2),x3=x(3),x4=x(4),x5=x(5),x6=x(6),x7=x(7),x8=x(8),x9=x(9),x10=x(10),x11=x(11),x12=x(12) %lamb1=x(13),lamb2=x(14),lamb3=x(15),lamb4=x(16),lamb5=x(17),lamb6=x(18),... %lamb7=x(19),lamb8=x(20),lamb9=x(21),lamb10=x(22),lamb11=x(23),lamb12=x(24); options = optimset('TolFun',1e-6,'TolCon',1e-8); options = optimset(options,'Algorithm','sqp'); uuE = fmincon(@(xx) myfun1(xx,x),[0;0],[],[],[],[],[-pi;-pi],[pi;pi],@(xx) cons1(xx,x),options); u3=uuE(1); u4=uuE(2); y=zeros(2,1); y(1)=u3; y(2)=u4; function q=myfun1(p,x) global TmE u3=p(1); u4=p(2); a=TmE*cos(u4)/x(8)*(x(15)*cos(u3)-x(8)*x(14)*sin(u3)); b=-TmE*x(14)*sin(u4)*cos(u3)+TmE/x(8)*(x(17)*cos(u4)/cos(x(9))-x(15)*sin(u3)*sin(u4)); q=a^2+b^2; function [c ceq]=cons1(p,x) global TmE u3=p(1); u4=p(2); H11=-TmE*cos(u4)/x(8)*(x(8)*x(14)*cos(u3)+x(15)*sin(u3)); H12=TmE*sin(u4)/x(8)*(x(8)*x(14)*sin(u3)-x(15)*cos(u3)); H22=-TmE*(1/x(8)*(x(15)*sin(u3)*cos(u4)+x(17)*sin(u4)/cos(x(9)))+x(14)*cos(u3)*cos(u4)); c=zeros(2,1); c(1)=H11+H22; c(2)=-(H11*H22-H12*H12); ceq=[]; 如上程序運行時報錯,請問哪位高手知道哪里出錯嗎? ??? Error using ==> sqpLineSearch at 20 Objective function is undefined at initial point. Fmincon cannot continue. Error in ==> fmincon at 816 [X,FVAL,EXITFLAG,OUTPUT,LAMBDA,GRAD,HESSIAN] = sqpLineSearch(funfcn,X,full(A),full(B),full(Aeq),full(Beq), ... Error in ==> controlE1 at 17 uuE = fmincon(@(xx) myfun1(xx,x),[0;0],[],[],[],[],[-pi;-pi],[pi;pi],@(xx) cons1(xx,x),options); Error in ==> statefunction at 14 uu1=controlE1(x); |

至尊木蟲 (著名寫手)
驃騎將軍

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 385分 生物學(071000)求調劑 +3 | qf626 2026-04-01 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京科技,085601總分305求調劑 +9 | 半生瓜! 2026-04-01 | 10/500 |
|
|
[考研] 江蘇科技大學招材料研究生 +4 | Su032713. 2026-04-01 | 5/250 |
|
|
[考研] 320分,材料與化工專業(yè),求調劑 +14 | 一定上岸aaa 2026-03-27 | 18/900 |
|
|
[考研] 332求調劑 +8 | Lyy930824@ 2026-03-29 | 8/400 |
|
|
[基金申請] esi高被引論文是不是能對中標有所加分和幫助呢 +3 | redcom 2026-04-01 | 3/150 |
|
|
[考研] 省雙一流重點一本大學招收調劑 +4 | wwwwffffff 2026-03-31 | 7/350 |
|
|
[考研] 生物學327,求調劑 +4 | 書上的梅子 2026-04-01 | 5/250 |
|
|
[考研] 生物學 296 求調劑 +9 | 朵朵- 2026-03-26 | 11/550 |
|
|
[考研] 267求調劑 +13 | uiybh 2026-03-31 | 13/650 |
|
|
[考研] 一志愿華南師范361分,化學求調劑 +4 | Nicole88888 2026-04-01 | 4/200 |
|
|
[考研] 339求調劑 +5 | zjjkt 2026-03-31 | 5/250 |
|
|
[考研] 材料工程085601數(shù)二英一335求調劑 +5 | 雙馬尾痞老板2 2026-03-31 | 5/250 |
|
|
[考研] 340求調劑 +4 | 希望如此i 2026-03-31 | 4/200 |
|
|
[基金申請] 面上5B能上會嗎? +8 | redcom 2026-03-29 | 8/400 |
|
|
[考研] 289求調劑 +3 | Acesczlo 2026-03-29 | 4/200 |
|
|
[考研] 085404 22408 315分 +5 | zhuangyan123 2026-03-31 | 6/300 |
|
|
[考研] 08工科求調劑286 +5 | tgs_001 2026-03-28 | 5/250 |
|
|
[考研]
|
鐘llll 2026-03-26 | 5/250 |
|
|
[考研] 292求調劑 +13 | 是妍子也是研子 2026-03-30 | 13/650 |
|