fminbnd循环求和

replace

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

clc;
clear all;
syms t x;
M=[2 4 8 16 32 64 128 256];
9 V+ I; c9 q1 L3 F. I7 Iv=10.5;
8 M/ x  V: p' L9 H, ~sigma=10^(-3);
z = sym('z');
/ b& U' d( ~7 d7 s) V0 Lx = sym('x');
  O, b- m" U5 s3 T! mP0=10^(-3); %目标误码率
t=exp(-1/sigma^2);
r=20; %平均SNR
R=10.^(r/10);
7 y9 T% W: x% n* h% Dw=0;
( K" Z$ Y9 A4 ufor i=1:7
  n( ], x9 h2 t0 r4 ~+ ]
        h(1)=sqrt(1/(2.*R))*qfuncinv(P0);
* N( x4 p$ F: P6 p( u6 c- t9 S, x1 T        h(i+1)=(1/sin(pi/M(i+1)))*sqrt(1/(2.*R))*qfuncinv(((log2(M(i+1)))*P0)/2);
3 @* h) r: |) l& [% 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))
        w=w+f(z)
! p4 m, w5 U; A. uend
2 |7 n6 u3 P. r5 L" Q8 p8 ^6 A& n
[z,feval] = fminbnd( w, 0, 1 )

someone

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

I_believe

解决如下，可以运行，不过结果看上去有些问题，去查查你的目标函数以及参数是不是都写对了
9 f3 u2 }3 z% @% b6 y8 K
clc;
2 R* v0 p: `9 b* t3 v: E' c) Sclear all;
' c8 z! u5 S' nsyms t x;
M=[2 4 8 16 32 64 128 256];
v=10.5;
  c$ ^( o' W- |4 z6 d7 {) Hsigma=10^(-3);
' s$ `6 r3 v0 z5 n$ p6 U' z% z = sym('z');
$ \, a' B3 j+ z% x = sym('x');
0 n3 S% }. X) C8 |0 N, m, R3 c# }; NP0=10^(-3); %目标误码率
3 d* O' z- i+ ]5 ]2 x5 C5 at=exp(-1/sigma^2);
3 l, d, U% }0 s  wr=20; %平均SNR
3 C7 g% v& S. p7 o5 \. ^R=10.^(r/10);
4 }! \! I" C$ P( d4 ?6 s, g7 V! Vw=0;
for i=1:7

* B2 M" C1 v0 _        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);
        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))
% j; p$ {/ U  y) |) q        if i==1
& y6 A- T8 H: C3 R& D9 b( U$ I% _7 M            w=@(z) f(z);
2 ~! ?$ \$ X0 m        else
9 q* w/ b" r2 T4 ^9 u            w=@(z) w(z)+f(z);
        end
$ J( Z" q2 ]0 T3 z# h        
end

[z,feval] = fminbnd( w, 0, 1 );

R_myself

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