fminbnd循环求和

replace

代码一直在报错，不知道问题在哪里
5 f/ k- |) ^1 V! ]& k( f7 E. R
clc;
clear all;
syms t x;
+ ?8 o' T0 M. q: ^: L  f9 |6 E7 e) ZM=[2 4 8 16 32 64 128 256];
8 M& K2 \: }* Z. i' r5 J/ J# N# a+ @v=10.5;
sigma=10^(-3);
' I8 n' M- A1 w8 I6 v/ Qz = sym('z');
2 y9 L5 O$ Y1 r& L! s& W! @; Rx = sym('x');
P0=10^(-3); %目标误码率
t=exp(-1/sigma^2);
r=20; %平均SNR
$ N+ Q; U6 m3 n) }R=10.^(r/10);
7 n+ Y; [6 z% u0 nw=0;
for i=1:7
: M- t7 z7 v3 b( L
% ]2 k+ J8 z# i        h(1)=sqrt(1/(2.*R))*qfuncinv(P0);
0 Y5 B" ?1 E- y) v) q" t        h(i+1)=(1/sin(pi/M(i+1)))*sqrt(1/(2.*R))*qfuncinv(((log2(M(i+1)))*P0)/2);
" ]' k9 Q$ ^; N0 P% X" o        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))
# u+ J$ z2 q) E- r1 V2 ]1 y* [        w=w+f(z)
$ t$ W3 G" v2 m( p" e) j0 jend

# U# P# q4 [6 c& K( \$ ]& U, t: B[z,feval] = fminbnd( w, 0, 1 )
$ Q0 W- K# A% M  ?

someone

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

I_believe

解决如下，可以运行，不过结果看上去有些问题，去查查你的目标函数以及参数是不是都写对了
3 j3 k4 s. p# O1 _' K7 ?
clc;
clear all;
syms t x;
M=[2 4 8 16 32 64 128 256];
! J8 S5 V, U+ q# N! }( Tv=10.5;
. @3 M* U' ~0 R, psigma=10^(-3);
/ f& F1 ^) `) C. Q% z = sym('z');
% x = sym('x');
6 D! [# o+ C) I: ^5 `" P$ S9 ~8 KP0=10^(-3); %目标误码率
t=exp(-1/sigma^2);
) P1 f! Z  C) @4 z7 mr=20; %平均SNR
R=10.^(r/10);
# q: T, @4 a' s! j' q! w% J0 dw=0;
for i=1:7

        h(1)=sqrt(1/(2.*R))*qfuncinv(P0);
% b" j# f" ^) w  I        h(i+1)=(1/sin(pi/M(i+1)))*sqrt(1/(2.*R))*qfuncinv(((log2(M(i+1)))*P0)/2);
; N$ \2 r( ^3 o6 c% \! q        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))
) }! u8 v2 B# u        if i==1
            w=@(z) f(z);
* \3 {) c; U/ g; O        else
( T. J& E' b! P* Y, h, o7 z            w=@(z) w(z)+f(z);
6 J/ G6 S' j& ^2 z5 y  e* q* p8 g        end
        
2 g5 S  `, L. O: V; T/ L* h  rend
8 @9 ]. N) ^' w8 @
[z,feval] = fminbnd( w, 0, 1 );

R_myself

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