fminbnd循环求和

replace

代码一直在报错，不知道问题在哪里

clc;
clear all;
syms t x;
" d0 p! q2 i, B, i8 n% iM=[2 4 8 16 32 64 128 256];
v=10.5;
; f" S0 j: P: q6 \" }6 psigma=10^(-3);
z = sym('z');
  B! q$ Y0 j0 C4 a1 M" kx = sym('x');
P0=10^(-3); %目标误码率
t=exp(-1/sigma^2);
r=20; %平均SNR
, z! E# b  b) J" CR=10.^(r/10);
w=0;
for i=1:7

) k& E) Q$ P6 K8 q: M! ~% m        h(1)=sqrt(1/(2.*R))*qfuncinv(P0);
        h(i+1)=(1/sin(pi/M(i+1)))*sqrt(1/(2.*R))*qfuncinv(((log2(M(i+1)))*P0)/2);
; T% J$ S9 M1 B$ ~' V/ i        f=@(z) integral(@(x) 0.5.*eRFc(0.5.*x.*v.*z.*t.^(z.^2.*x.^2).*sin(pi./M(i)).*(x./sigma.^2).*(t.^(-x.^2./2))),h(i),h(i+1))
        w=w+f(z)
  L( L# w3 U! ~0 n* M* Z. fend
/ O* U4 Y2 K1 k( p- Z3 w+ j6 ^% T! }
8 g+ |0 x6 J# U5 \, Q9 Q[z,feval] = fminbnd( w, 0, 1 )

someone

仅供参考，匿名函数f在定义之后，调用时f(z)里的z要是具体数值，否则无法计算匿名函数值

I_believe

解决如下，可以运行，不过结果看上去有些问题，去查查你的目标函数以及参数是不是都写对了

, o; S2 `' E3 U& eclc;
1 P& A0 j  m! e5 p: dclear all;
syms t x;
0 r# U" N& X( O/ @6 J3 fM=[2 4 8 16 32 64 128 256];
& J4 S6 ?5 }" q7 h0 W  Z; `( v( rv=10.5;
+ y- w8 |  m+ R& Q' F$ o; {' L, X8 usigma=10^(-3);
% z = sym('z');
% x = sym('x');
9 X/ @$ m2 A% f; D) ]1 P' gP0=10^(-3); %目标误码率
t=exp(-1/sigma^2);
r=20; %平均SNR
9 C0 g0 ?) x4 P6 M4 m8 JR=10.^(r/10);
w=0;
' `" y8 C6 F- J0 Tfor i=1:7
+ Q2 w, F1 V6 q0 O
        h(1)=sqrt(1/(2.*R))*qfuncinv(P0);
5 b7 |! N3 ?: Q% {8 b3 X+ E        h(i+1)=(1/sin(pi/M(i+1)))*sqrt(1/(2.*R))*qfuncinv(((log2(M(i+1)))*P0)/2);
        f=@(z) integral(@(x) 0.5.*erfc(0.5.*x.*v.*z.*t.^(z.^2.*x.^2).*sin(pi./M(i)).*(x./sigma.^2).*(t.^(-x.^2./2))),h(i),h(i+1))
. [) h* Y! ]2 g1 _1 V9 P0 M" e        if i==1
* f3 r) u# L6 d0 g            w=@(z) f(z);
, @7 _) b+ |: e) o0 }* t  G        else
& P" W9 I& S* U: Q            w=@(z) w(z)+f(z);
        end
7 @' f) A1 {3 d- S5 }        
( v0 V8 a( j+ N7 A; Y, \1 D1 N2 V3 gend
" ?7 l: R$ ]! g  L3 A7 \3 Q7 M- _
[z,feval] = fminbnd( w, 0, 1 );

R_myself

检查一下目标函数及参数设置