fminbnd循环求和

replace

代码一直在报错，不知道问题在哪里
9 w& e" L9 h/ ~2 q: a% `+ O; n" ?2 _
; h; [/ y: ~1 P6 [$ p3 H5 Kclc;
clear all;
syms t x;
M=[2 4 8 16 32 64 128 256];
) T" {: a' r  ?v=10.5;
sigma=10^(-3);
z = sym('z');
0 U3 [8 v/ q7 B" w( C8 J: [x = sym('x');
P0=10^(-3); %目标误码率
, ?5 i: V/ {! A4 i: ?2 tt=exp(-1/sigma^2);
) o! f. B4 g; T. |1 O% |, Vr=20; %平均SNR
8 L$ r; h- |  O& ~( m; w3 K7 XR=10.^(r/10);
w=0;
for i=1:7

- X' J; ?! ]  U% |0 A* D% t        h(1)=sqrt(1/(2.*R))*qfuncinv(P0);
0 M+ f# N3 \7 v  U+ X# |. N0 H' Z        h(i+1)=(1/sin(pi/M(i+1)))*sqrt(1/(2.*R))*qfuncinv(((log2(M(i+1)))*P0)/2);
6 o- Z0 X3 y8 v% t7 ?        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))
: |# ?& q  e, S: {) P        w=w+f(z)
" d" X/ R9 h: w; O, w2 yend
- J! x" l$ r) _
* o+ N9 s  S( R' s6 S1 A. [[z,feval] = fminbnd( w, 0, 1 )
, y+ P& O9 f- H+ b

someone

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

I_believe

解决如下，可以运行，不过结果看上去有些问题，去查查你的目标函数以及参数是不是都写对了
- Z7 @+ q, g1 Q0 Z$ Z* T1 v
& @5 h9 a- P) w  `clc;
( D9 _: y/ _! [1 \! Jclear all;
syms t x;
M=[2 4 8 16 32 64 128 256];
$ E2 n/ N* j2 Q& _+ hv=10.5;
# A: N/ R: u) ^1 jsigma=10^(-3);
# u5 w5 }; i$ t5 ]. c* c6 z% z = sym('z');
% x = sym('x');
$ F4 ^6 o) w, o- g! _P0=10^(-3); %目标误码率
% k6 |' F2 g2 P7 ^7 }$ d. A3 Dt=exp(-1/sigma^2);
r=20; %平均SNR
R=10.^(r/10);
( o& T. m# _, a9 n( Z9 Ow=0;
4 `. I% }/ ?+ v) qfor i=1:7
) z  r4 a) a4 r
& j3 P% S9 W, A! ^3 n, z        h(1)=sqrt(1/(2.*R))*qfuncinv(P0);
# ?% M* k4 V7 X8 _; x        h(i+1)=(1/sin(pi/M(i+1)))*sqrt(1/(2.*R))*qfuncinv(((log2(M(i+1)))*P0)/2);
" d3 A7 Y4 Q$ q; x) K% y        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))
        if i==1
            w=@(z) f(z);
" Z: R* @/ A, S; M+ {9 t" E        else
' q& g) B; ~# g' r7 F, j5 O1 Z9 f            w=@(z) w(z)+f(z);
) i; K* r' n5 D, h+ j, J4 a4 X8 _        end
        
end
  l0 I7 Q$ }( E) R
[z,feval] = fminbnd( w, 0, 1 );

R_myself

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