matlab优化求解

nkkopd

关于matlab优化求解的问题请教大家；程序报错截图如下，我分析报错的原因是迭代点（s,d1）改变时，（expect1，expect2）未改变，具体问题描述如下：
1.expect1，expect2是s,d1的函数，当fsolve函数按照优化算法找到新的点（s',d1'）时，我的问题里时需要expect1，expect2也变化，这样fsolve的目标函数就也变化了；
3 m4 B! Z1 S0 X- I2.根据上述1的描述，这个优化问题是否应该将（expect1，expect2）也纳入优化求解目标中去？这样该如何实现？
. |  T. l7 f1 f3 i关于这样的编程问题思考了很久也没有编程思路，请求大家给予指导，我在代码中相应部分也对问题进行了表述，如果我的问题没有表述清楚，十分欢迎追问，先谢过大家。
: M( w3 x7 l1 b# v" Z! aclc;
clear all;
  r* R; _1 _! L- X3 w( H8 r( J. w%%%  赋值外生给定的参数  %%%
wn=0.3;
wt=1-wn;
r=0.05;
beta=1/(1+r);
epislon=0.05;
phi=0.3;
y1=1;
( c* J( B$ P6 i+ syn2=1;
y2mean=1;
y3=1;
theta=0.9;

6 ^/ z& M& i/ J$ }% 赋值(y2,s,d1)
y2=0.95:0.01:1.05;
& l- ^1 O2 N- @7 R  ss=0.1; %给s,d1的任意初值
: `! \- D% ^; y% K* {6 Qd1=0.1;
( b. ?1 @* h! s6 ~: u' u! Dx0 =[s,d1];

; c/ J. C3 k8 p- q& ?" E% Z%在这个for循环里得到均值expect1,expect2，改变(s,d1),均值expect1,expect2也改变
2 ^- }% z# l5 S) |) }$ q, Z6 j: U7 a
for i=1:1:length(y2);
s=0.01;
2 E( d* _% a( |d1=0.1;
; J$ J" L& w; [5 B1 a
%%% 判断是否紧约束 %%%
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));
check(i)=d1/(1-s);
  q2 y- d6 V. |9 l) Fif check(i)>standard(i)
    ct2(i)=((1-s)*y2(i)-d1+phi*(1-s)*y2(i))/(1-phi*wn/wt);
& t5 g( E2 t: R0 ]7 Z    d2(i)=phi*(1+r)*(wn/wt*ct2(i)+(1-s)*y2(i));
    ct3(i)=(1-s)*y3-d2(i);
( g4 }9 f, y. R( U    mu(i)=wt/((1-s)*y2(i)-d1+d2(i)/(1+r))-beta*wt*(1+r)/((1-s)*y3-d2(i));

/ M8 F1 Q9 q1 d; ?2 c  B$ uelse
$ X. G+ ~  ?3 a- L  f) r8 x    ct2(i)=((1-s)*y2(i)-d1+(1-s)*y3/(1+r))/(1+beta);
. L3 w/ \  P; }$ T- {6 v    d2(i)=((1-s)*y3-beta*(1+r)*(1-s)*y2(i)+beta*(1+r)*d1)/(1+beta);
! g- {& d+ `8 ~% }6 N% [    ct3(i)=(1-s)*y3-d2(i);
4 e; n- u, v4 i' K
; S6 v2 K6 P0 G    mu(i)=0;
end
& n$ e2 j  n1 R) {e1(i)=wt/ct2(i);%每给定一个y2(i)，计算得到一个e1(i)，e2(i)
+ s8 g0 X1 h6 ?! p5 h, he2(i)=wt/ct2(i)*y2(i)+mu(i)*phi*y2(i)+beta*wt/ct3(i)*y3;
' O4 p' c  m. r8 e
%%%求均值%%%
expect1=mean(e1(1,);
5 z3 I( x% s. r7 w& a- y0 Nexpect2=mean(e2(1,);
, z3 m( G" c; A! J4 Fend
2 G; G- [1 _: v$ h! _
%我的问题中，expect1，expect2实质上是s,d1的函数
, X/ G2 D: ~* S! x& a+ n7 U3 c%在用fsolve优化求解的时候，每迭代一个新的（s,d1),都应该运行一次上述for这个循环，得到新的expect1、expect2再进行优化求解，
; E9 {; q5 ~" X' E%也即（s,d1)变化时，expect1、expect2也要变化，如此fsolve的目标函数也变化
%这样一个问题没有编程思路，请求大家给予指导，如果我的问题没有表述清楚，十分欢迎追问，谢谢。

[x,fval,exitflag,output,]=fsolve(@(x)foc(x,expect1,expect2),x0);
; o3 w0 u- c1 v- i
0 P! d( \) S: p3 h- e& W9 K. Zfunction f=foc(x,expect1,expect2)
5 k4 \' e: D% b. {/ r7 Vwn=0.3;
3 x$ C9 e8 K# }. c' K+ ?wt=1-wn;
r=0.05;
3 O! u& u7 c& mbeta=1/(1+r);
epislon=0.05;
; K. A6 s$ n6 `phi=0.3;
y1=1;
yn2=1;
* z* A5 y! [. z* M) [7 t# ?: gy2mean=1;
9 {" C, ]: Z; Q7 k4 X1 X) o7 ny3=1;
7 t4 r$ [/ c/ k& V' _" t: Btheta=0.9;
s=0.01;
d1=0.1;
& q+ v! Q4 R3 ~. |/ a  Ect1=s*(1-theta)*y1+d1/(1+r);% ct1也与s,d1相关；
1 C/ {5 z% {3 ~
2 i; P, X" c5 o( K) `1 zf(1)=wt/ct1-beta*(1+r)*expect1;
# `- i5 f; J, If(2)=wt/ct1*y1*(1-theta)-beta*expect2;
end
0 T2 `! f0 H  P' f+ b  d
" V. A- C3 N: J# W& q8 G
, @  h! y' J" g

pTDbn25

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

kekek

来学习一下

zaiyiaaaa