matlab优化求解的问题

Uifhjvv

关于matlab优化求解的问题请教大家；程序报错截图如下，我分析报错的原因是迭代点（s,d1）改变时，（expect1，expect2）未改变，具体问题描述如下：
1.expect1，expect2是s,d1的函数，当fsolve函数按照优化算法找到新的点（s',d1'）时，我的问题里时需要expect1，expect2也变化，这样fsolve的目标函数就也变化了；
" d, u! \- a2 Y* Z( R2.根据上述1的描述，这个优化问题是否应该将（expect1，expect2）也纳入优化求解目标中去？这样该如何实现？
clc;
* x( |5 t& X0 ~clear all;
%%%  赋值外生给定的参数  %%%
wn=0.3;
wt=1-wn;
; f  j& T5 _4 C- p! E$ K1 D8 Sr=0.05;
beta=1/(1+r);
epislon=0.05;
phi=0.3;
0 x+ W: j% t* v8 E4 t: C+ Cy1=1;
# v$ q/ `& }) f9 r8 \yn2=1;
y2mean=1;
4 ~! P/ [! h# ^3 d& W6 N& Ay3=1;
theta=0.9;


% 赋值(y2,s,d1)
y2=0.95:0.01:1.05;
s=0.1; %给s,d1的任意初值
d1=0.1;
x0 =[s,d1];

& G7 L0 S. \1 l* u$ U$ w% V9 S
%在这个for循环里得到均值expect1,expect2，改变(s,d1),均值expect1,expect2也改变
7 X( K' K. J  L- |
& n# E7 s- A& B1 L( n5 Z3 h
for i=1:1:length(y2);
s=0.01;
d1=0.1;


%%% 判断是否紧约束 %%%
6 a  ]1 |# ?$ ystandard(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);
  }# R9 u  v4 I* r* F0 oif check(i)>standard(i)
    ct2(i)=((1-s)*y2(i)-d1+phi*(1-s)*y2(i))/(1-phi*wn/wt);
    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));
   
else
7 _* y! N. @8 x# |, I    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);
    ct3(i)=(1-s)*y3-d2(i);

% L! D: \2 a" B( x1 x  W9 V
    mu(i)=0;
  f% W6 Z3 R* H6 p  I; Q( v2 E3 zend
e1(i)=wt/ct2(i);%每给定一个y2(i)，计算得到一个e1(i)，e2(i)
, Y" c0 N: g. D4 g( Ye2(i)=wt/ct2(i)*y2(i)+mu(i)*phi*y2(i)+beta*wt/ct3(i)*y3;


%%%求均值%%%
expect1=mean(e1(1,：));
! t# N3 T8 q5 e( K$ _expect2=mean(e2(1,：));
# F5 ]8 b/ ^0 eend

# n, t" p* m6 J
: j9 R% `" n! e9 d5 u7 j%我的问题中，expect1，expect2实质上是s,d1的函数
) h1 s; B) ?+ F1 O. |$ N, H%在用fsolve优化求解的时候，每迭代一个新的（s,d1),都应该运行一次上述for这个循环，得到新的expect1、expect2再进行优化求解，
%也即（s,d1)变化时，expect1、expect2也要变化，如此fsolve的目标函数也变化
%这样一个问题没有编程思路，请求大家给予指导，如果我的问题没有表述清楚，十分欢迎追问，谢谢。
* L3 J0 b) l! H
$ L' o' W- D0 T  v& L
[x,fval,exitflag,output,]=fsolve(@(x)foc(x,expect1,expect2),x0);
4 v. T- ]( X. [5 X

: m$ q; x% V* \: ]: U' Nfunction f=foc(x,expect1,expect2)
: j% Y( ~  O" Nwn=0.3;
( ~% q7 W9 o2 G  j/ q. V9 ^0 Lwt=1-wn;
7 f* r6 [/ ]- i" B1 S5 i8 }r=0.05;
beta=1/(1+r);
epislon=0.05;
phi=0.3;
) e8 O: L1 c8 K% g8 Hy1=1;
yn2=1;
/ [  J4 h8 K: j/ ny2mean=1;
  ^0 I( A, B* u2 J6 H" g* xy3=1;
0 s* F$ O5 I' P3 }& d  Btheta=0.9;
# G+ A. w( H; s/ M& Ys=0.01;
d1=0.1;
5 i6 I+ v* t, q7 Bct1=s*(1-theta)*y1+d1/(1+r);% ct1也与s,d1相关；

: p9 G; m# l  ?3 T) H  M# b3 P
f(1)=wt/ct1-beta*(1+r)*expect1;
f(2)=wt/ct1*y1*(1-theta)-beta*expect2;
" f7 ^2 v9 d- G: r) e1 B% t' bend

1 @$ s0 M' @1 r" @- a1 X
' `" P% B3 i3 G

pTDbn25

帮你顶一下

zzz.dan

' {7 X9 Y) J8 X4 U( N你这个写的啥，你foc的函数的变量f1和f2里面x的变量呢，你这个方程写的有很大问题

nkkopd