fminbnd循环求和

replace

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

clc;
clear all;
9 z  n* n" G- w; Ksyms t x;
M=[2 4 8 16 32 64 128 256];
v=10.5;
9 w( g8 I) k+ W4 Nsigma=10^(-3);
z = sym('z');
x = sym('x');
P0=10^(-3); %目标误码率
  r- L- J7 G: ?  Zt=exp(-1/sigma^2);
r=20; %平均SNR
4 M& M5 }8 e( e7 s# {) FR=10.^(r/10);
0 L; {, J; E4 Ew=0;
; |5 s9 Q% X2 F9 i& wfor i=1:7

  ~7 |4 V" y8 X- i$ O1 c9 i2 V        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))
+ F) M1 Y$ z/ ]0 ]- h- l' d        w=w+f(z)
end
9 V( ^1 R9 z; L& @  D2 H" `; l
[z,feval] = fminbnd( w, 0, 1 )

R_myself

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

I_believe

解决如下，可以运行，不过结果看上去有些问题，去查查你的目标函数以及参数是不是都写对了
* U# S+ K- p: N6 {  z
9 Y7 O7 C# S" E$ j  e2 Kclc;
1 A9 @1 T2 F) H! L. R, K& uclear all;
syms t x;
M=[2 4 8 16 32 64 128 256];
v=10.5;
sigma=10^(-3);
% z = sym('z');
% x = sym('x');
P0=10^(-3); %目标误码率
& `5 C% j: T% D& v( y- rt=exp(-1/sigma^2);
4 o! Q* o7 L% L3 F; I. Vr=20; %平均SNR
$ k4 |0 U$ p4 {& w  hR=10.^(r/10);
$ U# t9 Z# i1 L: m% y1 ]( xw=0;
# E4 O$ `3 \8 ^: k! l& {for i=1:7

        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))
$ K  a, Z/ e$ m# r$ H+ w        if i==1
            w=@(z) f(z);
        else
5 h( Z1 h) M7 y8 U, B3 {1 |            w=@(z) w(z)+f(z);
        end
        
4 U, n& }1 T" m0 N* Z; kend

9 y$ _5 s8 m0 L! z# w0 F[z,feval] = fminbnd( w, 0, 1 );

someone

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