matlab优化求解的问题

Uifhjvv

关于matlab优化求解的问题请教大家；程序报错截图如下，我分析报错的原因是迭代点（s,d1）改变时，（expect1，expect2）未改变，具体问题描述如下：
1 R$ j. _  a# u4 u1.expect1，expect2是s,d1的函数，当fsolve函数按照优化算法找到新的点（s',d1'）时，我的问题里时需要expect1，expect2也变化，这样fsolve的目标函数就也变化了；
$ z* @8 |+ N" n5 x2.根据上述1的描述，这个优化问题是否应该将（expect1，expect2）也纳入优化求解目标中去？这样该如何实现？
clc;
/ v* Y3 M2 o$ U& Q/ lclear all;
%%%  赋值外生给定的参数  %%%
wn=0.3;
wt=1-wn;
3 b. n) K( N; Q" J- I$ ^* H) u! Hr=0.05;
+ ~( V; g% K; s) j+ G' L9 lbeta=1/(1+r);
( D. }4 W& D0 s5 x4 sepislon=0.05;
( Y0 \- O* u" C! E% P  q1 ophi=0.3;
y1=1;
: C6 S! r: v# B; \& Zyn2=1;
$ G: M6 K% ~# j" s  [- i! b8 qy2mean=1;
  R' h8 I" B6 k" r7 [y3=1;
1 j* C% K6 a( k3 [3 x) Jtheta=0.9;


) t, S5 u4 J6 ^! ^, c% 赋值(y2,s,d1)
y2=0.95:0.01:1.05;
" d9 }; r6 [1 ?9 g  j# bs=0.1; %给s,d1的任意初值
d1=0.1;
x0 =[s,d1];

4 O! B/ R# v6 D
%在这个for循环里得到均值expect1,expect2，改变(s,d1),均值expect1,expect2也改变


for i=1:1:length(y2);
# [: I$ _  o  j. Qs=0.01;
! q: u) {% v5 hd1=0.1;
7 Y+ k8 L1 G+ r$ D2 d8 d+ {

%%% 判断是否紧约束 %%%
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));
7 P8 B: [( j7 T7 ]" fcheck(i)=d1/(1-s);
if check(i)>standard(i)
5 S' f) t4 S  k" t; v9 q    ct2(i)=((1-s)*y2(i)-d1+phi*(1-s)*y2(i))/(1-phi*wn/wt);
! `' `( d! v9 T( Y, K    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));
   
" y6 x" L6 u6 x- _" P. gelse
    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);
7 k  S7 k) l0 u8 T    ct3(i)=(1-s)*y3-d2(i);
+ E" g( ?/ C" q3 o

    mu(i)=0;
& }1 y* U, H* b/ \" f  ^# R* yend
e1(i)=wt/ct2(i);%每给定一个y2(i)，计算得到一个e1(i)，e2(i)
e2(i)=wt/ct2(i)*y2(i)+mu(i)*phi*y2(i)+beta*wt/ct3(i)*y3;


%%%求均值%%%
expect1=mean(e1(1,：));
expect2=mean(e2(1,：));
end


%我的问题中，expect1，expect2实质上是s,d1的函数
%在用fsolve优化求解的时候，每迭代一个新的（s,d1),都应该运行一次上述for这个循环，得到新的expect1、expect2再进行优化求解，
8 R8 ?( z" k! `. Q; V%也即（s,d1)变化时，expect1、expect2也要变化，如此fsolve的目标函数也变化
%这样一个问题没有编程思路，请求大家给予指导，如果我的问题没有表述清楚，十分欢迎追问，谢谢。

/ ]* ]5 b4 ?$ }4 E0 Y* t7 p" p
" V/ h9 e3 Z  P' C7 f  V7 r* ?[x,fval,exitflag,output,]=fsolve(@(x)foc(x,expect1,expect2),x0);


function f=foc(x,expect1,expect2)
wn=0.3;
; [) }- `# F5 |8 M' P6 J2 Gwt=1-wn;
5 A2 _5 [* v$ M) ~r=0.05;
* q: ~  j. w. p: S0 \- @2 Cbeta=1/(1+r);
epislon=0.05;
; w2 h$ `6 d* R) h3 s: ophi=0.3;
; L( N$ q7 q0 Uy1=1;
4 G) ?7 r3 E+ F+ D, X4 _yn2=1;
$ V  O/ c) n$ K) e. |# I4 j! @" Ky2mean=1;
y3=1;
/ ]2 y+ T! l  N$ b2 _theta=0.9;
s=0.01;
4 m# v6 `& ~2 R" V9 [d1=0.1;
! L  J8 k; Y; j7 |- u; ?& Gct1=s*(1-theta)*y1+d1/(1+r);% ct1也与s,d1相关；

6 K6 E; w* ]1 o5 W4 Z  @
5 j0 V6 e- ]/ a' |0 z% {$ Sf(1)=wt/ct1-beta*(1+r)*expect1;
f(2)=wt/ct1*y1*(1-theta)-beta*expect2;
9 _2 A2 C) P/ r5 Aend

pTDbn25

帮你顶一下

zzz.dan

) Z  y) B+ T/ H/ \2 l1 k你这个写的啥，你foc的函数的变量f1和f2里面x的变量呢，你这个方程写的有很大问题

nkkopd