matlab优化求解

nkkopd

关于matlab优化求解的问题请教大家；程序报错截图如下，我分析报错的原因是迭代点（s,d1）改变时，（expect1，expect2）未改变，具体问题描述如下：
1.expect1，expect2是s,d1的函数，当fsolve函数按照优化算法找到新的点（s',d1'）时，我的问题里时需要expect1，expect2也变化，这样fsolve的目标函数就也变化了；
4 n8 T# Y& l# ^2.根据上述1的描述，这个优化问题是否应该将（expect1，expect2）也纳入优化求解目标中去？这样该如何实现？
关于这样的编程问题思考了很久也没有编程思路，请求大家给予指导，我在代码中相应部分也对问题进行了表述，如果我的问题没有表述清楚，十分欢迎追问，先谢过大家。
1 h$ p! W1 q5 l7 c9 m' Hclc;
clear all;
3 k. I# ?6 O! X( x5 K9 ^%%%  赋值外生给定的参数  %%%
3 v; T3 I! W) m+ |4 @7 B% ]wn=0.3;
wt=1-wn;
. J! J! k' u5 R6 Mr=0.05;
- W9 _5 r8 |+ w( J' k2 y) sbeta=1/(1+r);
; i: @8 P# Z. D; o$ \- A+ zepislon=0.05;
phi=0.3;
! _" P& F- Z# F0 T6 I4 r6 x# ]6 Uy1=1;
, d% \) }: ]3 p  I9 tyn2=1;
$ t7 f! l4 q! ]y2mean=1;
y3=1;
+ C. [, v/ C, M, @* z" e. dtheta=0.9;

3 z* y  C4 p, U7 U( h% 赋值(y2,s,d1)
2 i: p' v- e0 L/ _; S3 [y2=0.95:0.01:1.05;
+ G; m6 g0 s5 r! ~s=0.1; %给s,d1的任意初值
d1=0.1;
x0 =[s,d1];

! c+ ?0 A% b, L: n5 E/ f%在这个for循环里得到均值expect1,expect2，改变(s,d1),均值expect1,expect2也改变

for i=1:1:length(y2);
- J6 {6 A% k, a0 o; M% qs=0.01;
& `& n8 s, g5 Q' Qd1=0.1;

%%% 判断是否紧约束 %%%
standard(i)=y2(i)*((beta/(1+beta)+phi/wt/(1-phi*wn/wt))-y3/((1+beta)*(1+r)))/(beta/(1+beta)+phi*wn/wt/(1-phi*wn/wt));
  C# B$ z! D5 B" E8 D4 ucheck(i)=d1/(1-s);
if check(i)>standard(i)
3 ~* Q4 ~" H9 |6 ?# u8 j$ w. Q    ct2(i)=((1-s)*y2(i)-d1+phi*(1-s)*y2(i))/(1-phi*wn/wt);
6 ~- [4 k/ x8 r7 K  s# S  p  N: Y/ q    d2(i)=phi*(1+r)*(wn/wt*ct2(i)+(1-s)*y2(i));
    ct3(i)=(1-s)*y3-d2(i);
    mu(i)=wt/((1-s)*y2(i)-d1+d2(i)/(1+r))-beta*wt*(1+r)/((1-s)*y3-d2(i));
% Q) L  p# C2 p2 G
1 m2 H6 z# o$ [, F9 }1 s: ~  r1 g- \else
    ct2(i)=((1-s)*y2(i)-d1+(1-s)*y3/(1+r))/(1+beta);
    d2(i)=((1-s)*y3-beta*(1+r)*(1-s)*y2(i)+beta*(1+r)*d1)/(1+beta);
1 v; a0 A9 h3 L! t: Z    ct3(i)=(1-s)*y3-d2(i);

    mu(i)=0;
end
. }  c2 R, o$ c. \! g. u( te1(i)=wt/ct2(i);%每给定一个y2(i)，计算得到一个e1(i)，e2(i)
* O) i+ Y& L% V( w2 n" Pe2(i)=wt/ct2(i)*y2(i)+mu(i)*phi*y2(i)+beta*wt/ct3(i)*y3;

# p7 T+ O8 T  H( K9 Q4 S8 x7 M2 x%%%求均值%%%
expect1=mean(e1(1,);
expect2=mean(e2(1,);
end

+ [: e) {" U( i6 Z6 S%我的问题中，expect1，expect2实质上是s,d1的函数
$ ~& `# A! ]: D% P5 N! g%在用fsolve优化求解的时候，每迭代一个新的（s,d1),都应该运行一次上述for这个循环，得到新的expect1、expect2再进行优化求解，
7 F- F  P, \4 Y" N0 U# {/ h4 N%也即（s,d1)变化时，expect1、expect2也要变化，如此fsolve的目标函数也变化
# g+ B& B0 O/ Y%这样一个问题没有编程思路，请求大家给予指导，如果我的问题没有表述清楚，十分欢迎追问，谢谢。
3 F6 z7 L9 S) _" h
[x,fval,exitflag,output,]=fsolve(@(x)foc(x,expect1,expect2),x0);

' s" L$ f9 U9 R- ufunction f=foc(x,expect1,expect2)
wn=0.3;
wt=1-wn;
: t; S6 A9 I3 s- e1 o6 l" @) lr=0.05;
8 ]) y/ C5 B1 S  f2 dbeta=1/(1+r);
epislon=0.05;
phi=0.3;
: g3 c5 n# M6 G2 R3 P5 {y1=1;
4 e4 i0 s, g0 h# n5 q4 Zyn2=1;
7 `1 {/ A/ }- Oy2mean=1;
# B, |$ a+ E! Hy3=1;
theta=0.9;
s=0.01;
d1=0.1;
) c# `4 a3 k" h) Gct1=s*(1-theta)*y1+d1/(1+r);% ct1也与s,d1相关；

f(1)=wt/ct1-beta*(1+r)*expect1;
f(2)=wt/ct1*y1*(1-theta)-beta*expect2;
end


2 j" \. e0 F# r3 L9 O

pTDbn25

你这个写的啥，你foc的函数的变量f1和f2里面x的变量呢，你这个方程写的有很大问题

kekek

来学习一下

zaiyiaaaa